Mon. Not. R. Astron. Soc. 000, ??-?? (2012) Printed 11 September 2012 (MN JATfcjX style file v2.2) 



A Suzaku survey of Fe K lines in Seyfert 1 AGN 



A.R. Patrick 1 , J.N. Reeves 1 ' 2 , D. Porquet 3 , A.G. Markowitz 4 ' 5 ' 6 , V. Braito 7 ' 8 , A.P. Lobban 1 ' 7 

1 Astrophysics Group, School of Physical Sciences, Keele University, Keele, Staffordshire, ST5 5BG, UK 

2 Department of Physics, University of Maryland Baltimore County, MD 21250, USA 

3 Observatoire astronomique de Strasbourg, Universite de Strasbourg, CNRS, UMR 7550, 11 rue de VUniversite, F-67000 Strasbourg, France 
4 Dr. Karl Remeis-Sternwarte and Erlangen Centre for Astroparticle Physics, Frederic-Alexander Universitdt Erlangen-Nurnberg, 

7 Sternwartstrasse, 96049 Bamberg, Germany 
. . 5 Alexander von Humboldt Fellow 

£Nj , 6 Center for Astrophysics and Space Sciences, University of California, San Diego, M.C 0424, La Jolla, CA, 92093-0424, USA 
T— H 1 7 X-Ray Astronomy Observational Group, Department of Physics and Astronomy, Leicester University, Leicester LEI 7RH 
Z ~_~) • S INAF - Osservatorio Astronomico di Brera, via E. Bianchi 46, 23807, Merate (LC) 



11 September 2012 

ABSTRACT 

We construct full broadband models in an analysis of Suzaku observations of nearby 
Seyfert 1 AGN (z < 0.2) with exposures > 50 ks and with greater than 30000 counts 
in order to study their iron line profiles. This results in a sample of 46 objects and 
1 84 observations. After a full modelling of the broadband Suzaku and SVift-BAT data 

CU (0.6-100 keV) we find complex warm absorption is present in 59% of the objects in 

this sample which has a significant bearing upon the derived Fe K region parameters. 
Meanwhile 35% of the 46 objects require some degree of high column density partial 
coverer in order to fully model the hard X-ray spectrum. We also find that a large 
number of the objects in the sample (22%) require high velocity, high ionization out- 
flows in the FeK region resulting from FeXXV and FeXXVI. A further four AGN 
feature highly ionized Fe K absorbers consistent with zero outflow velocity, making a 
total of 14/46 (30%) AGN in this sample showing evidence for statistically significant 
absorption in the Fe K region. 

Narrow FeKa emission from distant material at 6.4 keV is found to be almost 
ubiquitous in these AGN. Examining the 6-7 keV FeK region we note that narrow 
emission lines originating from FeXXV at 6.63-6.70 keV and from FeXXVI at 6.97 keV 
are present in 52% and 39% of objects respectively. 

Our results suggest statistically significant relativistic Fe Ka emission is detected 
in 23 of 46 objects (50%) at > 99.5% confidence, measuring an average emissivity index 
of q = 2.4 ± 0.1 and equivalent width EW = 96 ± 10 eV using the RELLINE model. 
When parameterised with a Gaussian we find an average line energy of 6.32±0.04kcV, 
Owidth = 0.470±0.05keV and EW = 97± 19 eV. Where we can place constraints upon 
the black hole spin parameter a, we do not require a maximally spinning black hole 
in all cases. 
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1 INTRODUCTION 

The analysis of the X-ray spectra of AGN can reveal infor- 
mation regarding the inner regions of the accretion disc, the 
AGN environment as a whole and subsequently the super 
massive black hole (SMBH) at its heart. It was suggested 
by Fabian et al. (1989) that emission occurring from the 
very inner regions of the accretion disc may be visible and 
subsequently broadened by Doppler motions and relativistic 
affects. The majority of AGN spectra show narrow line emis- 
sion from neutral iron at 6.4 keV (FeKa) likely originating 



from distant material e.g. the torus or the outer regions of 
the accretion disc (Krolik & Kallman 1987; Nandra 2006), 
particularly strong due to the high abundance and fluores- 
cent yield of iron. If Fe Ka emission additionally arises from 
material close to the central SMBH it will become relativis- 
tically broadened (Fabian et al. 1989; Laor 1991), producing 
a both blue and red-wings to the traditionally narrow line 
profile. 

In some AGN spectra this relativistic FeKa emission 
may be strong enough to be observed allowing its shape and 
strength measured using disk-line emission models such as 
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LAOR, kyrline, kerrdisk and relline (Laor 1991; Dovciak 
et al. 2004; Brenneman & Reynolds 2006; Dauser et al. 
2010). The application of these models allows properties 
such as the inclination and emissivity index of the disc to be 
measured in addition to the typical inner radius of emission 
and in some cases the spin of the central SMBH (Nandra 
et al. 2007; Patrick et al. 2011a). Gaining information re- 
garding the distribution of SMBH spins is an essential tool 
in aiding our understanding of galaxy evolution and dis- 
tinguishing between models such as prolonged or chaotic 
accretion and also the effect of mergers upon the SMBH 
spin (Hughes & Blandford 2003; Volonteri et al. 2005; King 
& Pringle 2007; Rezzolla et al. 2008). A spin distribution 
skewed towards higher values (a ~ 0.998) would suggest 
prolonged accretion, whereas low SMBH spin (a ~ 0) would 
indicate chaotic accretion models are more appropriate. In 
addition to this, the magnetic extraction of BH rotational 
energy through the Blandford-Znajek effect (Blandford & 
Znajek 1977) could cause a reduction in the spin (i.e. to- 
wards zero) of the SMBH in some AGN (Berti & Volonteri 
2008). 

Recent publications have made steps towards making 
spin estimates of SMBHs in a variety of AGN, including 
those which feature low levels of intrinsic absorption thereby 
offering the simplest spectrum to analyse, avoiding com- 
plications regarding the degree of spectral curvature intro- 
duced with warm absorbing zones (Miniutti et al. 2009; 
Schmoll et al. 2009; Patrick et al. 2011a; Emmanoulopoulos 
et al. 2011). More complex AGN spectra have also been anal- 
ysed and revealed further spin estimates e.g. MCG-06-30-15 
(Miniutti et al. 2007; Patrick et al. 2011b), Mrk 79 (Gallo et 
al. 2011) and NGC 3783 (Brenneman et al. 2011; Patrick et 
al. 2011b). However, as discussed in Patrick et al. (2011a), 
the estimated SMBH spin is highly model dependent and 
strongly related to the treatment of features such as the soft 
excess or any intrinsic absorbing zones. Assuming a Comp- 
tonization origin of the soft excess results in a range of low 
to intermediate spins, whereas using a high degree of rela- 
tivistic blurring to smooth the discrete soft emission lines 
into a continuum typically forces the spin to near maximal 
values and requires very high disc emissivities. 

This paper includes a sample of AGN from the public 
Suzaku archive of all observations of Seyfert 1 AGN with to- 
tal exposures > 50 ks and more than 30 000 counts in order 
to increase the likelihood of detection and broadened emis- 
sion from the inner regions if it is present. Suzaku is the 
ideal instrument with which to do this work since it allows 
us to gather both soft and hard X-ray data simultaneously 
using the X-ray Imaging Spectrometer (Koyama et al. 2007) 
and Hard X-ray Detector (Takahashi et al. 2007) detectors 
which, when combined with further non-simultaneous hard 
X-ray data from Swift BAT, gives a broad energy band- 
pass of 0.6-100.0 keV. The crucial difference being that we 
can obtain data regarding the strength of the Compton re- 
flection hump at ~ 30keV (George & Fabian 1991) which is 
beyond the capabilities of other current X-ray observatories. 
Only with hard X-ray data can the strength of the reflec- 
tion component be appropriately constrained and hence its 
contribution to the Fe K region assessed prior to attempting 
to determine broadening in the FeK region and eventually 
estimates upon SMBH spin. 

This is the final paper in a series of three in which a 



methodical and relatively uniform approach has been taken 
in an attempt to constrain accretion disc properties and 
SMBH spin from the FeK regions from an analysis of the 
X-ray spectra of Seyfert 1 AGN. In Patrick et al. (2011a) a 
small sample of six 'bare' Seyfert 1 AGN (i.e. those featur- 
ing low intrinsic absorption: Ark 120, Fairall 9, MCG-02-14- 
009, Mrk 335, NGC 7469 and SWIFT J2127.4+5654) was 
analysed, finding that narrow ionized emission lines such as 
FeXXVI are relatively common (4/6 objects) while the emis- 
sivity index of the accretion disc indicated that no strongly 
centrally concentrated emission was required to model any 
relativistic broadening in the Fe K region with an average of 
g~ 2.3 ±0.2. 

Patrick et al. (2011b) undertook and analysis of high 
quality, long exposure (> 200 ks) observations of Seyfert 1 
AGN with Suzaku (Fairall 9, MCG-06-30-15, NGC 3516, 
NGC 3783 and NGC 4051) making use of the full 0.6- 
100.0 keV bandpass in order to fully account for any warm 
absorber component to give the best possible opportunity 
to make estimates upon SMBH spin, also finding a low to 
moderate average emissivity index of q ~ 2.8±0.2. Other au- 
thors have also made some progress towards making spin es- 
timates on Seyfert 1 AGN e.g. Miniutti et al. (2009), Schmoll 
et al. (2009), Nardini et al. (2011), Gallo et al. (2011), Em- 
manoulopoulos et al. (2011) and Brenneman et al. (2011). 

The main aim of this paper is to assess the proper- 
ties and total percentage of AGN which have been observed 
with Suzaku that show evidence for broadened line emis- 
sion in the FeK region from the inner regions of the accre- 
tion disc resulting from an analysis of the broad-band X-ray 
spectrum. In this paper we expand our broad-band spec- 
tral analysis with Suzaku to all the currently archived type 
1-1.9 AGN, which have at least a 50 ks total exposure and 
30000 XIS band counts, sufficient for a broad band spectral 
analysis. This enables us to measure the iron line and re- 
flection properties of a more substantial sample of 46 type 
1 AGN, allowing the overall properties of the iron line and 
accretion disc to be investigated. We also aim to investigate 
Suzaku 's view of ionised emission and absorption lines in the 
FeK region and the occurrence or warm absorbers, highly 
ionised absorbers and partially covering absorbers and the 
subsequent effects upon the AGN X-ray spectrum. 



2 OBSERVATIONS & DATA REDUCTION 

2.1 Observations and sample selection 

The objects included within this sample are listed in Table 
[T]and are all the Seyfert 1-1.9 AGN with exposures > 50 ks 
and greater than 30 000 0.6-10.0 keV counts which have been 
observed with Suzaku with data publically available in the 
Suzaku data archiveQ as of 09/11. We also include data from 
some type 1 radio-loud (BLRGs - non-Blazar) AGN pro- 
vided they fit the above exposure and count criteria. High 
energy X-ray data from Swift /B AT from the 58 month BAT 
catalogue is also used in addition to that obtained from 
the HXD detector onboard Suzaku (but allowing the rela- 
tive cross-normalisation to vary), therefore the total energy 
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Table 1. The Suzaku Seyfert sample 



Object 


RA 


(J2000) 


Dec 


(J2000) 


Redshift 


N H (Gal) (10 22 


1H 0419-577 


04 


26 


00. 


.8 


-57 


12 


00.4 


o 


.1040 


0.0126 


3C 111 


04 


18 


21. 


.3 


+38 


01 


35.8 





.0485 


0.2910 


3C 120 


04 


33 


11. 


.1 


+05 


21 


15.6 


o 


.0330 


0.1060 


3C 382 


18 


35 


03. 


.4 


+32 


41 


46.8 


o 


.0579 


0.0714 


3C 390.3 


18 


42 


09 


.0 


+79 


46 


17.1 





.0561 


0.0347 


3C 445 


22 


23 


49. 


.5 


-02 


06 


12.9 





.0559 


0.0559 


4C 74.26 


20 


42 


37. 


.3 


+75 


08 


02.4 





.1040 


0.1160 


Ark 120 


05 


16 


11. 


.4 


-00 


08 


59.4 





.0327 


0.0978 


Ark 564 


22 


42 


39. 


.3 


+29 


43 


31.3 





.0247 


0.0534 


Fairall 9 


01 


23 


45. 


.8 


-58 


48 


20.5 





.0470 


0.0316 


IC 4329A 


13 


49 


19. 


.3 


-30 


18 


34.0 





.0161 


0.0461 


IRAS 13224-3809 


13 


25 


19. 


A 


-38 


21 


52.7 





.0658 


0.0534 


MCG-02-14-009 


05 


16 


21. 


2 


10 


33 


41.4 





.0285 


0.0924 


MCG-02-58-22 


23 i 


04 


43.i 


55 


-08 


11 


08.6 





.0649 


0.0291 


MCG-05-23-16 


09 


47 


40. 


.2 


-30 


56 


55.9 





.0085 


0.0870 


MCG-06-30-15 


13 


35 


53. 


.8 


-34 


17 


44.1 





.0077 


0.0392 


MCG+8-11-11 


05 


54 


53. 


.6 


+46 


26 


21.6 





.0205 


0.1840 


MR 2251-178 


22 


54 


05. 


.8 


-17 


34 


55.0 





.0640 


0.0640 


Mrk 79 


07 


42 


32. 


.8 


+49 


48 


34.8 





.0222 


0.0527 


Mrk 110 


09 


25 


12. 


.9 


+52 


17 


10.5 


o 


.0353 


0.0130 


Mrk 205 


12 


21 


•14. 


.0 


+75 


18 


38.5 





.0708 


0.0280 


Mrk 279 


13 


53 


03 


5 


+69 


18 


29.6 





.0305 


0.0152 


Mrk 335 


00 


06 


19. 


5 


+20 


12 


10.5 





.0258 


0.0356 


Mrk 359 


01 


27 


32. 


.5 


+19 


10 


43.8 





.0174 


0.0426 


Mrk 509 


20 


■14 


09 


.7 


-10 


43 


24.5 





.0344 


0.0344 


Mrk 766 


12 


18 


26 


5 


+29 


48 


46.3 





.0129 


0.0178 


Mrk 841 


15 


01 


01. 


.2 


+10 


26 


16.2 





.0364 


0.0222 


NGC 1365 


03 


33 


36. 


A 


-36 


08 


25.5 





.0055 


0.0134 


NGC 2992 


09 


45 


42. 


A 


-14 


19 


35.0 





.0077 


0.0487 


NGC 3147 


10 


16 


53. 


7 


+73 


24 


02.7 





.0093 


0.0285 


NGC 3227 


10 


23 


30. 


.6 


+19 


51 


54.2 





.0039 


0.0199 


NGC 3516 


11 


06 


47. 


.5 


+72 


34 


06.9 





.0088 


0.0345 


NGC 3783 


11 


39 


01. 


.7 


-37 


■14 


18.9 





.0097 


0.0991 


NGC 4051 


12 


03 


09 


.6 


+44 


31 


52.8 





.0023 


0.0115 


NGC 4151 


12 


10 


32. 


.6 


+39 


24 


20.6 





.0033 


0.0230 


NGC 4593 


12 


39 


39. 


.4 


-05 


20 


39.3 





.0090 


0.0189 


NGC 5506 


14 


13 


1 1. 


.9 


03 


12 


27.3 





.0062 


0.0408 


NGC 5548 


14 


17 


59. 


5 


+25 


08 


12.4 


o 


.0172 


0.0155 


NGC 7213 


22 


09 


16. 


.3 


-47 


09 


59.8 





.0058 


0.0106 


NGC 7314 


22 


35 


46. 


.2 


26 


03 


01.7 





.0048 


0.0150 


NGC 7469 


23 


03 


15. 


.6 


+08 


52 


26.4 





.0163 


0.0445 


PDS 456 


17 


28 


19. 


.8 


-14 


15 


55.9 





.1840 


0.1960 


PG 1211+143 


12 


14 


17. 


7 


+14 


03 


12.6 





.0809 


0.0274 


RBS 1124 


12 


31 


36. 


.4 


+70 


14 


14.1 





.2080 


0.0152 


SWIFT J2127.4+5654 


21 


27 


45. 





+56 


56 


39.7 





.0144 


0.7650 


TON S180 


00 


57 


19. 


9 


-22 


22 


59.1 





.0620 


0.0136 



range covered is 0.6-100.0 keV. 
included are listed in Table iBll 



Details of the observations 



2.2 Data reduction 



In some instances, objects may have been observed on 
numerous occasions, provided that there is little variation 
between the data sets they are combined and a single anal- 
ysis is performed. If the observations do indeed vary a sep- 
arate analysis is performed on each data set, although with 
similar model components and inferred geometries where 
possible, e.g. the inclination angle of the accretion disc would 
be linked between observations. 



The Suzaku data in this paper were reduced using version 
6.10 of the HEASOFT data reduction and analysis package. 
The XIS source spectra were extracted using 3.0'circular 
regions centred on the source. Background spectra were 
also extracted using 3.0'circular regions, this time centred 
on a region of the CCD not featuring any of the source 
or Fe55 calibration regions. Only data from the front- 
illuminated XIS 0, 2 and 3 detectors were used due to their 
greater sensitivity at FeK energies, however data from the 
back-illuminated XIS 1 remains consistent with the front- 
illuminated detectors. Observations since November 2006, 
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however, do not include data from XIS 2 as it is no longer 
operational. 

Objects which have been observed on multiple occasions 
but where the X-ray spectra (or spectral shape) does not 
appear have varied by a significant amount (e.g. other than 
simple changes in power law normalisation) are combined 
into a single data file and then a suitable model is fit to the 
data. This is the case for IC 4329A, MCG-06-30-15, Mrk 
509, Mrk 841, NGC 2992 and NGC 5548. However, if there 
are discernible differences between multiple observations of 
the same object (e.g. differing spectra shapes or changes in 
absorption), these data sets are analysed separately while 
retaining as many similar model components as possible. 
Such an analysis is performed for Fairall 9, NGC 1365, NGC 
3227, NGC 3516, NGC 3783 and NGC 4051. Both of the 
2006 observations of NGC 5506 are combined into a single 
spectrum while the 2007 observation is kept separate. Sim- 
ilarly, the later three observations of 3C 120 are combined 
whereas the first observation (OBSID: 700001010, see Table 
IB 1 [) is analysed separately. MCG-06-30-15 was observed on 
three occasions in January 2006, this paper makes use of all 
three observations and the time averaged spectrum is used 
in the main analysis. Analysis of the individual observations 
yields FeK line profiles which were consistent within errors 
and as such all three of the January 2006 observations have 
been coadded. Systematic errors are not included in datasets 
other than NGC 1365 in which systematic errors are set at 
2%. 

The HEASOFT tool xisrmfgen was used to generate 
the XIS redistribution matrix files (rmf) and the ancillary 
response files were generated using xissimarfgen. The data 
from each of the front-illuminated XIS detectors were then 
co-added using mathpha, ADDRMF and ADDARF in order 
to increase signal to noise. We ignore all XIS data below 
0.6 keV, above 10.0 keV and between 1.7-2.0 keV due to cal- 
ibration uncertainties of the detectors around the Si K edge. 
Some of the most recent observations with Suzaku (2011 on- 
wards) show evidence for contamination on the CCD of XIS 
which primarily affects the soft X-ray energies resulting 
in a divergence of XIS with XIS 3 below 1 keV. This is 
notable in the two observations of Fairall 9; both XIS and 
3 are consistent in the 2007 observation, however, the 2010 
observation displays a notable difference between XIS and 
3 (see Figure[T]). In this analysis, we do not co-add XIS and 
3 for the contaminated Fairall 9 observation, instead prefer- 
ring to include them as separate data sets although limiting 
XIS to > 1.0 keV, thus retaining some of the sensitivity at 
Fe K energies. Both detectors, however, remain consistent at 
FeK energies, see Figure [5] 

Spectra from the HXD were extracted from the cleaned 
HXD/PIN events files and subsequently corrected for in- 
strument dead time using the tool hxddtcor. The tuned 
HXD/PIN non X-ray background (NXB) events files were 
used for background subtraction (Fukazawa et al. 2009) and 
generated with 10 x the actual background rate in order to 
reduce photon noise, also with identical good time intervals 
(GTIs) as used in the source events. A simulated cosmic 
X-ray background was also produced for each observation 
using XSPEC v 12.6.0q with a spectral form identical to 
that used in Gruber et al. (1999), this was then added to 
the corrected NXB to form a single background file for each 
observation. The appropriate response and flat field files 




1 1.2 1.4 1.6 

Observed Energy (keV) 




1 1.2 1.4 

Observed Energy (keV) 

Figure 1. Ratio plots of 0.6-1.0 keV of XIS (red) and XIS 3 
(black) for the 2007 (upper panel) and 2010 (lower panel) ob- 
servations of Fairall 9. These show the contamination which has 
recently affected XIS at soft X-ray energies causing a diver- 
gence with XIS 3, note that this is not apparent in the first ob- 
servation. The soft X-ray spectrum was modelled with a basic 
(powerlaw+COMpTT)*wabs model. 




5 6 7 8 9 

Observed Energy (keV) 

Figure 2. Ratio plot of the 4-9 keV FeK regions of XIS (red) 
and XIS 3 (black) for the 2010 observation of Fairall 9. Indicating 
that while contamination has a significant effect at soft X-ray 
energies, the Fe K region remains entirely consistent between XIS 
and XIS 3. 
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were used from the Suzaku CALDB suitable for the respec- 
tive epochs of each observation. HXD data is used over the 
15.0-60.0 keV range along with the appropriate PLN/XIS- 
FI cross-normalisation according to the epoch and nominal 
pointing position^. 

While the BAT data from Swift are not necessarily si- 
multaneous with the Suzaku observations and is of relatively 
low statistical weight, we find throughout the paper that the 
BAT spectral indices are consistent with those obtained from 
the XIS and HXD/PIN data. That is, due to variability of 
the reflection fraction over time scales of years the spectral 
shape may vary and hence may not be equivalent to the 
spectral shape of the hard X-rays obtained with Suzaku. 
However, removing the BAT data and refitting still yields 
the original photon index Y as well as the other spectral 
parameters and there appears to be little evidence that the 
spectral shape > 15 keV varies over the timescales consid- 
ered here. It is therefore appropriate to consider that there is 
a simple change in the normalisation of the hard X-ray com- 
ponent (as accounted for by the cross-normalisation between 
XIS and BAT being allowed to vary within the models) and 
the inclusion of BAT data from Swift serves to improve our 
view of the X-ray spectral shape above lOkeV. The BAT 
data therefore improves the accuracy to which we can esti- 
mate the powerlaw continuum and the reflection component. 



3 ANALYSIS & RESULTS 

Spectral analysis and model fitting is performed from within 
XSPEC v 12.6.0q (Arnaud 1996), all models are modified 
by Galactic absorption which is accounted for by the WABS 
multiplicative model (Morrison & McCammon 1983). The 
respective Galactic column densities were obtained using 
the nh ftool for each source giving the weighted average 
Nil value of the LAB Survey of Galactic H I (Kalberla et al. 
2005), using abundances from Anders & Grevesse (1989). 
Data are fit over the full 0.6-100.0 keV range available, ex- 
cluding those regions affected by the uncertainties in the 
XIS calibration mentioned above. The x 2 minimisation tech- 
nique is used throughout, all errors are quoted at the 90% 
confidence level (A^ 2 = 2.71 for one interesting parame- 
ter) and include statistical and not instrumental systematic 
errors. Where the significance of components is quoted in 
terms of A^ 2 , the component in question has been removed 
from the model and the data refit to ensure that the order 
in which components are added to the model does not affect 
the quoted statistical significance. Throughout this paper, a 
positive Ax 2 corresponds to a poorer fit, whereas a negative 
Ax 2 corresponds to an improvement in the fit. 



3.1 Spectral models 

The data used in this paper are selected such that there are 
a sufficient number of counts for detailed spectral analysis 
to be conducted. The wide 0.6-100.0 keV energy bandpass 
as used in this sample means that a full treatment must 
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be given to all the main components of the X-ray spec- 
trum of these Seyfert 1 AGN, namely the soft excess, Comp- 
ton reflection, Galactic absorption, warm absorbers and any 
highly ionized outflows or relativistic broadening which may 
be present. In order to model the X-ray spectra in this sam- 
ple as consistently and uniformly as possible we construct 
models following the criteria in the following sections. 

Model components are added according the residuals in 
the data and the statistical requirement of such a compo- 
nent. For example, a model to account for the soft excess 
is considered to be required, provided that its inclusion ex- 
ceeds the 90% confidence level according to A\ 2 given the 
appropriate number of free parameters. For example, an xs- 
TAR grid typically has two: column density and ionisation 
parameter although some objects may require a third due 
to an outflow velocity. 



3.1.1 Modelling of the soft excess 

A number of AGN spectra show an excess over a simple 
powerlaw at lower X-ray energies, typically < 2.0 keV, this 
relatively common feature has been termed the 'soft excess'. 
The origin of the soft excess is as yet unknown, however, 
there are a number of different interpretations and methods 
with which to model it. Perhaps the most basic method of 
modelling the soft excess is as a black body to replicate di- 
rect thermal emission from successive annuli of the accretion 
disc (Malkan & Sargent 1982). The inferred constant tem- 
perature of the soft excess component (which scales with 
Mgy'' 4 ) is in disagreement with the typical accretion disc 
properties of a SMBH and is arguably more suitable for an 
intermediate mass BH (Gierlihski & Done 2004). 

A variation upon this concept is that the soft excess is 
produced by the Compton up-scattering of EUV seed pho- 
tons from the disc in a hot plasma or corona lying above the 
disc. For example, the COMPTT model (Titarchuk 1994) al- 
lows for a large variation in photon seed temperature while 
still producing a relatively constant output photon temper- 
ature as is required by many of the observed soft excesses in 
AGN. This method has been successful in modelling the soft 
excess of a large number of AGN, including a sample of PG 
quasars (Porquet et al. 2004) and in Mrk 509 by Mehdipour 
et al. (2011). This interpretation of the soft excess is also 
used throughout this paper, consistent with the method in 
Patrick et al. (2011a,b); similar results can also be obtained 
with a second steep power law component. 

The above methods assume a smooth shape to account 
for the soft excess, however alternatives have been sug- 
gested e.g. an atomic origin (Brenneman & Reynolds 2006). 
Soft X-ray emission lines emitted from regions close to the 
central SMBH are relativistically blurred and broadened, 
much in the same way as the often observed broad 6.4 keV 
FeKa emission line and its red-wing. A series of discrete 
emission line features are then relativistically blurred to 
such an extent that they merge to form a smooth contin- 
uum (Ross, Fabian & Ballantyne 2002). A number of AGN 
have been modelled in this fashion e.g. Ark 120, Fairall 9, 
Mrk 335 & RBS 1124 (Schmoll et al. 2009; Miniutti et al. 
2010; Nardini et al. 2011; Patrick et al. 2011a), although 
a significant amount of relativistic blurring is required to 
smooth the characteristically narrow features into a broad 
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continuum forces the accretion disc emissivity index to high 
values (q ?t 4.5) and SMBH spin towards near maximal 
(a ~ 0.998). Indeed Schmoll et al. (2009) note that ignoring 
the XIS data below 2keV (to avoid fitting the soft excess) 
relaxes the spin constraint to low values. These parameters 
are often at odds with those similarly derived from the asym- 
metric line profile from the Fe K region, for example Miniutti 
et al. (2010) find q = 4.1±f, | and a > 0.6 in RBS 1124 while 
the authors note that there is no evidence for a strong broad 
line in the Fe K region. However, given the assumption that 
the material responsible for both the soft emission lines and 
FeKa emission is located in the same region, it is logical to 
assume that estimated accretion disc and SMBH parameters 
should be consistent (Patrick et al. 2011a). 

Note that relativistic blurring of the reflection compo- 
nent at soft X-ray energies can still contribute towards the 
soft excess, however, the soft X-ray flux from relativistically 
blurred emission of the inner regions is unlikely to prove par- 
ticularly strong in comparison to the flux from the Comp- 
tonization component due to any residuals in the Fe K region 
driving the fit with the majority of the soft X-ray flux al- 
ready accounted for with the COMPTT model. Therefore, if 
the FeK line is to be used as a diagnostic for determining 
accretion disc parameters and constraining SMBH spin, the 
treatment of the soft excess has a significant role to play in 
ensuring that it is not driving the fit if blurred reflection 
models are used in the fitting process. As our aim is to pa- 
rameterise the properties of the disc and reflection based on 
the FeK line profile, independent of the nature of the soft 
excess, we retain this approach of the Comptonization origin 
of the soft excess here. 

3.1.2 Compton reflection and emission lines 

As the narrow iron Ka line is ubiquitous in virtually all 
these AGN, an unblurrcd REFLIONX (Ross & Fabian 2005) 
component, measured down to an ionisation parameter of 
£ = 1 and additional narrow fluorescent emission lines rep- 
resentative of reflection off distant material (e.g. the torus or 
outer regions of the accretion disc) are included in all AGN 
spectra showing evidence for a hard excess. Previous stud- 
ies of the iron line regions in Seyfert AGN by Bianchi et al. 
(2004), Nandraet al. (2007) and Patrick et al. (2011a,b) also 
suggest that ionized species of iron are relatively common at 
energies of 6.7 keV and 6.97 keV for FeXXV and FeXXVI 
respectively. Neutral narrow FeKa and the accompanying 
FeK/3 emission at 6.4 keV and 7.056 keV respectively are 
ubiquitous in AGN spectra (Nandra et al. 2007). 

As stated above, in this paper we use the reflionx 
model to account for the distant near neutral reflection con- 
tinuum with the input photon index T tied to that of the in- 
trinsic powerlaw. Although soft emission lines are included 
in reflionx, these are also added on an ad-hoc basis when 
required and as such are modelled using narrow Gaussians 
of fixed width a — 0.01 keV e.g. OVIII or NelX emission 
from distant photo-ionized gas e.g. the BLR or NLR. FeK/3 
emission at 7.056 keV with flux Fk/i = 0.13 x Fn a is not in- 
cluded self-consistently in reflionx which is therefore mod- 
elled using a narrow Gaussian of fixed width a — 0.01 keV 
and with flux fixed at the value obtained during an ini- 
tial parametrisation of the FeK region and the narrow 
FeKa flux when modelled with a Gaussian. Ionized lines 



in the FeK region (such as FeXXV at ~ 6.63 — 6.7 keV and 
FeXXVI at 6.97 keV) are accounted for with narrow Gaus- 
sians (<j = 0.01 keV) if and when required, although in some 
circumstances these can be indistinguishable from the blue- 
wing of a relativistically broadened Fe Ka line profile. 

While the iron abundance Zf c is left free to vary, in 
some AGN the strength of the narrow Fe Ka core may force 
Zf c to unfeasibly high values (e.g. Zf c > 3) while improp- 
erly modelling the hard X-ray reflection spectrum due to the 
greater number of counts in the 5-6 keV region. To avoid this 
scenario, we fix Zf c at Solar abundance and add an addi- 
tional narrow Gaussian of fixed width (a — 0.01 keV, to pre- 
vent interference with any broad component from the disc) 
at 6.4 keV to model the Fe Ka core while maintaining a good 
fit to the HXD and BAT data. Rather than simply being an 
ad- hoc solution, this could be representative of FeKa emis- 
sion additionally arising from Compton-thin matter such as 
the BLR or NLR as well as, for example scattering off a 
distant Compton-thick torus. 

3.1.3 Warm absorption 

The X-ray spectra of many AGN feature one or more zones 
of warm absorbing gas, while primarily affecting the spec- 
trum at soft X-ray energies, with higher column densities 
(e.g. > 10 22 cm~ 2 ) it can add subtle curvature above 2.5 keV. 
Previous studies (e.g. Miniutti et al. 2007; Nandra et al. 
2007) restrict their analysis to 2.5-10.0 keV to avoid compli- 
cations with the warm absorber, instead choosing to focus 
in upon the FeK region. However, as found in Reeves et al. 
(2004), Turner et al. (2005), Miller et al. (2009) and Patrick 
ct al. (2011b) effects of the warm absorber extend even to 
the FeK region, contributing significantly to the strength 
of the observed red-wing below 6.4 keV. For example, in 
MCG-06-30-15 Miniutti et al. (2007) model the time av- 
eraged January 2006 spectra without absorption, resulting 
in an apparent strong broad FeKa line and near maximal 
SMBH spin; while a full treatment of the warm absorber 
appears to reduce the strength of the broad component and 
subsequently BH spin to more intermediate values (Zycki et 
al. 2010; Patrick et al. 2011b). 

To model the soft X-ray warm absorber components in 
this paper, we use an xstar (Kallman et al. 2004) gener- 
ated grid illuminated by an X-ray photon index of Y — 2.0, 
in agreement with the mean values found in radio-quiet type 
I AGN (Scott et al. 2011). The abundances are fixed at Solar 
values (except Ni, which was set to zero), where xstar uses 
the values of Grevesse & Sauval (1998). The turbulent veloc- 
ity is set to 100 km s -1 . The grid is well suited to account- 
ing for a variety of absorption zones due to its wide range 
in column density (5 x 10 18 cm~ 2 < JV H < 5 x 10 24 cm~ 2 , 
in steps of ANu = 10 19 cm~ 2 ) and ionization parameter 
(0 < log£ < 5, with grid steps A(log£) = 0.5). An elec- 
tron density of n c = 10 10 cm~ 3 is assumed, although the 
data are largely insensitive to the density at CCD resolu- 
tion. The computed XSTAR spectra were generated over the 
energy range from 0.05 — 120 keV, so are suited to broad- 
band spectra, especially at hard X-ray energies. 

During the fitting process absorption zones are added 
as required, in some objects more than one zone may be 
statistically required and in exceptional circumstances a grid 
with a higher turbulent velocity may be required to model 
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highly ionized absorption in the FeK region (see below), 
as at low turbulences the iron K absorption can become 
easily saturated. Absorption zones are added until a good 
overall fit is found to the data and there are no clear residuals 
remaining in the soft X-ray spectrum. 

3.1.4 Highly ionized absorption 

Further to typical warm absorption, highly ionized absorb- 
ing zones of gas (where present) also play an important role 
in the FeK region and in the determination of the mea- 
sured line parameters and strength of the observed red- 
wing (Reeves et al. 2004; Turner et al. 2005). Absorption 
lines indicative of such zones are fairly common in X-ray 
spectra of AGN e.g. ls-2p resonance lines from FeXXV 
and/or FeXXVI at their rest-frame energies, i.e. 6.7 keV and 
6.97 keV respectively, although there is evidence for Wue- 
shifted absorption lines in many AGN (Braito et al. 2007; 
Tombesi et al. 2010a; Lobban et al. 2011). 

In the event that absorption features in the Fe K region 
are found, which are not adequately modelled by the above 
low turbulence grid, then these are accounted for using a 
model representative of a highly ionized absorption zone. 
We use an xstar generated grid with a turbulent velocity 
of lOOOkms -1 , with an input X-ray photon index of T = 
2.0 and using the input SED from Tombesi et al. (2011). 
This grid covers a range in column density and ionization 
parameter of 1 x 10 2o cm~ 2 < A H < 1 x 10 24 cm~ 2 (in steps 
A Ah = 6 x 10 19 cm" 2 ) and < log£ < 6 (with A(logf) = 
0.32) respectively. The electron density assumed is also n c = 
10 10 cm" 3 . 

3.1.5 Partial covering 

Some models used in the analysis of these AGN use par- 
tial covering geometries whereby a fraction of the observed 
X-rays are absorbed by a surrounding gas in the line-of- 
sight (in addition to typical fully-covering absorbers) as de- 
scribed in Section 3.1.3, while some fraction of the contin- 
uum 'leaks' or scatters through and is unattenuated by the 
partially covering material (e.g. Miller et al. 2008). The par- 
tial coverer used here takes the form of (powerlaw + xs- 
tar*powerlaw) with the photon index of both the power- 
law components tied and normalisations of both powerlaws 
free to vary. The parameters of the warm absorber (ion- 
ization and column density) are also allowed to vary. The 
column density of the partially covering medium can have 
a significant affect upon the spectrum. For example, high 
column partial coverers (Ah ~ 10 24 cm -2 ) predominantly 
affect the hard X-ray energies and can contribute towards 
a hard X-ray excess in cases where a hard excess remains 
(i.e. in addition to reflionx). Lower column density partial 
coverers may greatly affect spectral curvature at lower X-ray 
energies and in some cases can entirely remove any 'broad' 
residuals in the FeK region (Miller et al. 2009). While using 
a partial coverer in this way provides an alternate explana- 
tion for observed curvature or a red-wing in the Fe K region, 
we restrict the analysis to a maximum of one partially cov- 
ering zone in order to not arbitrarily model away spectral 
features present in the spectra. The absorption grid used 
is the same as used to model the warm absorber, with a 
turbulence of lOOkrns -1 . 



3.1.6 Relativistic line emission 

The majority of X-ray lines originate from material suffi- 
ciently distant to the SMBH (e.g. the ubiquitous narrow 
6.4 keV line) such that effects such as relativistic Doppler 
motions and gravitational redshift have a negligible effect 
upon the observed spectrum. If, however, emission comes 
from the very inner regions of the accretion disc it will of 
course be subject to these effects. As discussed above, this 
is one interpretation for the commonly observed red-wing in 
the Fe K regions of some AGN which may remain after a full 
modelling of the broad-band spectrum, taking into account 
warm absorbers which may be present. If an excess remains 
at ~ 5 — 6.4 keV, we initially model this using the relativistic 
line emission model relline (Dauser et al. 2010). This al- 
lows properties such as the emissivity index and inclination 
of the disc to be measured in addition to placing estimates 
upon the spin of the central SMBH (or the inner radius of 
emission) . A more comprehensive modelling of the inner re- 
gions of the accretion disc would include the accompanying 
blurred reflection spectrum. For example a convolution of 
the relline kernel (i.e. relconv, Dauser et al. 2010) with 
reflionx allows relativistic effects to be applied to both 
the hard and soft X-ray reflection spectrum as well as to 
the FeKa emission line at 6.4 keV, we investigate such an 
approach in Section 3.4. 

3.2 Baseline model 

The baseline model is intended to model the entire 0.6- 

100.0 keV spectrum, accounting for features such as the un- 
derlying continuum, soft excess, distant reflection and both 
warm and neutral absorption. With this model we aim to as- 
sess the remaining residuals in the Fe K region which may be 
attributed to relativistic line emission from the inner regions 
of the accretion disc, i.e. we use a combination of the mod- 
els and methods outlined in Section 3.1 to form a baseline 
model which does not include broadened line emission. The 
baseline model acts as the null hypothesis whereby the X-ray 
spectrum of these AGN can be described and adequately fit 
by emission or reflection from purely distant material. Fig- 
ure [3] therefore shows the Fe K regions of these AGN prior 
to any modelling of broad emission, distant emission lines 
or highly ionised absorption. 

3.2.1 Analysis of 'bare' Seyfert spectra 

As noted in Patrick et al. (2011b) the warm absorber be- 
low 2.5 keV has a significant affect upon the X-ray spectrum 
at higher energies, particularly important when attempting 
to test for the strength or indeed presence of a relativisti- 
cally broadened red-wing in the FeK region. Subsequently 
the broadband 0.6-100.0 keV analysis conducted here sug- 
gests that the majority of Seyfert 1 AGN in this sample 
show evidence for at least one warm absorber zone or addi- 
tional neutral absorption over the simple Galactic absorbing 
column. 

Only 11/46 objects in this sample can be consid- 
ered 'bare' in that no absorption (either neutral, warm 
or partially covering) whatsoever is required: 3C 390.3, 
Ark 120, Fairall 9, MCG-02-14-009, Mrk 110, Mrk 335, 
Mrk 359, NGC 3147, NGC 7213, NGC 7469 and RBS 
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Table 2. Model applicability to each AGN. Comparison of components and features use to model each AGN; COMPTT is used to model 
the soft excess, an unblurred REFLIONX is used to model the reflection component if statistically required. Only an increase in the 
neutral Galactic column is required. 2 No XSTAR grid is required in MCG-05-23-16, however additional neutral absorption is used to 
absorb only the intrinsic POWERLAW. Residuals in the FeK region are considered to be present if there is an improvement of Ax 2 > 6.3 
with the addition of a broad Gaussian (see Table 131 . 3 Only required in the 2006 observation of Mrk 766. 







1 V.IL.ll.M.l^ Llt_Jll 


1 til Lldl ^_.l_JVl_.llll£i 


V V dl 111 CLlJoUl UC1 


llltill C ctUoUl l) 11U11 


KriQiniialc at Pr> T\ 

1 V\^oll_l Lldlo CI L 1_ l_. Iv 


1H 0419-577 


/ 


/ 


/ 


/ 


X 


X 


3C 111 


X 


X 


X 


X 1 


/ 


/ 


3C 120 


/ 


/ 


X 


X 1 


X 


/ 


3C 382 


/ 


/ 


X 


/ 


X 


/ 


3C 390.3 


/ 


/ 


X 


X 


X 


/ 


3C 445 


X 


/ 


/ 


/ 


/ 


/ 


4C 74.26 


X 


/ 


X 


/ 


X 


/ 


Ark 120 


/ 


/ 


X 


X 


X 


/ 


Ark 564 


/ 


/ 


/ 


/ 


X 


X 


Fairall 9 


/ 


/ 


X 


X 


X 


/ 


IC 4329A 


X 


/ 


X 


/ 


X 


/ 


IRAS 13224-3809 


/ 


X 


/ 


/ 


X 


/ 


MCG-02-14-009 


X 


/ 


X 


X 


X 


X 


MCG-02-58-22 


/ 


/ 


X 


/ 


X 


X 


MCG-05-23-16 


/ 


/ 


X 


X 2 


X 


/ 


MCG-06-30-15 


/ 


/ 


/ 


/ 


/ 


/ 


MCG+8-11-11 


X 


/ 


X 


/ 


X 


/ 


MR 2251-178 


/ 


X 


X 


/ 


/ 


/ 


Mrk 79 


X 


X 


X 


/ 


X 


/ 


Mrk 110 


/ 


/ 


X 


X 


X 


X 


Mrk 205 


/ 


/ 


/ 


X 


X 


X 


Mrk 279 


X 


/ 


/ 


/ 


X 


X 


Mrk 335 


/ 


/ 


X 


X 


X 


/ 


Mrk 359 


/ 


/ 


X 


X 


X 


/ 


Mrk 509 


/ 


/ 


X 


/ 


X 


/ 


Mrk 766 


/ 


/ 


/ 


/ 


/ 


/ 3 


Mrk 841 


/ 


/ 


X 


/ 


X 


/ 


NGC 1365 


X 


/ 


/ 


/ 


/ 


/ 


NGC 2992 


/ 


/ 


X 


X 1 


X 


/ 


NGC 3147 


X 


/ 


X 


X 


X 


/ 


NGC 3227 


/ 


/ 


/ 


/ 


/ 


/ 


NGC 3516 


X 


/ 


/ 


/ 


/ 


/ 


NGC 3783 


/ 


/ 


X 


/ 


/ 


/ 


NGC 4051 


/ 


/ 


/ 


/ 


/ 


/ 


NGC 4151 


/ 


/ 


/ 


/ 


/ 


X 


NGC 4593 


X 


/ 


X 


/ 


X 


/ 


NGC 5506 


X 


/ 


/ 


/ 


X 


/ 


NGC 5548 


X 


/ 


X 


/ 


/ 


X 


NGC 7213 


/ 


/ 


X 


X 


X 


X 


NGC 7314 


/ 


X 


X 


X 1 


X 


/ 


NGC 7469 


/ 


/ 


X 


X 


X 


/ 


PDS 456 


/ 


X 


/ 


X 


/ 


X 


PG 1211+143 


/ 


X 


X 


/ 


/ 


X 


RBS 1124 


/ 


/ 


X 


X 


X 


X 


SWIFT J2127.4+5654 


X 


/ 


X 


X 1 


X 


/ 


TON S180 


/ 


/ 


/ 


/ 


X 


X 


Fraction 


31/46 


39/46 


16/46 


27/46 


14/46 


32/46 



1124. These AGN are therefore straight forward to 
model, without any complications due to curvature from 
warm absorbing components, representing the most fidu- 
cial of AGN X-ray spectra simply consisting of (power- 
law+compTT+reflionx)*wabs. Slightly more complex 
are those objects which are free from warm absorption, how- 
ever, requiring an additional neutral absorbing column over 
and above that from the standard Galactic absorption cal- 
culated by the nh ftool. This is accounted for by multiplying 



the typical 'bare' AGN model by a single zphabs fixed at 
the redshift of the object with column density free to vary, as 
used in 6/46 of AGN in the sample: 3C 111, 3C 120, MCG- 
05-23-16, NGC 2992, NGC 7314 and SWIFT J2127.4+5654. 
Note that the higher neutral absorption in 3C 111 is likely 
due to the presence of a giant molecular cloud in the line 
of sight, see Bania et al. (1991) and Rivers, Markowitz & 
Rothschild (2011b). MCG-05-23-16 does not require an ad- 
ditional totally covering zone of neutral absorption, instead 
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Figure 3. Ratio plots of the 4-9 keV residuals modelled with the baseline model minus REFLIONX, narrow lines and FeK absorbers, i.e. 
the Fe K region is left totally unmodelled with only warm absorption at soft X-ray energies being taken into account in addition to the 
intrinsic powerlaw and COMPTT (where required). The data has been refit after the removal of components. The vertical dashed lines 
represent 6.4 keV in the observed frames. Red data points indicate those from further observations. 



requiring a neutral zphabs geometry whereby only the pow- 
erlaw component is absorbed. Two further AGN (Mrk 205 
and PDS 456) require a partial coverer, but not a warm ab- 
sorber. Hence a total of 19/46 objects in this sample do not 
indicate the presence of a warm absorbing component. 



3.2.2 Absorbed Seyfert spectra - fully covering 

The remaining objects in the sample (27/46) feature some 
degree of complex warm absorption, which is modelled with 
successive zones of a multiplicative XSTAR grid (see Table 
[2| . The majority of these AGN are well modelled using one 
or more fully covering XSTAR grids with turbulent velocity 



Vturb = 100 km s 1 fixed at the redshift of the object (i.e. 
from Table[T]) as described in Section 3.1.3. 

Rather using the typical iwb = 100 kms -1 XSTAR grid, 
some AGN in this sample require the use of an XSTAR gener- 
ated grid with iwb = 1000 kms -1 i.e. the turbulent velocity 
exceeds the local thermal velocity of the absorbing ion, prov- 
ing to be the dominant factor for absorption line broadening, 
see Figure [4] This scenario is required when modelling the 
soft X-ray absorber in only 3/28 of objects featuring com- 
plex absorption: 3C 445, Ark 564 and IC 4329A. IC 4329A 
requires two low i^turb XSTAR grids in addition to the single 
high wturb XSTAR grid. The warm absorber properties are 
parameterised in Table [B2l 
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Table 3. Broad Gaussian parametrisation of FeK region in addition to the baseline model. a Flux given in units (10 — 5 phcm~ 2 s~ 1 ). b 
Parameters other than normalisation tied between multiple observations, all are tied in NGC 5506. — Indicates that no broad Gaussian 
can be fit i.e. no discernible reduction in x 2 whatsoever. 



Object 








Broad line 








FeXXV 


Fe XXVI 




LineE (keV) 


°"Broad (keV) 


SVK Br0 ad (eV) 


Flux a 


A 2 


emission 


emission 


1H 0419-577 


_ 






_ 






_ 


X 


X 


3C 111 


6.28±g"gg 




< 0.136 


19+ 23 


0.55 


+0.67 
-0.27 


-11 


/ 


/ 


3C 120 


6.38l°;°i 





i r> .+0.025 
iz4 -0.022 


76+J 


4.04 


+0.35 
-0.49 


-75 


/ 


X 


3C 382 


6.45t°^ 





1 Oft ' u 1 
• lao _().0|)0 


42+lf 


2.45 


+ 1 07 
-1.01 


-20 


X 


X 


3C 390.3 


6.49loT 





onA + . 256 




3.29 


+ 1.17 
-0.78 


-40 


X 


X 


3C 445 


6.04+°-; 2 

— 0.11 




< 0.436 


23+?o 

— Id 


1.07 


+0.58 
-0.46 


-7 


/ 


X 


4C 74.26 


6.12+°-°? 




< 0.186 




1.07 


+0.58 
-0.46 


-20 


/ 


X 


Ark 120 


6.36tn 8 9 





oon+O.llO 
ci/U -0.090 


105+.|> 


0.70 


+0.32 
-0.32 


-12 


X 


/ 


Ark 564 


fi 42 +0.09 
u — 08 




< 0.161 




0.34 


+0.25 
-0.22 


-3 


/ 


X 


Fairall 9 


c qi+0^4 
o.yj._Q 2i 





dud -0.173 


49±?° 


1.45 


+0.58 
-0.55 


-36 


/ 


/ 


IC 4329A 


6-34±g : ig 


0. 


cci+0.119 
ooi -0.100 




7.94 


+ 1.96 
-1.91 


-53 


X 


/ 


IRAS 13224-3809 


6.19+S-Sf 

^ ^ — (l.Oo 


1 


nl .+0.613 
ui4 -0.346 


632+29?, 

— 2o0 


0.37 


+0.18 
-0.15 


-9 


/ 


X 


MCO— 02-1 4-009 


6 33+ 021 
o.oo_ 24 





oQfi+0. 280 
zou -0.130 


09+59 

az -56 


0.49 


+0.31 
-0.30 


.3 


x 


/ 


MCG-02-58-22 
















X 


X 


MCG-05-23-16 


6-34+°'?n 

— U. 1U 





■°- LZ -{).()78 


108+1 7 




-15 


/ 


X 


MCG-06-30-15 


5.93+HI 





_ . „_i_n otS 

'•84oiS:« 


149+ 21 


1.31 


+0.68 
-0.66 


-302 


X 


/ 


MCG+8-11-11 


6.35+» : °3 

— U.U4 


0. 


- n +o. ur>6 
170 -0.042 




5.30 


+ 1.02 
-0.97 


-73 


X 


/ 


MR 2251-178 


6.46+^0 





0(17+0.199 
ou ' -0.100 


4a+ 14 


2.65 


+0.76 
-0.75 


-26 


X 


X 


Mrk 79 


6.17tf 2 2 7 





„___l_n 208 
Jza -0.159 




2.83 


+0.91 
-0.93 


-51 


/ 


/ 


Mrk 110 
















/ 


X 


Mrk 205 


— 






— 






— 


X 


X 


Mrk 279 


6.59l^f 9 




< 0.333 


20+?* 


0.15 


+ 0.19 
-0.14 


-6 


X 


X 


Mrk 355 


6.27±g-i? 





Knn+0.130 
500 -0.110 


134lg 


2.28 


+0.72 
-0.65 


-53 


/ 


/ 


Mrk 359 


6 40+ ' 06 
°-^ u _0 06 




< 0.177 


88+. 3 ? 


0.52 


+0.23 
-0.25 


-10 


/ 


X 


Mrk 509 


6 - 6 °to?4 





finn+i-363 
690 -0.151 


120+ 2 7 


6.32 


+ 1.50 
-1.42 


-57 


/ 


X 


Mrk 766 


6.66l° : ^ 





.+0.109 
144 -0.115 


7 q+27 
' — 25 


0.81 


+0.30 
-0.28 


-14 


/ 


X 


Mrk 841 


5.89i° :22 





. n9 +0.288 
4U -0.173 


80i 5 3 ° 7 


1.47 


+0.92 
-0.68 


-20 


/ 


X 


NGC 1365 


6.48t^ 2 




< 0.028 


< 13 


< 0.29 


-63 


/ 


/ 


NGC 2992 


6-51±g-it 





090+0.153 

ozo -0.100 


■70 + 28 

' ° — 25 


1.39 


+0.53 
-0.47 


-17 


X 


X 


NGC 3147 


6.45l» : ^ 




< 0.066 


no±|l 


0.21 


+0.10 
-0.08 


-15 


X 


/ 


NGC 3227 6 







7n7 +0.100 
lu '-0.087 


80±2| 


4.50 


+ 1.33 
-1.28 


-124 


X 


/ 










176+; 32 


5.38 


+0.97 
-0.94 


















1.17 


+ 1.13 
-1.13 
















28llg 


5.71 


+0.79 
-0.77 
















<55+ 23 


2.22 


+0.92 
-0.90 
















43+ 32 


1.35 


+ 1.00 
-0.98 








NGC 3516 b 


6-32±g;12 





, _ , in 110 


151+2? 


6.95 


+ 1.01 
-0.97 


-84 


/ 


X 










8ll 2 2 9 9 


1.89 


+0.69 
-0.68 








NGC 3783 6 


6.07+!!'?° 

— (J . 1 o 


0. 


+0.221 

761 -0.104 


97+oS 

— 32 


6.38 


+2.51 
-2.10 


-68 


X 


/ 










57±1| 


4.89 


+ 1.25 
-1.25 








NGC 4051 6 


6-23±g;15 





7. 9 + 0.151 
-0.124 


139+36 


1.80 


+0.50 
-0.47 


-70 


/ 


X 










74+^ 


2.12 


+0.68 
-0.66 
















112+. 3 , 


2.48 


+0.83 
-0.81 








NGC 4151 
















X 


X 


NGC 4593 


6-65±g"il 





o^o+0 289 
368 -0:i36 




1.26 


+0.58 
-0.51 


-12 


/ 


X 


NGC 5506 6 


6 50+° °I 

Ul -0.09 





317+0-074 
01 ' -0.060 


19+5 


7.83 


+ 1.91 
-1.51 


-32 


/ 


/ 












20+i 








NGC 5548 


6.36+°;°* 




< 0.098 


26+ 3 2 


0.60 


+0.72 
-0.52 


-5 


X 


X 


NGC 7213 


6 74+°' 16 
D - '^-0.16 




< 1.294 


19+30 


0.30 


+0.74 
-0.22 


-4 


/ 


/ 


NGC 7314 


6.36l u: ?| 




< 0.057 


co+37 
08_34 


0.66 


+0.43 
-0.39 


-7 


X 


/ 


NGC 7469 


fi 09+O.O6 





'■i5ot°:& 


fi9 +928 


1.55 


+0.70 
-1.10 


-24 


X 


X 


PDS 456 
















/ 


/ 


PC 1211+143 


6 49+ 07 
u -^ a -0.06 





1.9+0. 068 

-o.oei 




0.87 


+0.17 
-0.26 


-5 


/ 


X 


RBS 1124 
















X 


X 


SWIFT J2127.4+5654 


fi OO+0.24 


0. 


oqn+0.450 
Oau -0.260 


71+53 


3-43+11 


-12 


/ 


/ 


Ton S180 
















X 


X 



© 2012 RAS, MNRAS 000, ??-?? 



Suzaku observations of iron lines 



11 



MCG-02-1 4-009 




■ 


t , : 


4 5 6 7 8 9 
Observed Energy fkeV) 


MCG-6-30-15 

I 




1 5 6 7 8 9 
Observed Energy (keV) 



Mrk 79 










Observed Energy (keV) 



Observed Energy (keV) 



Observed Energy (keV) 



Mrk 335 










Observed Energy (keV) 



Observed Energy (keV) 



Observed Energy (keV) 



Figure 3 - continued Ratio plots of the 4-9 keV residuals without any modelling of the reflection component i.e. the FeK region is left 
totally unmodelled with only warm absorption at soft X-ray energies being taken into account in addition to the intrinsic powerlaw and 
COMPTT (where required). The data has been refit after the removal of components. The vertical dashed lines represent 6.4 keV in the 
observed frames. Red data points indicate those from further observations. 



3.2.3 Absorbed Seyfert spectra ~ partially covering 

A reasonable fraction of the AGN in this sample (16/46) re- 
quire the application of a partial covering geometry whereby 
some fraction of the underlying powerlaw is subject to an 
extra zone of absorbing gas (parameters tabulated in Table 
IB2|) . The majority of the partial covering scenarios here in- 
volve high column density warm absorbing zones which have 
their largest effect at harder X-ray energies. For example, 
some AGN still exhibit a hard X-ray excess at high energies 
> 10 keV after the use of a single reflionx unblurred dis- 
tant reflection component and require a partial coverer with 
iVH > 10 24 cm" 2 to fully model the HXD and BAT data, see 
Figure [5] In some cases extra modelling of the hard X-ray 
spectrum may be due to a second blurred reflector e.g. Mrk 



205, see Section 3.4. Due to the nature of a partial coverer 
with such a high column, there is not expected to be much 
interference with residuals in the FeK region except at the 
FeK edge (see discussion in Patrick et al. 2011b, Section 
3.3.1). In some objects, residuals remain above lOkeV even 
after applying a reflection component such as an unblurred 
reflionx can be significant, see Figure also see Turner et 
al. (2009) as example of the presence of a strong Compton- 
thick partial coverer. 

Five AGN (NGC 1365, NGC 3227, NGC 3516 and NGC 
4151) are best fit with a lower column density partial cov- 
erer with Nn ~ 10 23 cm -2 which adds curvature (i.e. a small, 
shallow bump) at > 2keV. Of course, in this region of the 
X-ray spectrum, there is likely to be some degree of inter- 
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Figure 3 - continued Ratio plots of the 4-9 keV residuals without any modelling of the reflection component i.e. the FeK region is left 
totally unmodelled with only warm absorption at soft X-ray energies being taken into account in addition to the intrinsic powcrlaw and 
COMPTT (where required). The data has been refit after the removal of components. The vertical dashed lines represent 6.4 keV in the 
observed frames. Red data points indicate those from further observations. 



play between the partial covering parameters (e.g. covering 
fraction, column density and ionization) and the strength 
of 'broad' residuals in the FeK region. Note that the par- 
tial coverer properties in NGC 1365 vary between each of the 
three Suzaku observations, in particular the second observa- 
tion (OBSID: 705031010) requires an additional neutral ab- 
sorber i.e. POWERLAW + XSTAR*ZPHABS*POWERLAW. Only 
IRAS 13224-3809 features a partial coverer with a lower col- 
umn density of Nn ~ 10 22 cm~ 2 which models small residu- 
als remaining in the 2-4 keV spectrum. The 2-10 keV Suzaku 
spectrum of IRAS 13224-3809 only has a small number of 
counts and subsequently no great deal of information can be 
gathered regarding any partial covering scenario. It should 



be noted that the majority of the partially covering zones in 
this sample are of low to neutral ionisation. 

3.2.4 Highly ionized absorption and outflows 

A total of 14/46 (30%) of the AGN in this sample fea- 
ture highly ionised absorption at a minimum of the 90% 
confidence level. A substantial proportion (10/46) of AGN 
in this sample are found to possess highly ionized zones 
which are somewhat blue-shifted with respect to absorp- 
tion lines from FeXXV and FeXXVI at rest-frame ener- 
gies of 6.7keV and 6.97keV (3C 111, 3C 445, MCG-06-30- 
15, MR 2251-178, Mrk 766, NGC 1365, NGC 3227, NGC 
3516, NGC 3783, NGC 4051, NGC 4151, NGC 5548, PDS 
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Figure 3 - continued Ratio plots of the 4-9 keV residuals without any modelling of the reflection component i.e. the FeK region is left 
totally unmodelled with only warm absorption at soft X-ray energies being taken into account in addition to the intrinsic powerlaw and 
COMPTT (where required). The data has been refit after the removal of components. The vertical dashed lines represent 6.4 keV in the 
observed frames. Red data points indicate those from further observations. 



456 and PG 1211+143). All of these zones (bar those in 
NGC 3783 and NGC 5548) are outflowing i.e. the absorp- 
tion lines are blue-shifted to higher energies with typical 
velocities of a few thousand kms -1 and column densities 
1 x 10 21 cm~ 2 < N H < 1 x 10 24 cm- 2 , see Table lB3l We 
note that the absorption zones in NGC 3227, NGC 3516, 
NGC 3783 and NGC 5548 have outflow velocities which are 
consistent with zero, i.e. 10/46 (22%) of the AGN in this 
sample require highly ionised outflows at greater than the 
90% confidence level. 



3.2.5 Soft excess and reflector properties 

The majority of the AGN in this sample show both soft 
and hard X-ray excesses above the basic intrinsic powerlaw. 
A substantial 31/46 AGN indicate an observed soft excess 
below ~ 2 keV (modelled as described in Section 3.1.1), how- 
ever this is lower than the fraction obtained in Porquet et al. 
(2004) in a selection of PG quasars who found that 19/21 ob- 
jects exhibited a significant soft excess. It should be noted, 
however, that nearly all of the objects without a soft excess 
feature complex warm or neutral absorption. This can lead 
to any weak soft excess which may be present being 'hidden' 
or simply absorbed, thereby leave little trace or indication of 
its presence. Only MCG-02-14-009 and NGC 3147 are with- 
out additional absorption (either neutral or ionized) and 
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Figure 3 - continued Ratio plots of the 4-9 keV residuals without any modelling of the reflection component i.e. the FeK region is left 
totally unmodellcd with only warm absorption at soft X-ray energies being taken into account in addition to the intrinsic powcrlaw and 
COMPTT (where required). The data has been refit after the removal of components. The vertical dashed lines represent 6.4keV in the 
observed frames. Red data points indicate those from further observations. 



still show no indications of a soft excess, whereas SWIFT 
J2127.4+5654 has sufficiently strong neutral absorption (no 
warm absorption) to mask the presence of a soft excess. 

A large number of AGN in this sample also exhibit 
strong excesses at hard X-ray energies, indicative of reflec- 
tion off distant material e.g. the putative torus or outer 
accretion disc. 39/46 of objects in the sample are mod- 
elled using the reflionx reflection model, with proper- 
ties such as ionization and iron abundance Zf b left free to 
vary. The vast majority of AGN modelled using reflionx 
here are well fit with a neutral or close to neutral reflector 
(£ < GOergcms -1 ) with only Ton S180 featuring a moder- 
ately ionized reflector (£ ~ 270 erg cm s -1 ). The remainder 
of the objects without a strong reflection component still 
feature neutral FeKa emission, however, IRAS 13224-3809 
features neither a reflection component nor 6.4 keV FeKa 
emission. 

In the majority of these AGN, the reflection component 
is well modelled with Solar iron abundance and in 23/39 
of objects with evidence for strong reflection it is fixed at 
Zp c = 1.0. This fraction includes the scenario in which some 
proportion of the observed narrow FeKa flux additionally 
arises from material which may be in the broad line region 
(BLR) i.e. the Fe Ka core also originates from both reflection 
off distant matter and the Compton-thin BLR. We should 



note that if we assume FeK emission from the BLR, this 
does not necessitate the addition of a further reflection com- 
ponent such as reflionx. This geometry is most evident 
when the strength of the narrow core is disproportionately 
stronger than the strength of other reflector features such 
as the Compton hump, often indicated by what may be an 
unfeasibly high iron abundance e.g. Zf c > 3.0 being forced 
into the fits. We account for this by freezing Z? e = 1.0 and 
including an additional narrow Gaussian ~ 6.4 keV (such 
as in MCG-05-23-16, Mrk 110, Mrk 509, NGC 2992, NGC 
3147, NGC 3516, NGC 4151, NGC 4593 and NGC 7213). We 
note, however, that if this additional FeKa component orig- 
inates from the BLR, it might be expected that the cr w idth 
be higher e.g. up to a ~ 0.10 keV. While the exact origin 
of the line may be unknown, we keep a — 0.01 keV for con- 
sistency with the other additional narrow lines used in the 
models, noting that some moderate broadening consistent 
with BLR widths cannot be ruled out. 

Only 4/39 of the remaining objects still require a 
slightly super-Solar iron abundance; namely 3C 120, Fairall 
9, Mrk 335 and Ton S180. Slightly sub-Solar iron abun- 
dance reflectors are found in 12/46 of objects and a good 
fit is found to both the Compton hump in the HXD and the 
flux of the narrow Fe Ka core with few residuals remaining. 
NGC 3227 and NGC 4051 (Lobban et al. 2011) are best fit 
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Table 4. Summary of rclline fits to objects in which some statistical improvement is made. Line energy fixed at 6.4 keV. * Denotes a 
frozen parameter. a Flux given in units (10 -5 phcm -2 s -1 ). b Inner radius of emission is quoted as an alternative to the spin parameter 
a, i.e. spin is frozen at maximal (a = 0.998) and the inner radius is allowed to vary away from -Risco- c Relline component is only 
employed in the 2006 of Mrk 766 observation. d Emissivity frozen at best-fit value from Patrick et al. (2011b). 
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with a sub-Solar iron abundance in conjunction with a nar- 
row 6.4 keV Gaussian. This is in a similar fashion to above 
whereby the narrow FeKa flux exceeds that expected from 
a simple Solar abundance reflector, however, both of these 
objects appear to still have weaker relative reflection compo- 
nents i.e. a sub-Solar reflector plus narrow Gaussian is pre- 
ferred. This may represent the scenario whereby the FeKa 



line either additionally or predominantly arises in Compton- 
thin material such as the BLR or NLR rather than the outer 
regions of the accretion disc or torus. The baseline model pa- 
rameters for each object and observation are summarised in 
Table [B2l and Table [B3l Of the objects which do not require 
a reflection component, only IRAS 13224-3809 does not fea- 
ture a narrow neutral Fe K line of any kind, while the narrow 
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Figure 4. Comparison of the differences between the soft X-ray 
modelling of the spectrum when using an XSTAR grid with tur- 
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Figure 5. Examples of the hard excesses which sometimes remain 
after modelling the full broad-band spectrum with a reflection 
component. In the baseline model we have accounted for these 
using high column density partial covering geometries. The black 
data points are XIS, the red are HXD data and the green are 
BAT data from Swift. 



iron line in both 3C 111 and MR 2251-178 is near neutral 
at ~ 6.5 keV. 



3.2.6 Additional distant ionized emission 

We examine the FeK regions of these AGN to determine 
whether emission lines due to ionized H-like and He-like 
species of iron are present, making up the final piece in 
our baseline models for each object. Lines at 6.63-6.70 keV 
and 6.97 keV are (likely due to FeXXV and FeXXVI re- 
spectively) emitted from highly ionized gas existing at large 
distances from the central SMBH. The presence of such lines 
has been well documented in many Seyfert 1 spectra and it 
should come as no surprise that they appear relatively com- 
mon in this sample (Bianchi et al. 2004; Nandra et al. 2007; 
Bianchi et al. 2009; Patrick et al. 2011a). Approximately 
half (24/46) of the objects in this sample show evidence for 
FeXXV emission at ~ 6.7 keV, while fewer (18/46) show 
evidence for FeXXVI at 6.97 keV (Table |B3|) . However, the 
lower fraction of H-like Fe detections may be due to its prox- 
imity to FeK/3 at 7.056 keV which is included self consis- 
tently in all objects in addition to the narrow 6.4 keV FeKa 
core. 



3.3 Broad residuals in the FeK region 

We first attempt to parameterise the strength of any broad 
residuals which may remain in the Fe K regions of these ob- 
jects by including a broad Gaussian with line energy, cr w i dth 
and normalisation left free to vary. If this broad Gaussian 
component proves to be statistically significant at the 90% 
level (three additional parameters i.e. an improvement to 
the fit of Ax 2 > 6.3), we then proceed to add more physical 
models of the line emission form the inner regions of the ac- 
cretion disc. Figure shows ratio plots of the 4-9 keV region 
prior to any model of reflection or disc emission components 
i.e. the FeK band is left unmodelled. 

Indeed, not all of the objects in the sample require any 
further modelling and an excellent fit is found to the data. 
For example, no reasonable statistical improvement is made 
when adding a broad Gaussian component to the follow- 
ing AGN in this sample: 1H 0419-577, Ark 564, MCG-02- 
14-009, MCG-02-58-22, Mrk 110, Mrk 205, Mrk 279, NGC 
4151, NGC 5548, NGC 7213, PDS 456, PG 1211+143, RBS 
1124 and Ton S180 (see Table |3J. The remaining AGN in the 
sample show at least some indications of a broad red-wing 
in the FeK region and therefore perhaps emission from the 
inner regions of the accretion disc (32/46 objects), however, 
it should be noted that this is only a simple parametrisation 
and the true fraction of AGN may differ from this when 
a more physical model is used. From Table [3] we can esti- 
mate the mean parameters of the simple broad Gaussian 
at LineE= 6.32 ± 0.04 keV, tT width = 0.470 ± 0.051 keV and 
equivalent width EW = 97 ± 19 eV. Of course, attempting 
to model an asymmetric line emission profile with a simple 
Gaussian can lead to a mis-modelling of the Fe K region, for 
example, broader Gaussians (such as those found here) are 
heavily influenced by the way in which the narrow FeKa 
core and Compton shoulder are modelled. Since the use of a 
broad Gaussian is essentially simply a means of adding cur- 
vature to the 5-7 keV region, in a sample of complex AGN 
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featuring warm absorption/ highly ionized outflows there 
is likely to be a reasonable amount of interplay between the 
width and strength of the broad line with the absorber prop- 
erties. 

We next replace the broad Gaussian with the relline 
relativistic line emission model (Dauser et al. 2010). This 
produces a more physical asymmetric line profile taking into 
account effects such as relativistic Doppler broadening pro- 
ducing both red and blue-wings. The shape of this line pro- 
file allows properties such as the inclination and emissivity 
index q of the accretion disc to be measured. We assume 
the emission line has a rest frame energy of 6.4 keV and 
that the inner radius of the accretion disc extends down to 
Risco with a single uniform emissivity index throughout the 
disc i.e. the inner and outer q values are tied. The emissiv- 
ity index scales as R~ q where q ~ 3 would be expected from 
the inner regions of the accretion disc and q > 5 would sug- 
gest that emission is very centrally concentrated, invoking a 
significant amount of light bending e.g. Miniutti & Fabian 
(2004). In some AGN, the broad line may not be sufficiently 
strong enough to constrain the emissivity index and in such 
cases we fix q = 3.0. The relline model also allows esti- 
mates to be placed upon the black hole spin parameter a 
(assuming R la — ifeco), which varies between a — —0.998 
for a maximally rotating retrograde SMBH and a — 0.998 
for a maximally spinning prograde BH. Objects with multi- 
ple observations are modelled allowing for the normalisation 
of the relline component to vary, but all other parameters 
tied e.g. in NGC 1365, NGC 3227, NGC3516 and NGC 5506. 
The relline component in 3C 120, Fairall 9, NGC 3783 and 
NGC 4051 is consistent with a constant flux held between 
observations (i.e. the line does not vary) and the normalisa- 
tion is tied as such, see Table [4] 

After applying the relline model a good fit is found in 
all the remaining objects, some of which are significantly im- 
proved with the introduction of a relativistic line emission 
model. Emission from the inner regions of the disc is not 
formally statistically required at the 90% confidence level in 
20/46 objects (see Table U for a summary of relline fits for 
objects which required a broad Gaussian in Table O note 
that not all objects listed in the table formally require the 
relline model). Of the remaining 26/46 objects which do 
require emission from the inner regions of the accretion disc 
(23/46 of the objects in this sample require the relline 
model at > 99.5% confidence), we find a moderate average 
strength of EW = 96 ± 10 eV and a low accretion disc emis- 
sivity index of q = 2.4 ±0.1 at an inclination of i — 33 ± 2° 
(see Figure [6|. Alternatively fixing a = 0.998 and allowing 
Ri n to vary yields an average R[ n — 21 ± 6 R s . 

With the advent of high quality broadband X-ray data, 
tentative steps can be made towards placing estimates upon 
black hole spin. In particular, objects with multiple or deep 
observations may possess spectra with a sufficient number of 
counts and temporal information with which to form models 
allowing for consistent spin estimates and the variation in 
spectral shapes between observations, e.g. see FigureO Here 
we make a total of 11 tentative constraints upon the spin 
parameter a (upper and lower bounds in 5 AGN), see Section 
4.3.3. 



3.4 Dual reflector fits 

If strong emission arises from reflection off the inner regions 
of the accretion disc producing a relativistically broadened 
Fe Ka line profile, the same degree of broadening and rela- 
tivistic affects can also be applied to the entire reflection con- 
tinuum. In this scenario we form dual reflector fits (for those 
objects which do formally require a relline component, see 
Table [4]) consisting of a distant unblurred reflionx and a 
second inner reflionx convolved with the relconv kernel 
(Dauser et al. 2010). However, we still include a COMpTT 
component to model the soft excess continuum. For example, 
the analysis of a small sample of 'bare' Seyferts in Patrick 
et al. (2011a) suggests that in order to produce the signifi- 
cant amount of blurring required to smooth the discrete soft 
emission lines into a continuum, both the spin parameter a 
and emissivity of the disc are forced to high and extreme val- 
ues (a = 0.998 and q > 4). Retaining the use of the COMpTT 
to model the soft excess ensures that the main feature driv- 
ing the fit of the blurred reflector are the broad residuals in 
the FeK region regardless of the interpretation of the soft 
excess. 

Due to the additional hard X-ray flux produced by a 
dual reflector model, we test this as an alternative to high 
column partial covering scenarios, i.e. partial covering ge- 
ometries are avoided during fitting unless no good fit can be 
found without a partially covering model. The iron abun- 
dance is tied between the inner and distant reflectors for 
consistency. The ionisation of the two reflectors are initially 
allowed to vary but tied if both are approximately equal, 
this is found to be the case in all objects. 

Simply replacing the relline component in each model 
with the relconv* reflionx convolution produces a good 
fit to the data in objects without partial covering, with es- 
timated accretion disc parameters consistent with those ob- 
tained in the previous model (Table [B4)l . This should come 
as no surprise since the Fe K region remains the main driver 
behind the fit. The dual reflector approach as an alterna- 
tive to partial covering provides a reasonable fit to Mrk 
766 (Ax 2 = 1109.6/1002), however, this is notably worse 
than when the use of a partial covering geometry is retained 
( X l = 1055.1/995). A poor fit is obtained to all three ob- 
servations of NGC 4051 without invoking partial covering 
(xl = 3321.9/2943 versus xl = 3177.1/2939) while keeping 
the iron abundance of both reflectors tied (Zf c — 1.0), how- 
ever, disc parameters are forced to extreme values (i.e. q > 7 
and a = 0.998). In addition to this, the 2-100 keV blurred 
reflector flux of the 2005 observation is approximately 2-3 
times that of the 2008 observations despite being in a far 
more absorbed state (see Figure [7J), we therefore retain the 
use of a partial coverer due to both physical and statistical 
arguments (see relevant section in the appendix for more 
details on NGC 4051). 

A poor and statistically inferior fit to the broadband 
data is found in some objects when a high column partial 
coverer is replaced with a dual reflector as an alternate mech- 
anism for accounting for additional hard X-ray flux above 
lOkeV. For those objects which feature a broad iron line, a 
partial covering geometry is still statistically required in Mrk 
766, NGC 1365, NGC 3227, NGC 3516 and NGC 4051 (Ta- 
ble lB4|l . Treating these AGN as above and replacing relline 
with relconv*reflionx while still including a partial cov- 
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erer again produces a good fit to the data with consistent 
accretion disc parameters to those in Table [4] (see Table 
IB4[I . Removing the partial coverer and accounting for the 
hard excess with the inner reflector yields an improved fit to 
MCG-06-30-15 and NGC 5506 (i.e. 2/16), however, the pa- 
rameters obtained are consistent with those estimated with 
a partial covering geometry. 



4 DISCUSSION 

This sample of 46 objects includes all the Seyfert 1 AGN 
matching our selection criteria outlined in Section 2.1 with 
observations which are publically available in the Suzaku 
archive. The aim of this paper is to form detailed broad- 
band models for each object, fully modelling and account- 
ing for any absorption where required which then allows us 
to parameterise the FeK regions and ultimately see in how 
many objects we can estimate or constrain SMBH spin. The 
main motivation for conducting this analysis is to produce 
an in-depth study of the FeK regions of these AGN and to 
assess the strength and prevalence of broad emission from 
the inner regions of the accretion disc. In order to accu- 
rately examine both emission and absorption properties in 
the FeK region, it is essential that high energy data is used 
to measure the strength of the Compton hump and hard X- 
ray excess as part of the formation of a broad-band model; 
the Suzaku X-ray observatory is unique its ability to provide 
such data simultaneously with soft X-ray data. The data in- 
cluded in this sample have been selected with sufficient S/N 
such that, if broadened line emission from the inner regions 
of the accretion disc exists (down to tens of eV in equivalent 
width), we should be able to detect it. 

It should be noted that the broad-band X-ray spectra of 
these AGN are complex with many of them featuring multi- 
ple absorption zones which have a significant role in adding 
spectral curvature and hence have a significant affect upon 
and 'broad' residuals in the Fe K region. Nandra et al. (2007) 
also note that the continua are complex in a study of the 
2.5-10.0 keV energy band with data obtained with XMM- 
Newton, however, here we take this study a step further 
by also considering a detailed modelling of the 0.6-10.0 keV 
soft X-ray data, allowing us to form more complete mod- 
els of the entire X-ray spectrum up to 100 keV. Low energy 
X-ray data below 2.5 keV is essential to perform a detailed 
analysis of the FeK region of any object since first forming 
an appropriate baseline model, with soft X-ray warm ab- 
sorption included, is an important prerequisite. This is due 
to the warm absorption zones which are required to model 
the X-ray data below 2.5 keV can add significant spectral 
curvature above 2.5 keV and therefore influence the FeK re- 
gion and residuals which may otherwise be interpreted as a 
broad red-wing from the very inner regions of the accretion 
disc, for example see MCG-06-30-15 (Patrick et al. 2011b). 

4.1 Ionized emission in the FeK region 

Narrow ionized emission from distant material is found 
to be relatively common (in the baseline model) amongst 
these Seyfert 1 AGN with 24/46 (52%) of objects featur- 
ing FeXXV emission at 6.63-6.7 keV and 18/46 (39%) fea- 
turing FeXXVI emission at ~ 6.97 keV. Meanwhile 10/46 



objects feature both FeXXV and FeXXVI emission (Table 
IB3|) . These fractions are much higher than those obtained 
by Nandra et al. (2007) who only find significant emission 
at 6.7 keV in NGC 5506 and significant 6.97 keV emission 
in both NGC 3783 and NGC 4593. Comparing to Fukazawa 
et al. (2011), we also find a higher fraction here, however, 
this may be due to our more stringent selection criteria (e.g. 
> 30000 counts and exposure > 50 ks) and hence an in- 
creased ability to detect a weaker ionised line in the data 
if it is indeed present. The findings here, however, concur 
with those found by Bianchi et al. (2009) and Patrick et 
al. (2011a) who similarly find that narrow ionized emission 
lines are a common feature of a number of Seyfert 1 spectra. 

After the inclusion of a component such as relline to 
account for any broad iron line residuals, however, these 
fractions can be somewhat reduced due to interplay with the 
blue-wing of a relativistically broadened line profile falling 
in the 6.6-7.0 keV region in some AGN. For example, a fea- 
ture at ~ 6.7 keV can often be modelled as either narrow 
ionized emission or as part of a broad line profile with fairly 
typical emissivity and inclination parameters. It is likely, 
therefore, that the true fraction of Seyfert 1 spectra with 
narrow ionized emission (FeXXV in particular) is lower af- 
ter the inclusion of a broad line component. In 4 objects 
the feature at ~ 6.7 keV is preferentially described by a rel- 
ativistic line profile for the above reasons, leaving the final 
number of AGN in this sample with distant FeXXV emis- 
sion as 18/46 objects. Calorimeter resolution spectra from 
Astro-H will help to resolve these issues. 

4.2 Absorption in the FeK region 

In addition to absorption at soft X-ray energies, absorption 
in the Fe K region due to ls-2p resonance lines from Fe XXV 
and Fe XXVI are found to be relatively common amongst the 
AGN in this sample. Evidence for highly ionized absorption 
is found 14/46 (30%) of objects: 3C 111, 3C 445, MCG-6- 
30-15, MR 2251-178, Mrk 766, NGC 1365, NGC 3227, NGC 
3516, NGC 3783, NGC 4051, NGC 4151, NGC 5548, PDS 
456 and PG 1211+143; consistent with the 36% detected 
by Tombesi et al. (2010b) in a larger sample of 101 XMM- 
Newton observations searching for ultra-fast outflows. The 
majority (10/14) of these lines are blue-shifted, suggesting 
that they originate from an outflowing zone of absorbing 
gas. When modelled with a high turbulent velocity xstar 
grid («turb = 1000 kms -1 ) we find a mean column density 
N H = (1.74 ± 0.85) x 10 23 cm~ 2 , mean ionization log(£) = 
3.97 ± 0.13 and a wide range in outflow velocity ranging 
between 400kms _1 < v out < 84600 km s -1 with no clear 
peak in the distribution of v out values (as shown in Figure[8]). 
More details regarding highly ionised and outflowing zones 
will appear in a forthcoming paper (Gofford et al. 2012, in 
prep), which performs a more exhaustive statistical search 
for blue shifted iron K absorption lines. 

We must note that these highly ionised zones can influ- 
ence the strength of any broad residuals in the FeK region, 
for example, Reeves et al. (2004) found that no broad residu- 
als remained in an XMM-Newton observation of NGC 3783 
after the application of the required xstar grid to model 
the highly ionized absorption. The measured strength of 
any 'broad' residual is therefore reduced given an appro- 
priate high £ absorption zone and hence the detection and 
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Figure 6. v Fv plots of the baseline+RELLlNE model in the 4-9 keV region. Note the strong absorption in the Fe K region in many objects 
and relatively weak red-wing. The solid black line shows the total model. 



modelling of such zones is essential if robust statistics re- 
garding the FeK region are to be formed. Based upon the 
work by Tombesi et al. (2010a & 2010b) the actual fraction 
of AGN with statistically significant outflows may indeed be 
higher than the 30% found here due to the presence of UFOs 
(Ultra Fast Outflows) with absorption lines blue-shifted to 
> 7keV since we have not systematically searched for these 
from 7-10 keV (see Gofford et al. 2012, in prep). The fraction 
of AGN with highly ionized zones presented here may also 
in fact be larger. It is possible that the presence of narrow 
emission lines from ionized species of iron may reduce our 
ability to detect absorption lines in CCD resolution spectra 
from the same species of iron with current instruments. 



4.3 Broad emission in the FeK region 

After the formation of a complete baseline model, we in- 
troduce the physically motivated relativistic line emission 
model relline to account for any remaining broad resid- 
uals in the FeK region. The criteria for the presence of a 
broad component is A\ 2 > 7.8 with the introduction of four 
additional free parameters (q, a, i and normalisation), this 
is equivalent to the 90% confidence level. Based upon this 
requirement, we find that 26/46 (i.e. 57%, Table [5} of the 
AGN (by object, not observation) in this sample formally 
statistically require relativistic emission from the inner re- 
gions of the accretions disc. Furthermore, 23/46 (50%, Table 
13} require broadened line emission at > 99.5% confidence 
(Ax 2 > 15). This fraction is somewhat lower (although not 
substantially so) than that obtained by Nandra et al. (2007) 
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Figure 6 - continued v Fv plots of the baseline+RELLlNE model in the 4-9 keV region. Note the strong absorption in the Fe K region in 
many objects and relatively weak red-wing. The solid black line shows the total model. 



who find that 18/26 (69%) of objects in their sample show 
evidence for broad line emission at > 99% confidence. The 
slight disparity between the two surveys may be due to the 
energy range over which the analysis has been conducted: 
0.6-100.0 keV here, versus 2.5-10.0 keV. This study also takes 
into account partial covering where required. However, the 
slight difference may be due to the smaller number statistics 
in the Nandra et al. (2007) sample. In a flux limited sam- 
ple, de la Calle Perez find 11/31 of objects (36%) feature 
relativistic iron lines. 



4-3.1 Average parameters of the disc 

Taking the parameters obtained with relline for the 26 
objects featuring broadened line emission (Table [4} we can 
estimate the typical parameters of the accretion disc (Table 



[BJ. The average inclination of the disc is i = 33° ± 2° for 
21 objects, consistent with Nandra et al. (2007) who find 
i — 38° ± 6°. The average emissivity index of the disc is 
measured at a low to moderate q = 2.4 ± 0.1 (for the 20 
objects in which it can be constrained) and is consistent with 
an analysis of 6 'bare' Seyferts by Patrick et al. (2011a) and 
an XMM-Newton survey by de la Calle Perez et al. (2010) 
who find q = 2.4 ± 0.4 and i = 28° ± 5°. This value is much 
lower compared to an often high emissivity (q > 5) assumed 
by many light bending models (Miniutti & Fabian 2004; 
Miniutti et al. 2009; Brenneman et al. 2011; Gallo et al. 2011; 
Nardini et al. 2011), suggesting that strong GR effects may 
not be present in the X-ray spectra of these Seyfert AGN. 
Figure [9] shows the distribution of the emissivity indices for 
the sample, note that there is a relatively small dispersion 
in q with the majority of values centred around q ~ 2.4. 
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Figure 6 - continued v Fv plots of the baseline+RELLlNE model in the 4-9 keV region. Note the strong absorption in the Fe K region in 
many objects and relatively weak red-wing. The solid black line shows the total model. 



The average strength of the relline component is mea- 
sured at EW = 96 ± 10 eV for the total of 26 objects, con- 
sistent within errors with Patrick et al. (2011a) and Nan- 
dra et al. (2007). The typical equivalent widths of the line 
component are also largely consistent with Fukazawa et al. 
(2011) who parameterise all of the 6.4 keV narrow and broad 
emission with a single Gaussian in a large Suzaku sample of 
Seyfert 1 and 2 AGN. As can be seen from Figure [5] there 
is a wide range in the measured equivalent widths of the 
relline component, revealing a possible bi- modal nature 
with a group centering around EW ~ 70 eV and a second 
stronger though less populated group at EW ~ 150 eV. Al- 
though there are insufficient statistics in this sample to draw 
any conclusions from this, indeed the distribution may be 
more continuous when taking error bars and different confi- 
dence levels into account. It is also interesting to note that 



in the cases where BH spin can be constrained, maximal 
spin is ruled out in all objects with the spin parameter a 
preferring to take low to intermediate values, similarly the 
typical inner radius of emission is at tens of R s with many 
clustered around _R ln ~ 20 — 30 7? g (see Figure 0. 

When parameterised with a Gaussian, we find an av- 
erage line centroid energy of LineE = 6.32 ± 0.04 keV (see 
Figure [TP)) . <T widt h = 0.47±0.05 and equivalent width EW = 
97 ± 19 eV (Table [6}, both of which are also consistent with 
a similar parametrisation with a broad Gaussian by Nandra 
et al. (2007) who find an average LineE = 6.27 ± 0.07 keV, 
EW = 91 ± 13 eV respectively. In an analysis of the 2- 
lOkeV spectrum of 149 radio-quiet Type 1 Seyfert AGN 
with XMM-Newton, de la Calle Perez (2010) estimate a 
mean broad line equivalent width of EW = 143 ± 27 eV, 
which is higher than that obtained here, however, the au- 
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Figure 6 - continued v Fv plots of the baseline+RELLlNE model in the 4-9 keV region. Note the strong absorption in the Fe K region in 
many objects and relatively weak red-wing. The solid black line shows the total model. 
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Figure 8. The distribution of ionization and outflow velocity of 
absorption zones at FeK energies. Note that the outflow veloc- 
ity covers a wide range velocities with no clear peak within the 
parameter space. 
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Figure 10. Histograms showing the distribution of line centroid 
energy and <r w jdth when broadened emission in the Fe K region is 
modelled with a Gaussian. 
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Figure 7. v Fv plots of the full broad-band models of some objects displaying significant variability. Note that the 2005 (black) 
observation of NGC 3516 is higher flux, yet appears to be more absorbed whereas in the majority of these examples it intuitively appears 
that simply the covering fraction has varied between observations. 



thors note that this may be an overestimate due to an in- 
ability to detect weak lines and a bias towards higher values. 

Note that when obtaining measurements for the emis- 
sivity index of the disc and subsequently the spin of the 
SMBH, the way in which the soft excess is modelled can 
have a significant effect. For example, if the soft excess 
is not modelled according to a soft Comptonization origin 
(e.g. compTT, or to a lesser extent diskbb) and instead an 
atomic origin is preferred, the disc parameters are drastically 
altered, see Section 3.1.1 . 

4-3.2 Shape of the line profile 

The emissivity index obtained here is based primarily upon 
any remaining residuals in the FeK region i.e. relline is 
modelling a relativistic line profile which takes into account 
both red and blue-wings of a line profile rather than sim- 
ply adding a broad 'hump' to the data below 6.4 keV (as, 



for example, a Gaussian might). It should be noted that for 
moderate to high inclinations (i > 20° ) the blue- wing of the 
line profile peaks at a higher flux than the red-wing peak, 
due to a boosting effect from material travelling at relativis- 
tic velocities in orbit 'close' to the central black hole. This 
effect is accentuated with lower emissivities since with in- 
creasing q (and hence more centrally concentrated emission) 
a greater proportion of the emitted photons are present in 
the red-wing of the line profile due to their proximity to the 
very inner regions and increased gravitational broadening. 
At particularly high emissivities, the blue-wing diminishes 
and the majority of the flux is contained within the red- 
wing of the line profile, forming a single peaked profile and 
an increasingly strong red-wing below 6.4 keV (Fabian et al. 
1989; Laor 1991; Dauser et al. 2010). Further still, at near 
maximal spin and high q the entire line profile can produce 
an asymmetric featureless hump which simply creates addi- 
tional curvature in the ~ 1 — 7keV region. 
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Figure 9. Figures of the distribution of emissivity index, relline equivalent width, SMBH spin and inner radius of emission. Note that 
the typical inner radius of the disc for the sample centres around tens of R g . 



Despite the narrow nature of the blue-wing at low a 
and q, its strength relative to the red-wing (and continuum) 
is typically insufficient to provide an alternative explana- 
tion for narrow emission lines commonly attributed to dis- 
tant FeXXV (at ~ 30°) or FeXXVI emission (at ~ 40°). 
However, there should be some caution to account for the 
presence of narrow ionized emission lines prior to applying a 
relativistic line model to ensure that the correct inclination 
of the disc is measured. 

The driving factor behind the estimated emissivity and 
spin values is therefore primarily the strength and width of 
the red-wing. Figure[lO]indicates that when modelled with a 
broad Gaussian, the typical line centroid energy of the broad 
line peaks at 6.3—6.4 keV which also suggests that any broad 
emission in the FeK region is not strongly gravitationally 
redshifted. The low average emissivity index (q = 2.4 ± 0.1) 
of the disc found in this sample of 46 objects is as a result of 
a full and thorough modelling of the broad-band spectrum 
i.e. after including xstar grids for the required absorption 
zones which dominate the spectrum in many objects be- 
low 2.5 keV, we note that the associated spectral curvature 
these zones add at higher energies (Miller et al. 2008; Zycki 
et al. 2010) leave little room for a significantly broadened 
relativistic line profile with high q and maximal spin. The 
extremely centrally concentrated emission that emissivities 
of q > 4 suggest are simply no longer required and relativis- 
tic emission from the disc (if present) originates from regions 
more distant to the BH (typically tens of R s , see Table [3] 
and Figure [SJ. 

The strength of the iron line profile has been noted to 
vary in a number of objects over long timescales and also 
over the course of a single observation (e.g. Matsumoto et 



al. 2003; Miyakawa et al. 2009; Risaliti et al. 2009). It is 
indeed possible that the analysis of the time-averaged spec- 
trum may instead produce an artificial broad line profile 
which is simply an average line profile which could therefore 
result in inaccurately derived disc properties. It is likely, 
however, that the changes in the observed line profile are in 
fact a consequence of changes in the warm absorber or par- 
tial covering properties. For example, over the six Suzaku 
observations of NGC 3227 the EW of the broad iron line is 
seen to vary greatly when the absorber properties are fixed 
between each observation (i.e. by hundreds of eV), however, 
allowing the ionisation of the partial coverer to vary sig- 
nificantly reduces the requirement for a strong broad line, 
instead resulting in typical EW ~ 50 eV. Allowing for such 
absorber changes on medium to long term timescales to ac- 
count for spectral variability can therefore drastically alter 
the measured Fe K line profile (Miyakawa et al. 2009; Risaliti 
et al. 2009). 



4.3.3 SMBH spin 

With the use of broad-band data from observatories such as 
Suzaku we can make tentative steps towards estimating the 
spin of the central supermassive black hole. The primary ef- 
fect of the spin parameter a upon the line profile is to set the 
degree to which the red-wing extends to the soft X-ray ener- 
gies, i.e. a maximal Kerr black hole would have the smallest 
inner radius of the accretion disc, thereby subjecting emit- 
ted photons to a stronger red-shift and sharply extending the 
red-wing to well into the soft X-ray regime. A Schwarzschild 
black hole (a = 0), however, truncates the red-wing of the 



© 2012 RAS, MNRAS 000, ??-?? 



Suzaku observations of iron lines 25 



Table 5. Table of main model components and the percentage 
of objects in the sample of 46 objects which feature those com- 
ponents based upon > 90% confidence levels. 



Component Percentage 



Soft excess 


67 


± 


6% 


Reflection 


85 


± 


6% 


Warm absorber 


59 


± 


5% 


Partial covering 


35 


± 


4% 


Fe K outflow 


30 


± 


4% 


Fe XXV emission 


52 


± 


5% 


FeXXVI emission 


39 


± 


4% 


Broad Fe Ka line 


57 


± 


5% 



line profile at ~ 4keV due to a minimum inner radius of 
6R g and less pronounced red-shifting of emitted photons 
(indeed a maximally retrograde BH pushes the inner radius 
further out to R[ n ~ 9R S ). It is clear to see, therefore, that 
if the spectral curvature associated with various warm ab- 
sorbing zones is not taken into account, a relativistic line 
emission model such as relline or kerrdisk may be forced 
to higher near maximal spin values to 'take up the slack' of 
an improperly modelled continuum. 

Many of the prime candidates for measuring SMBH spin 
with Suzaku have already been presented in previous papers, 
firstly the 'bare' Seyferts (Patrick et al. 2011a) and a col- 
lection of deep Suzaku observations (Patrick et al. 2011b), 
hence few new BH spin estimates are presented here. In 
Fairall 9, NGC 3516 and NGC 3783, however, we conduct for 
the first time an analysis of each of the multiple Suzaku ob- 
servations of each object, finding a baseline+RELLlNE model 
which can accurately describe the objects in different flux 
states with only subtle variations to the model over time 
e.g. changes in continuum flux. As stated in Section 3.3, all 
parameters of the relline component are tied between ob- 
servations except the normalisation i.e. the flux of the line 
is allowed to vary appropriately with the flux of the other 
primary components of the baseline model. Finding a model 
which can adapt over time to fit each observation while re- 
taining the same basic properties (e.g. emissivity, inclina- 
tion, spin) arguably increases the robustness of the results. 
Encouragingly, the fits to Fairall 9 and NGC 3783 are con- 
sistent with the analysis of the individual observations by 
Schmoll et al. (2009), Emmanoulopoulos et al. (2011) and 
Patrick et al. (2011a, 2011b), finding a = O.GOto.ll in Fairall 
9 ruling out maximal spin at > 95% confidence with an emis- 
sivity index of q = 2.9to'|. 

In addition to the spin estimates from Patrick et al. 
(2011a, 2011b), the results of this paper suggest further 
low-to- moderate SMBH spin values on MCG-05-23-16 {a < 
0.50), Mrk 79 (a < 0.80) and NGC 3516 (a < 0.30). Out of 
the total of 46 objects we find that spin constraints can be 
placed on 11/46 AGN, while 5 of these allow us to place up- 
per and lower bounds upon a and hence constrain the spin 
to some degree. Albeit with this small number of SMBH spin 
estimates, we can take tentative steps towards starting to de- 
velop a basic picture of the spin distribution of these Seyfert 
1 AGN (see Figure [9]). A maximal prograde SMBH is ruled 
out in all of these objects at a minimum of 90% confidence 
and > 99.5% confidence in 5/11 (MCG-6-30-15, NGC 3516, 
NGC 3783, NGC 7469 and SWIFT J2127.4+5654; Figure 



[9]). The dual reflector fits in Table [B4l also appear to suggest 
that 4 AGN in this sample may have retrograde SMBHs (i.e. 
a < 0; Mrk 79, NGC 3227, NGC 3516 and NGC 3783), how- 
ever, a more likely scenario is the possibility of a prograde 
rotating BH (or even Schwarzschild) with an accretion disc 
truncated short of the ISCO. 

4.3.4 NGC 3783 - high or low spin? 

In an analysis of the long 210 ks 2009 Suzaku observation 
of NGC 3783, Patrick et al. (2011b) (hereafter P2011b) 
concluded that the data appeared to rule out a maximal 
black hole in NGC 3783 (constraining spin to a < 0.31 in a 
dual reflector model), while the fit could be achieved with 
approximate Solar abundances for the reflector. However, 
Brenneman et al. (2011) (hereafter Br2011) came to the op- 
posite conclusion, appearing to require near maximal black 
hole spin from fits to the iron line data in NGC 3783 (with 
a > 0.98 at > 90% confidence) and obtaining a high iron 
abundance of Zf c — 3.7to'g and a steep inner emissivity 
law for the innermost disc. Subsequently Reynolds et al. 
(2012b) (hereafter Rey2012) also appeared to confirm the 
high spin, high abundance scenario in a re-analysis of the 
2009 dataset, suggesting that the iron abundance of the re- 
flector and the black hole spin may be degenerate upon each 
other and that a statistically preferred fit can be obtained 
with higher abundances and higher spin. 

In this section we attempt to discuss the differences 
between these works, by considering the long 2009 Suzaku 
observation of NGC 3783. Firstly the model of Br2011 was 
reconstructed, in order to understand the difference in the 
spectral modelling. The main difference between the Br2011 
and P2011b models is the construction of the warm ab- 
sorber; in Br2011 all 3 zones of the warm absorber only 
partially cover the AGN, where a fraction of ~ 17% of the 
direct continuum is unabsorbed by the warm absorber, or 
alternatively is scattered back into the line of sight. On the 
other hand in P2011b and in this paper, the warm absorber 
fully covers the X-ray continuum emission. Indeed we note 
that previous analyses have not needed to invoke partial 
covering in order to model the warm absorber in NGC 3783 
(Reeves et al. 2004; Yaqoob et al. 2005) , including high reso- 
lution X-ray spectroscopy from XMM-Newton RGS (Blustin 
et al. 2002) and during a 900 ks Chandra HETG observation 
(Kaspi et al. 2002). 

The other main difference in the model construction is 
that the pexmon neutral reflection model (Nandra et al. 
2007) is used for the distant (narrow) reflector in Br2011. In 
P2011b and in this paper, the reflionx ionised reflection 
model (Ross & Fabian 2005) is adopted, allowing the ionisa- 
tion state to reach a low value of £ = 1 erg cms -1 appropri- 
ate for low ionisation iron. For the disc (i.e. blurred) reflec- 
tion, both analyses use the reflionx table, convolved with 
a relativistic blurring function such as kerrconv (Brenne- 
man & Reynolds 2006) or relconv (Dauser et al. 2010). For 
simplicity in this section we use the kerrconv model, with 
spin allowed to vary between a — — 0.998 for a prograde 
black hole. The emissivity index is modelled as a broken 
powerlaw function (where qi is the inner emissivity and 52 
the outer emissivity and usually q\ > q2), breaking at a 
disk radius of r^ in units of R g . The inner disk radius is set 
equal to the ISCO, while the outer radius is set to 400i? g . 
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We adopt here the Solar abundances of Anders & Grevesse 
(1989) for the Galactic column, noting that there is little 
difference to the blurred reflector model parameters if the 
Wilms et al. (2000) ISM abundances are used instead. 

In order to test the Br2011 model, we use the identi- 
cal energy ranges for the XISFI spectrum and HXD/PIN 
adopted by Br2011, from 0.7-45 keV, ignoring the 1.5- 
2.5 keV band in the XIS. We also used the same warm ab- 
sorber model and tables as per the Br2011 and Rey2012 
papers. The Br2011 model has a steep index for the inner 
emissivity law, a high black hole spin consistent with max- 
imal and a high iron abundance of the inner reflector of 
■Z'Fe, inner ~ 3, while the iron abundance of the distant re- 
flector is initially fixed equal to Solar, i.e. Zf c , outer = 1. 
This model as detailed in Br2011 initially gives a poor fit 
(xt = 1501/1234), however subsequently refitting the model 
parameters then gives an excellent fit to the data, with 
xt = 1329/1234. Nonetheless the model parameters are in 
good agreement with those obtained in Br2011, where we 
obtain a formal 90% lower limit of a > 0.82 to the black 
hole spin (for 1 interesting parameter), a disc inclination 
of i — 22 ± 4, while the emissivity indices have values of 
qi = 4.4 ± 1.2 and q-2 = 2.8 ± 0.2, with a break radius of 
~6i? g . 

However upon taking the same model and fixing the 
spin to a — 0, then the fit statistic obtained is only slightly 
worse, where xt ~ 1340/1237. In this case only a single 
emissivity is required, where qi = q2 — 3.0 ± 0.5, while the 
disc abundance is consistent with Solar. Given that the high 
black hole spin model only shows a marginal improvement 
in fit statistic (Ax =11 for 3 fewer degrees of freedom), 
then the claim of high black hole spin cannot be confirmed 
at a high confidence level, compared to the case where a ~ 0. 
Furthermore if the Br2011 model is altered such that all 3 
soft X-ray warm absorber zones fully cover the AGN, then 
this leads to a lower (and less constrained) value for the black 
hole spin, of a = 0.35±g;5jj and is then formally consistent 
with the upper-limit of a < 0.31 from the dual reflector 
model in P2011b. 

The lower spin value may be due to the fact that the 
fully covering absorber adds greater spectral curvature to 
the model due to bound-free absorption, compared to the 
partial coverer. In the latter case, the absorption is diluted 
by an unabsorbed powerlaw, reducing the amount of spec- 
tral curvature, while the breadth of the iron line profile can 
increase to compensate. Finally we note that replacing the 
pexmon model with a reflionx table for the distant re- 
flector made little difference to any of the blurred reflection 
parameters and this appears unaffected by the parametrisa- 
tion of the distant reflection component. 

We also attempted to recreate the maximal spin sce- 
nario subsequently presented by Rey2012, but applied to 
the model as constructed in P2011b and in this paper. Thus 
in this case we assumed that all 3 zones of the warm ab- 
sorber in NGC 3783 fully cover the AGN as per P2011b 
and use the same warm absorber grids as in this paper, 
although the choice of a particular absorption model grid 
appears to not effect the reflection models. For ease of com- 
parison we also adopt the energy ranges used by Br2011 as 
above. We focus in particular on the iron abundances, which 
Rey2012 suggest is critical for determining the blurred re- 
flection parameters. The values for the disc reflector param- 



eters from Table 1 in Rey2012 are used, namely the iron 
abundances, spin, emissivities and break radii, all with the 
inner and outer radii set to the ISCO and 400i? g respec- 
tively. In summary we test three scenarios for the 2009 data: 

(i) ^Fc, inner = ^Fe, outer = 3.3 with near maximal spin; (ii) 

^Fe, inner = ^Fe, outer = 1 With a = 0; (ifi) ^Fe, inner = 4.2, 

^Fe.outer = 1 with near maximal spin. These correspond 
closely to the models A, B and C respectively, as presented 
in Table 1 of Rey2012. The kerrconv model is used to blur 
the disc reflection spectra, while an unblurred low ionisation 
reflionx grid is maintained to model the distant reflector. 
Aside from the blurred reflection parameters as per Rey2012, 
the fit parameters are allowed to vary. 

All three scenarios give a good fit to the 
data, with reduced chi-squared values of xt = 
1300/1227,1280/1227,1284/1227 for models (i), (ii) 
and (iii) respectively. Thus all the spin and abundance 
scenarios are statistically acceptable, while the a = and 
Solar abundance case formally gives the better fit. Allowing 
the disc parameters to vary for the Solar abundance case 
then gives a constraint on the spin of a < 0.45 at 90% 
confidence, while the disk emissivity is qi = q2 — 3.0 ± 0.5 
and no break radius is required. Indeed the same Solar 
abundance model when fitted over the energy ranges used 
in this paper also gives identical results and an equally 
good fit, with xt = 1414/1374, see Table lB4l 

We also compare the fit focusing in on the iron K band 
in Figure rm where the top panel shows the best fit dual re- 
flector model presented here (with a — 0) and the lower two 
panels show the data/model ratios for the scenarios (i) and 

(ii) above, with high spin and low spin respectively. Clearly 
there is little difference in the residuals between either the 
high spin scenario with super Solar abundances and the low 
spin scenario with Solar abundance. The only marginal dif- 
ference is the narrow core of the Fe Ka line at 6.4 keV is 
slightly over-predicted in the high spin, high Fe abundance 
scenario. This is in contrast to the plots shown in Figure 5 
of Rey2012, where the a = scenario strongly over-predicts 
the red-wing of the line between 5-6 keV. This could be due 
to the fact that the blurred reflector has a higher normalisa- 
tion in the high spin scenario, which when applied directly 
to the a = case could appear to over-predict the red-wing, 
unless the normalisation of the reflector is refitted accord- 
ingly. 

Finally we note that if the relline model is used in- 
stead of a blurred reflection component, then a lower spin 
value is usually preferred. For example, the upper-limit on 
spin to NGC 3783 with this model is a < 0.24 (also see Table 
4). Indeed P2011b also noticed some tendency (although not 
statistically significant) for the line to favour a retrograde 
black hole spin. The lower spin value may be due to the fact 
that the model only fits the broad iron line and not the re- 
flected continuum, where the latter can often be blurred to 
such an extent that it is hard to distinguish from the direct 
continuum. 

Thus the overall conclusion would appear to be that it is 
very difficult to definitively determine the black hole spin in 
NGC 3783 with the present data, given the complexities of 
the models involved and the deep warm absorber present in 
this AGN. For instance as discussed above, the construction 
of overall the model and the warm absorber in particular can 
have an effect. In addition the spin value is indeed degener- 
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ate upon the iron abundance, as discussed in Rey2012 and 
thus the high spin cases always require high centrally con- 
centrated iron abundances. Nonetheless a statistically good 
fit is obtained with a simple Solar abundance reflector and 
no black hole spin. The evidence for a more complex picture 
with high abundances, complex emissivity profiles and high 
spin, would appear not to formally be required at high con- 
fidence, although neither can the high spin case be ruled out 
at present. 

4-3.5 Truncation of the disc 

As can be seen from Table 2] and Figure O the majority 
of AGN in this sample do not feature strong emission from 
the very inner regions of the accretion disc, both the emis- 
sivity index and spin of the central BH take typically low 
values. This suggests that the dominant regions producing 
the typical broad iron line profile are regions close to but 
more distant from the central BH than has been estimated 
in many analyses (Reynolds et al. 2005; Brenneman et al. 
2006, 2011; Miniutti et al. 2007; Reynolds et al. 2012b), orig- 
inating at typically tens of R s and hence with a significant 
amount of flux present in the blue-wing of the line profile 
as opposed to the red-wing. The relline fits in Tableware 
under the assumption that the inner radius of the accretion 
disc extends down the inner-most stable circular orbit (i.e. 
^isco), if the majority of the line profile originates from 
the outer regions of the disc, it is plausible that in some 
cases the accretion disc may be truncated at a few to tens of 
R g . For example, Lobban et al. (2010) suggest a truncated 
accretion disc in NGC 7213, stating a lack of a significant 
broad component in the Fe K region and weak reflection. If 
the inner edge of the accretion disc does not extend down 
to the innermost stable circular orbit, it would make the 
determination of SMBH spin through spectral fitting very 
difficult. Essentially at i£j n > 9 R s the spin could take any 
value —0.998 < a < 0.998 with the disc truncated at any 
Ri n with little difference upon the relativistic line profile. 

4-3.6 Alternative mechanisms 

Recent work by Tatum et al. (2012) suggests that an al- 
ternative mechanism to reflection from the inner regions of 
the accretion disc is emission from a Compton thick disc 
wind based upon a disc wind model by Sim et al. (2010). 
The authors take the cleanest sources possible i.e. the 'bare' 
Seyfert sample from Patrick et al. (2011a) and assuming a 
relatively face-on viewing angle whereby only the scattered 
light from the wind is observed rather than any absorption 
features with a more edge-on viewing angle. The line profiles 
are therefore subject to velocity broadening due to Doppler 
velocity shear across the wind and Compton down-scattering 
of the Fe Ka flux, these can combine to imitate a broad line 
profile from reflection off the inner regions of the accretion 
disc. Good fits are made to the broadband spectra for this 
model without the requirement for any further broadened 
emission from reflection off the inner regions of the disc. In 
this scenario, the full FeK line profile is reproduced with 
Solar iron abundances at typically tens to hundreds of R g , 
however, in some objects an additional cold neutral reflector 
and neutral FeKa line component is still required in order 
to model the FeK edge and hard X-ray excess. 
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Figure 11. Upper panel vFv 4 — 9keV plot of the Patrick et 
al. (2011b) best- fit model (with a = 0) to the 2009 NGC 3783 
Suzaku data, note that the Fe K region is well modelled in the fit 
as is the broad red-wing below 6.4 keV. Middle panel Residuals 
to the best fit dual reflector model, as shown in the upper panel 
with a = 0. Lower panel Residuals from the model with inner 
disc parameters as obtained by Reynolds et al. (2012b), i.e. with 
maximal spin, high emissivity and Zp c > 3. Both of the residual 
plots in the lower two panels show no noticeable difference in the 
iron K region, thus it appears difficult to discriminate between 
the two scenarios. 
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While the disc wind model alone can start replicate 
many of the features in the broad-band spectrum, it is likely 
that the true scenario is a combination of absorption, reflec- 
tion and scattered components off both the disc and a wind, 
a scenario which may not be resolved with current X-ray 
observatories without calorimeter resolution spectra in the 
FeK band. It is, however, interesting to note that the tenta- 
tive bi-modality of the broad line equivalent widths in this 
sample (Figure 0) may represent such a scenario whereby a 
fraction of the weaker observed broad lines are as a result 
of disc wind geometries. Meanwhile objects which display 
stronger broad lines (EW > 150 eV) may be the only ob- 
jects in which we can see the inner tens of R s of the disc. It 
must be stressed, however, that the apparent (albeit weak) 
bi-modal nature of the equivalent widths in this sample may 
simply due to insufficient statistics. 

4.4 Distant reflection 

The use of the HXD instrument on Suzaku and BAT on- 
board Swift allows us to examine the hard X-ray spectra 
of the objects in this sample and to appropriately parame- 
terise the distant reflection component. This may arise, for 
example, from the reflection of continuum X-rays off the 
cold (T < 10 6 K) outer regions of the accretion disc or via a 
parsec scale torus. 

The majority of AGN in this sample show evidence for 
a reflection component, only 7/46 (see Table [5} do not re- 
quire the addition of the unblurred reflionx model (3C 
111, IRAS 13224-3809, MR 2251-178, Mrk 79, NGC 7314 
and PG 1211+143). Neither are these objects best described 
with a high column density partially covering absorbing zone 
with A^h ~ 10 24 cm -2 which can replicate a hard X-ray ex- 
cess to some extent i.e. there is no apparent hard excess in 
these objects. Table IB5I also indicates that the large ma- 
jority of AGN in this sample do feature excess emission in 
the 15-50 keV region over the extrapolated 2-10 keV intrin- 
sic powerlaw. A number of the objects in the sample fea- 
ture particularly strong reflection with 15-50 keV reflector 
flux Freflector > 3 x 10" 11 erg cm -2 s" 1 , namely IC 4329A, 
MCG+8-11-11, NGC 3227, NGC 3783, NGC 4151 and NGC 
5506. In addition to this, a large number of the AGN in 
this sample feature particularly hard X-ray spectra i.e. those 
with a high hardness ratio Jis-sokev/^b-iokeV ^ 3: 3C 445, 
Mrk 279, NGC 1365, NGC 3227 (Obs 2, 4 & 6), NGC 3516 
(Obs 1) and NGC 4151 (see Table |B5|) . 

Reflected flux 
15 50 Full model flux — Re fleeted flux 

Taken from values in Table IB5I we also calculate a re- 
flection fraction over the 15-50 keV range using Equation [T] 
note that this J?i5_so is not the same as the reflection frac- 
tion R included in models such as pexrav (Magdziarz & 
Zdziarski 1995). For ease of comparison, we calculate the 
equivalent i?i5_so value for an R — 1 pexrav versus a sim- 
ple unabsorbed r = 2 powerlaw, this yields Risso ~ 0.8. 
The reflection fraction is particularly high in 8 observations 
in the sample and four objects (MCG-02-14-009, NGC 3227, 
NGC 3783 and NGC 4151). In these objects the reflector flux 
in the 15-50 keV band contributes approximately half of the 
total 15-50 keV flux, however, this appears to bear no cor- 
relation with properties such as covering fraction although 



each of these four objects do feature relatively high hardness 
ratios (ratio of 15-50 keV to 2-10 keV flux). The majority of 
AGN, however, feature reflection fractions fitting a smooth 
distribution between J?i5_5o = and _Ri5_5o = 0.7 (see 
Figure I12[) . The measured strength of the reflection com- 
ponent can be somewhat degenerate with the method in 
which the partial coverer is modelled. However, the reflec- 
tion continuum below lOkeV can help to ensure that the 
correct amount of hard X-ray flux is produced by reflionx 
by measuring the strength of the narrow 6.4 keV core, for 
an appropriate iron abundance, the self-consistent reflec- 
tion strength is then modelled above lOkeV. Any remain- 
ing residuals > lOkeV must then be due to a high column 
density partial coverer. That is, the strength of the narrow 

6.4 keV core can give a reasonable estimate of the contribu- 
tion of the reflector above 10 keV. 

A significant number of AGN in this sample are fit very 
well assuming a simple Solar iron abundance i.e. Zp e = 1.0 
when allowed to vary as a free parameter within reflionx. 
Although the narrow FeKa core is sufficiently strong in at 
least 9 objects that further neutral iron emission is required 
from material distant to the central SMBH, e.g. the broad 
line region in addition to reflection off the torus. Four objects 
strongly require sub-Solar iron abundance (IC 4329A, NGC 
3227, NGC 5506 and SWIFT J2127.4+5654) while 4/39 sta- 
tistically prefer a slight super-Solar iron abundance (typi- 
cally Zfb < 3). This suggests that the hypothesis of Solar 
abundances throughout the disc and central engines of these 
AGN is a reasonable assumption in contrast to super-Solar 
Zpe based upon inner disc reflection (e.g. Brenneman et al. 
2011; Fabian et al. 2012). The analysis of the 46 objects in 
this sample indicates that the reflection in an overwhelming 
majority of the Seyfert 1 type AGN in the Suzaku archive 
is well described by reflection of a cold, distant, neutral ma- 
terial with Solar abundances throughout, also see Nandra & 
Pounds (1994) and Rivers, Markowitz & Rothschild (2011a) 
for similar findings regarding the recurrent presence of the 
Compton hump. We note that in higher 2-10 keV luminosity 
AGN we find that the equivalent width of the narrow Fe Ka 
emission line is low in comparison to other lower luminosity 
AGN, although the sample in this paper lacks the sufficient 
statistics to support the evidence for the X-ray Baldwin ef- 
fect as suggested by Jiang et al. (2006) and Fukazawa et al. 
(2011). 

4.5 Warm absorbers and partial covering 

The way in which absorbing zones are modelled in X-ray 
spectra can have a significant affect upon the FeK re- 
gion and the estimation of properties relating to the cen- 
tral SMBH and the surrounding accretion disc. Within this 
Suzaku sample, a large proportion of AGN feature complex 
absorbers, requiring one or more zones of an xstar grid. 
Only 12/46 (26%) of objects in this sample can be con- 
sidered 'bare' for the purposes of this paper - featuring no 
statistically significant additional absorption zones in ad- 
dition to the neutral Galactic column, these are 3C 390.3, 
Ark 120, Fairall 9, MCG-02-14-009, Mrk 110, Mrk 335, Mrk 
359, NGC 3147, NGC 7213, NGC 7469 and RBS 1124. A 
few objects in the sample (6/46, 9%) feature an additional 
neutral zone of absorbing gas at the redshift of the object, 
while the remaining 27/46 (59%) all require the application 
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Figure 12. The distribution of reflection fractions for each obser- 
vation in the sample with values obtained from Tablc lBSl whereby 
the fraction is the ratio of the 15-50 keV reflected flux to the 15- 
50keV continuum flux (i.e. full model minus reflected flux). Note 
that the majority of AGN feature reflection fractions between 
and 0.7. The final bin from 1.1 upwards contains only three ob- 
servations of AGN with i?i5_5o = 1.6, 3.5 and 7.2. Ri5_5o = 0.8 
is equivalent to a R = 1 PEXRAV versus a simple unabsorbed 
POWERLAW. 



Table 6. Table of average parameters of the sample, values 
quoted for a broad Gaussian and relline are for those objects in 
which those components were statistically significant at > 90% 
confidence level. 



Parameter 



Average value 



POWERLAW 
T 1.82 ±0.03 

Broad Gaussian 
Energy (keV) 6.32 ± 0.04 
EW (eV) 97 ± 19 

^width (keV) 0.47 ± 0.05 

RELLINE 

EW (eV) 96 ± 10 

q 2.4 ±0.1 
R in (R g ) 21 ±6 

Inclination (°) 33 ± 2 



of complex warm (i.e. ionized) zones of gas to model the 
broad-band spectrum. 

Just over a third of the sample (16 AGN) also statisti- 
cally require a partially covering geometry (i.e. 35% of the 
total sample, Table |3J. In this scenario one of the absorbing 
zones of gas (typically high column Nu > 5x 10 23 cm -2 ) par- 
tially obscures the nucleus in addition to one or more lower 
column density fully-covering absorbing zone. The majority 
of the partial covering models in this sample have a suf- 
ficiently high column density such that their predominant 
effect is to supplement and increase the hard X-ray flux. 
The distribution of covering fractions for the partial cover- 
ing fits here cover the full range from low to high covering 
fractions (Figure [13)) . 

We note that only single layer partial covering geome- 
tries are used (where formally required) in this analysis since 
multiple layers can be somewhat ambiguous and care must 
be taken to ensure that layers are not added on an ad-hoc 
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Figure 13. The distribution of partial covering covering frac- 
tions. Note that the full range is covered with a relatively even 
distribution. 



basis until all features are removed from the spectrum. None 
of the objects in this sample require more than a single layer 
of partial covering in order to obtain an acceptable fit to the 
broad-band data. 

Partial covering scenarios in some previous analyses 
have been used to model spectral curvature in the FeK re- 
gion which may otherwise be attributed to relativistic line 
emission e.g. MCG-6-30-15 by Miller et al. (2009). However, 
careful examination of the effects of high column partial cov- 
erers such as those used here suggests that the measured 
strength of broad line emission in the FeK region is sim- 
ply reduced (i.e. lower EW) while the estimated accretion 
disc parameters and BH spin estimates remain relatively un- 
changed (for example MCG-06-30-15 in Patrick et al. 2011b, 
Section 3.3.1). During the analysis of this sample we have 
also noted that the application of high column density par- 
tial covering scenarios in fact has some bearing upon the 
iron abundance of the reflector derived from reflection mod- 
els such as reflionx (or indeed pexmon). Since the par- 
tial coverer produces additional flux at hard X-ray energies, 
it is possible that the reflection continuum is subsequently 
under-predicted and the iron abundance is increased to com- 
pensate for the lack of narrow Fe Ka core flux. However, the 
majority of the objects here in which we require additional 
narrow Fe Ka emission through the use of an additional nar- 
row Gaussian do not feature partial covering and therefore 
it seems likely that the additional narrow line flux in these 
particular objects originates from distant material, as op- 
posed to being a consequence of improper modelling of the 
spectrum above 10 keV. 



4.6 Inner reflection versus partial covering 

The baseline model constructed in the main models of this 
paper makes use of a high column partial coverer to account 
for any additional hard X-ray flux which may be present 
in the X-ray spectra of these AGN. An alternate scenario 
is to use a blurred reflection component from the inner re- 
gions of the disc in addition to an unblurred reflector from 
the outer regions of the disc can provide a means of supple- 
menting the total hard X-ray flux at > 10 keV. As noted in 
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Section 3.4, this scenario can provide a good description of 
the broad-band spectrum for MCG-06-30-15 and NGC 5506 
whilst retaining relatively typical accretion disc parameters 
i.e. low to moderate emissivity index and inclination, similar 
to those obtained from an analysis of the broad Fe Ka line 
profile (see Section 3.3 and Table 

The fit to the remaining objects (Mrk 766, NGC 1365, 
NGC 3516, NGC 3783 and NGC 4051) is worse without a 
partial covering geometry. This may be due to the inability 
of the inner reflector geometry to successfully describe the 
long term variability (e.g. variability in spectral curvature) 
of these AGN, all of which are noted to have varied between 
each of the observations analysed here. In each of these AGN 
the disc parameters are either poorly constrained or are 
forced to extreme values, for example, in order to describe 
the broad-band spectrum and long term variability of NGC 
4051 the emissivity of the disc is increased to q = 6.llg'j 
and spin to near maximal a > 0.99, while producing an 
overall worse fit. This suggests that emission is required to 
be very centrally concentrated with the inner reflector prov- 
ing the dominant component in the broad-band model. For 
example, the inner reflector for 2005 observation of NGC 
4051 has 15-50 keV flux ~ 8.8 x higher than the outer re- 
flector (F2005 inner = 0-.49±°-%l) x 10" 11 erg cm* 2 s~ 1 versus 
Faoosouter = (0.17±^m) x lO-^ergcm-^s- 1 ). Comparing 
this to the alternative (i.e. with partial covering explaining 
the long term spectral variability, Table [B4| . the inner reflec- 
tor 15-50 keV flux is significantly reduced when partial cover- 
ing is invoked: F2005 inner = (0.29±g;g|) x 10 _11 erg cm" 2 s" 1 
compared to F 20 os outer = (0.37^5) x lO-^ergcm" 2 s" 1 . 

We find here that in the majority of cases, the dual re- 
flector scenario does not provide a reasonable alternative to 
a partial covering geometry particularly in objects which fea- 
ture significant long term spectral variability (see Table [B4| . 
NGC 1365 is very poorly fit with xt ~ 2.2 above 2keV and 
xl ~ 5.1 over the 0.6-100 keV range; NGC 3227 and NGC 
3516 also featuring poor fits with xt ~ 1-2 and xt ~ 1-8 
respectively when fitted without partial covering. It should 
be noted that many of the iron line parameters in Tables [B4l 
are poorly constrained, this is due to the interplay between 
the inner reflector and outer reflector (and/ or partial cov- 
ering) at hard X-ray energies and the COMPTT component 
at soft X-ray energies. 



5 CONCLUSIONS 

Based upon an analysis of all publically available Suzaiu 
observations of Seyfert 1 AGN with observations longer than 
50 ks and greater than 30000 XIS counts, we conclude the 
following: 

(i) The majority (59%) of AGN in this sample feature 
complex warm absorption which has a significant affect upon 
the FeK region and any accretion disc parameters derived 
from it. The use of the full 0.6-100.0 keV broad-band data 
is therefore essential prior to any attempt to use relativistic 
line profile models as a diagnostic for the inner regions of 
the accretion disc. The mean photon index of the 46 objects 
on the sample is T = 1.82 ± 0.03. 

(ii) Absorption in the FeK region due to highly ion- 
ized gas producing absorption features from FeXXV and 
FeXXVI are relatively common in AGN (30%), most of 



which (71%) are outflowing at high velocities. While a large 
fraction of the detected highly ionised winds are outflowing, 
there may be a larger number of low velocity winds which are 
not detected due to the possible presence of ionised Fe XXV 
and FeXXVI emission lines. The additional curvature added 
to the region through modelling with an appropriate xstar 
grid, while subtle, has a notable effect upon the strength of 
any broad residuals which may remain below 6.4 keV and be 
interpreted as strong relativistic emission. 

(iii) A partial covering geometry is required in 35% of all 
objects in the sample. These high column density zones pri- 
marily affect the hard X-ray spectrum above 7 keV although 
reducing the strength of broad residuals in the Fe K region 
rather than removing them entirely. 

(iv) Narrow ionized emission in the FeK region from 
FeXXV and FeXXVI are relatively common in these AGN, 
featuring in 24/46 and 18/46 of objects respectively. Of these 
AGN, 10/46 feature both FeXXV and FeXXVI emission. 
These lines are found to be much more common compared 
to an XMM-Newton survey by Nandra et al. (2007) de- 
spite possible interplay with the blue-wing of more sophis- 
ticated relativistic line emission models which could reduce 
the number of narrow line detections. 

(v) Examining the Fe K region after a complete modelling 
of the broad-band 0.6-100.0 keV spectrum and all required 
absorption zones yields a range of weak to moderate strength 
broad residuals below 6.4 keV. We find that 26/46 (56%) of 
the objects in this sample require some degree of relativistic 
line emission in the Fe K region at 90% confidence and 23/46 
(50%) at > 99.5% confidence. 

(vi) These broad residuals are well fit with the relline 
model and yield an average broad line strength of EW = 
96 ± 10 eV for the total of 26 objects. The line energy and 
Cwidth of the broad residuals when modelled with a Gaussian 
are consistent with Nandra et al. (2007). 

(vii) We estimate an average emissivity index of the ac- 
cretion disc of q = 2.4 ± 0.1, suggesting that emission from 
the accretion disc responsible for relativistic lines is not ex- 
tremely centrally concentrated when purely the line profile 
in the FeK region is used as a diagnostic. The majority of 
the line flux therefore occurs from the blue-wing of the line 
profile with emission being insufficiently close to Risco as 
to redshift a significant proportion of the X-ray flux into a 
strong red-wing. We also measure an average disc inclination 
of i = 33°±2° and inner radius of emission R ln — (21±6) i? g . 

(viii) With the assumption that the inner radius of the 
accretion disc (Rm) extends down to the innermost stable 
circular orbit (-Risco), loose constraints upon the SMBH 
spin parameter a can be made. The relativistic line emis- 
sion profiles are sufficiently distinguished in 11/46 objects to 
place upper or lower bounds on the spin. After a broad-band 
analysis we make the following estimates: Ark 120, a < 0.94; 
Fairall 9, a = 0.60t£;^; MCG-02-14-009, a < 0.88; MCG- 
05-23-16, a < 0.50; MCG-06-30-15, a = 0.49±°; 2 °; Mrk 79, 
a < 0.80; Mrk 335, a = 0.70±oJ 2 ; NGC 3516, a < 0.30; 
NGC 3783, a < 0.24; NGC 7469, a = 0.69^9 and SWIFT 
J2127.4+5654 with a = 0.70±g;^- Under the assumption 
that Ri n — -Risco , a maximally rotating SMBH is ruled out 
in each of these 11 objects. 
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APPENDIX A: SELECTED INDIVIDUAL 
OBJECTS 

Al Fair all 9 

The FeK region of the 'bare' Seyfert Fairall 9 has been well 
studied (Schmoll et al. 2009; Emmanoulopoulos et al. 2011; 
Patrick et al. 2011a, 2011b), here a simultaneous analysis of 
both the 2007 and 2010 observations provides a good fit to 
the data simply allowing the normalisation of the compo- 
nents to vary. In line with each of the previous analyses, we 
find an emissivity q = 2.9l '4 although a some what more 
poorly constrained spin parameter a = 0.60l o 'g| i- e - consis- 
tent with zero spin in this analysis of both Fairall 9 Suzaku 
observations 



A2 Mrk 205 

The best-fit model for Mrk 205 presented here features par- 
tial covering to account for additional flux at X-ray ener- 
gies > lOkeV. The subsequent addition of a broad Gaus- 
sian or relline component has little affect upon the overall 
fit (an improvement of A\ 2 ~ 4) due to the spectral cur- 
vature introduced as a result of a partial coverer with col- 
umn density iV H = (S-ljl^s) x 1q2S cm~ 2 , ionisation log(£) = 
2.8lg'5 erg cms -1 and covering fraction Cf rac = 13%. How- 
ever, this AGN can be equally well described with a very 
broad iron line with equivalent width EW = 254lgg eV and 
an increase of A^ 2 ~ +23 when the relline component 
is removed and the model refit without a partial coverer 
in place. This yields a relatively high emissivity index of 
q = 3.4lJ'g, an inclination of i° — 30lg° and a lower limit 
placed upon the spin parameter a > 0.1. 

A3 NGC 3227 

The six observations of NGC 3227 amount to nearly 500 000 
2-10 keV counts and a relatively simple model is formed to 
describe each component and the variations between ob- 
servations. The best-fit model for NGC 3227 consists of a 
high-column partial covering component of column A^h = 
(2.9±°;^) x 10 23 cm with a variable covering fraction of 
Cfrac ~ 32 — 89% in addition to a single warm absorbing 
zone and both hard and soft excesses fitted with reflionx 
and COMPTT respectively. Examining the Fe K region there 
are clear absorption lines due to FeXXV and FeXXVI and 
small residuals indicative of weak to moderate emission from 
the inner regions of the accretion disc (with EW generally 
higher in a lower continuum flux spectrum, also noted by 
Fukazawa et al. 2011). A good fit to each of the six ob- 
servations is found simply by letting the strength of the re- 
flionx and COMPTT to vary as well as the covering fraction 



of the partial coverer. All absorber properties such as ioniza- 
tion and column density are tied between observations other 
than the covering fraction and ionization of the partial cov- 
erer which is required to vary between log£ ~ 0.50 — 1.53 
with each observation. Using the relline model to account 
for the broad residuals with parameters tied (but again al- 
lowing normalisation to vary) yields an improved fit to the 
data with xl = 4465.6/4189. The relline model yields a 
relatively typical emissivity of q = 2.71q'| and an inclination 
of i = 33±1 ° and R in = 11±| R g (Tableland Figure[7j). 

A4 NGC 3516 

The two observations of NGC 3516 differ in their shape sig- 
nificantly (Markowitz et al. 2008), the 2005 Suzaku obser- 
vation 2-10 keV flux is a factor of 1.8 higher than in the 
2009 observation. Here we have formed a model in which a 
good fit is found to both observations, allowing for changes 
in the absorber/ partial coverer properties and strength of 
reflection between data sets, although the basic geometry 
is maintained and a good overall fit is found. As noted in 
Patrick et al. (2011b) and Turner et al. (2011), the 2009 ob- 
servation showed no strong indication of a broad red-wing, 
however, the Markowitz et al. (2008) analysis of the earlier 
2005 data do suggest the presence of such a feature. The 
2005 observation is more absorbed than the 2009 observation 
(i.e. C*f rac = 83% compared to 18% in the 2009 observation), 
yet it has a much higher continuum flux level throughout 
the spectrum above ~ 2 keV (Figure [7} . This simultaneous 
analysis yields consistent results with previous work, retriev- 
ing the broad feature found by Markowitz et al. (2008) in 
the 2005 data. Fitting a relline component to both data 
sets with all parameters (other than normalisation) tied sug- 
gest a fairly typical and relatively weak broad line with 
EW = 58l^eV and significantly weaker EW = 14^ 2 eV 
in the 2005 and 2009 observations respectively. An upper 
limit to the spin parameter is found a < 0.30 and inclina- 
tion i < 41° at a disc emissivity of q = 3.11q2- 



A5 NGC 4051 

The three Suzaku observations of NGC 4051 (as used in Lob- 
ban et al. 2011) included in this analysis vary significantly 
(see Figure [S] and 0. In particular, in the 2005 observa- 
tion (OBSID: 700004010) the source dips into an extended 
period of low flux versus the two 2008 observations which 
show the object in a period of high flux (-F2-10 kevzoos ~ 
3.8 x F 2 _io fcei/2005, see Figure [7] Table [B5)l . Similarly to 
Lobban et al. (2011), we can successfully describe each ob- 
servation and the long term spectral variability with a par- 
tial covering scenario whereby parameters such as the col- 
umn density and ionization of the absorbing zones of gas 
remain approximately constant (Nn ~ 9 x 10 22 cm~ 2 , over 
these time scales at least). The differences in the broad-band 
spectrum of each observation can be accounted for simply 
by allowing the normalisation of the intrinsic powerlaw, soft 
excess and distant reflection component to vary, in addition 
to the covering fraction of the partially covering absorption 
zone. Here we find covering fractions of 67%, 11% and 27% 
for observations 1, 2 and 3 respectively albeit with a much 
flatter intrinsic powerlaw (r ~ 1.88 here versus F ~ 2.49 in 
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Lobban et al. 2011), however, this is likely due to the treat- 
ment of the soft excess and reflection components. Here (for 
consistency with the analysis of other objects in the sample) 
we account for the soft excess with a COMPTT component 
which is akin to a second soft powerlaw, this could in part 
explain the discrepancies between the two different intrinsic 
powerlaw components. 

Similarly to Lobban et al. (2011), we find strong 
evidence for blue-shifted absorption in the FeK region 
indicative of an outflowing highly ionized zone of gas. 
When described as such, we find an outflow velocity of 
58001 

Kioo km s 1 in each observation, consistent with Lob- 
ban et al. (2011). Both the Lobban et al. (2011) and the 
baseline models in this analysis describe the spectral vari- 
ability with absorption dominated models. Replacing the 
partial coverer with a reflection component representing 
emission from the inner regions of the accretion disc as an 
alternate means of accounting for the hard excess and spec- 
tral curvature (i.e. a dual reflector) yields a fit worse by 
Ax 2 ~ +137 for 4 fewer de grecs of freedom. Parameters such 
as SMBH spin and the emissivity index of the disc are forced 
to high values i.e. a ~ 0.996 and q = 6.ltg"J. This is in order 
to smooth the reflection continuum to the extent to which it 
can successfully account for the long term spectral variabil- 
ity between observations i.e. accounting for both the hard 
excess and subtle continuum curvature changes. The partial 
covering scenario is therefore preferred both statistically and 
in terms of physical implications i.e. extreme parameters and 
relativistic blurring are not required which would otherwise 
significantly deviate from the sample norm. Note that there 
is still a evidence for an highly ionized outflow in the FeK 
region, regardless of the application of a dual reflector or 
partial covering based model. 

Based upon this, it is clear that the long term spec- 
tral variability of NGC 4051 cannot be produced purely by 
varying the reflected flux in a reflection dominated model; 
instead a partial covering scenario must be invoked to some 
extent in order to accurately reproduce the differences be- 
tween each of the three Suzaku observations. Reintroducing 
a partial covering geometry (in addition to a dual reflec- 
tor, Section 3.4) restores a fit similar to the baseline model. 
It is, therefore, perhaps more feasible that the broad-band 
spectrum and variability of NGC 4051 (in addition to other 
AGN) is primarily a result of variations in covering fraction 
in an absorption and not reflection dominated spectrum. 
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Table Bl. Summary of observations for the objects in the sample. 1 The observed 2— lOkeV flux for XIS, 15-50 keV flux for HXD and 
20-100 keV flux for BAT, in units 10 — 11 ergcm — 2 s — 1 from the baseline model. The XIS count rates listed are per XIS. We use BAT data 
from the 58 month BAT catalogue (Baumgartncr ct al. 2010). 



UDJCCt 


Mission 


Instrument 


uatc 


Exposure (s) 


Count rate 


r lux 


LJDS. LL) 


No. Counts 




SuzaJcu 


XIS 
HXD 


2007-07-25 


205863 
142600 


1.251 ±0.002 
0.052 ± 0.002 


1.75 
2.93 


702041010 


519055 
45303 


1x1 U41y— O ( ( 


SuzaJcu 


VTQ 

HXD 


2010-01-16 


IZZoou 

104900 


n QKfi -U n no 
0.028 ±0.001 


1 /in 
2.18 


704064010 


oi i a a n 

ZL I 44U 

22240 




OWlit 


"R AT 






(o a -\- n q*\ \s in - 4 
yZ.Q ± yj.o) X 1U 


"1 7 1 
1. I 1 




y4u 


3C 111 


SuzaJcu 


XIS 
HXD 


2008-08-22 


122378 
101900 


0.610 ±0.002 
0.068 ±0.002 


1.95 
3.29 


703034010 


152164 
36150 




Swift 


BAT 






(10.2 ±0.3) x 10~ 4 


7.43 








SuzaJcu 


XIS 
HXD 


2006-02-09 


41932 
31870 


q 082 + OO'i 
0.099 ± 0.004 


4.63 
6.24 


700001010 


389949 
13996 




SuzaJcu 


XIS 
HXD 


2006-02-16 


42555 
34540 


2.464 ± 0.004 
0.144 ±0.005 


4.01 
5.72 


700001020 


309300 
21832 


3C 120 


Suzaku 


XIS 
HXD 


2006-02-23 


40907 
36200 


2.53 ± 0.004 
0.110 ±0.004 


4.04 
6.07 


700001030 


312550 
15360 




SuzaJcu 


HXD 


2006-03-02 


A none; 
37870 


o -U n nn/i 
Z.oOl ± U.UU4 

0.093 ±0.003 


o.yo 
5.59 


700001040 


OGn^fin 
zyuooU 

14533 




Swift 


BAT 


— 




(8.5 ± 0.2) x 10~ 4 


5.82 






3C 382 


SuzaJcu 


AID 

HXD 


2007-04-27 


loUooU 
114300 


O KVkK -I- n flflQ 

0.097 ±0.002 


4.UD 

5.11 


702125010 


OooZoo 

47714 




Swift 


BAT 


- 




(7.6 ± 0.2) x 10~ 4 


5.11 






3C 390.3 


SuzaJcu 


AID 

HXD 


2006-12-14 


1 / yoUU 
92060 


i qqi _i_ n nnQ 
l.ool it U.uUo 

0.096 ±0.003 


q in 
o. 1U 

5.81 


701060010 


oooo i o 
37078 




Swift 


BAT 


— 




(9.09 ± 0.18) x 10~ 4 


6.82 






3C 445 


SuzaJcu 


AID 

HXD 


2007-05-25 


loy / oy 
109500 


n t c:q -U n nm 
U.loy it U.uUl 

0.049 ± 0.002 


n in 
2.69 


702056010 


/i anos 
4oUyo 

37749 




Swift 


BAT 


— 




(3.7 ±0.2) x 10~ 4 


2.35 






4C 74.26 


SuzaJcu 


XIS 
HXD 


2007-10-28 


91583 
87340 


1.493 ±0.003 
0.084 ± 0.002 


3.13 
4.00 


702057010 


276179 
34810 




Swift 


BAT 






(4.6 ± 0.2) x 10~ 4 


3.37 






A 1 on 
Ark IzU 


SuzaJcu 


XIS 

T4"VTl 


2007/04/01 


100864 
oy4 / U 


1.896 ± 0.003 

U.114 ± U.UUo 


3.05 

O.40 


702014010 


384821 




Swift 


BAT 


_ 


2453000 


(6.9 ± 0.4) x 10~ 4 


4.89 




1690 


Ark 546 


SuzaJcu 


XIS 
HXD 


2007-06-26 


99978 
81330 


2.277 ±0.003 
0.021 ±0.002 


1.84 
1.54 


702117010 


458819 
15343 




SuzaJcu 


XIS 
HXD 


2007/06/07 


167814 
127310 


1.718 ±0.002 
0.089 ± 0.002 


2.32 
2.97 


702043010 


581331 
46809 


Fairall 9 


SuzaJcu 


XIS 
HXD 


2010/08/26 


229296 
162200 


1.780 ± 0.002 
0.039 ±0.002 


2.17 
3.51 


705063010 


742007 
65230 




Swift 


BAT 






(4.7 ± 0.2) x 10~ 4 


3.34 




1619 
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Table Bl - continued Summary of observations for the objects in the sample. 1 The observed 2—10 keV flux for XIS, 15-50 kcV flux for 
HXD and 20-100 keV flux for BAT, in units 10 -11 ergcm -2 s _1 from the baseline model. The XIS count rates listed are per XIS. We use 
BAT data from the 58 month BAT catalogue (Baumgartncr ct al. 2010). 



Object 


^Mission 


Instrument 


Date 


Exposure (s) 


Count rate 




Flux 1 


Obs. ID 


No. Counts 




Suzaku 


XIS 
HXD 


finfi7 no ni 
ZUU / -Uo-Ul 


25453 
20050 


5.212 ±0.010 
0.289 ± 0.006 




10.64 
15.96 


I UzlloUlU 


266796 
11647 




Suzaku 


XIS 
HXD 


2007-08-06 


30623 
24080 


/ ' 1 OO 1 / \ ill / 1 

6.138 ± 0.010 
0.340 ± 0.005 




12.58 
17.97 


702113020 


377880 
15349 


IC 4329A 


Suzaku 


XIS 
HXD 


9007 OS 1 1 


26896 
22110 


5.554 ±0.010 
0.345 ± 0.006 




11.89 
18.77 




300624 
15541 




Suzaku 


XIS 
HXD 


2007-08-16 


24219 
18750 


5.261 ±0.010 
0.300 ± 0.006 




10.88 
16.47 


702113040 


256495 
11923 




Suzaku 


XIS 
HXD 


2007-08-20 


24026 
17520 


3.221 ±0.008 
0.242 ±0.006 




7.03 
13.69 


702113050 


156141 
9755 




Swift 


BAT 


- 


- 


(26.8 ±0.3) x 10 


-4 


19.02 






IRAS 13224-3809 


Suzaku 


XIS 
HXD 


2007-01-26 


197938 
158500 


0.085 ±0.001 
0.004 ± 0.002 




0.06 
0.015 


701003010 


39148 
86214 


MCG-02-14-009 


Suzaku 


XIS 
HXD 


2008/08/28 


142152 
120028 


0.216 ± 0.001 
0.017 ±0.002 




0.43 
0.61 


703060010 


63436 
40150 


MCG-02-58-22 


Suzaku 


XIS 
HXD 


2009-12-02 


138969 
97980 


2.921 ± 0.003 
0.163 ±0.002 




4.87 
8.73 


704032010 


817420 
44009 




Swift 


BAT 






(9.8 ± 0.2) x 10" 




7.48 






MCG-05-23-16 


Suzaku 


XIS 
HXD 


2005-12-07 


95677 
79690 


2 822 + 003 
0.313 ±0.003 




8.93 
14.52 


700002010 


815975 
55913 




Swift 


BAT 






(19.8 ±0.2) x 10 


-4 


14.27 




4526 




Suzaku 


XIS 
HXD 


2006/01/09 


143196 
118900 


2.862 ± 0.003 
0.094 ± 0.002 




4.38 
4.70 


700007010 


1240090 
53689 


MCG-06-30-15 


Suzaku 


XIS 
HXD 


2006/01/23 


98483 
76800 


2.461 ±0.003 
0.103 ±0.003 




3.81 
5.06 


700007020 


734189 
36695 




Suzaku 


XIS 
HXD 


2006/01/27 


96691 
83660 


2.708 ± 0.003 
0.104 ± 0.002 




4.16 
4.98 


700007030 


792389 
40408 




Swift 


BAT 


- 




(6.8 ±0.3) x 10~ 


-4 


3.86 






MCG+8-11-11 


Suzaku 


XIS 
HXD 


2007-09-17 


98748 
82900 


2.863 ± 0.004 
0.196 ± 0.003 




6.48 
10.44 


702112010 


568883 
43023 




Swift 


BAT 


- 




(8.8 ±0.2) x 10~ 


-4 


6.27 




2082 


MR 2251-178 


Suzaku 


XIS 
HXD 


2009-05-07 


136924 
103800 


2.089 ± 0.003 
0.101 ± 0.002 




4.23 
5.47 


704055010 


576608 
40562 




Swift 


BAT 


- 




(9.1 ±0.2) x 10- 


-4 


6.02 






Mrk 79 


Suzaku 


XIS 
HXD 


2007-04-03 


83704 
76920 


0.768 ± 0.002 
0.027 ± 0.002 




1.47 
1.79 


702044010 


130639 
19777 




Swift 


BAT 


- 




(4.1 ±0.2) x 10- 


-4 


2.87 






Mrk 110 


Suzaku 


XIS 
HXD 


2007-11-02 


90871 
80370 


1.220 ± 0.003 
0.053 ± 0.002 




2.13 
2.99 


702124010 


224098 
26958 




Swift 


BAT 


- 




(5.4 ±0.3) x 10- 


-4 


3.95 




2128 


Mrk 205 


Suzaku 


XIS 
HXD 


2010-05-22 


100961 
85280 


0.501 ± 0.002 
0.019 ± 0.001 




0.93 
1.94 


705062010 


103464 
13647 




Swift 


BAT 


- 




(1.4 ±0.2) x 10- 


-4 


1.05 






Mrk 279 


Suzaku 


XIS 
HXD 


2009-05-14 


160351 
139800 


0.263 ± 0.001 
0.024 ± 0.001 




0.49 
2.42 


704031010 


88122 
29520 




Swift 


BAT 






(5.2 ±0.3) x 10- 


-4 


3.93 




2617 


Mrk 335 


Suzaku 


XIS 
HXD 


2006/08/21 


151296 
131744 


1.324 ±0.002 
0.012 ±0.001 




1.49 
1.31 




606927 




Swift 


BAT 






(2.5 ±0.3) x 10- 


-4 


1.81 




830 


Mrk 359 


Suzaku 


XIS 
HXD 


2007-02-06 


107507 
96130 


0.329 ±0.001 
0.007 ±0.001 




0.52 
0.98 


701082010 


73004 
13136 




Swift 


BAT 






(1.2 ±0.1) x 10- 


-4 


0.92 
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Table Bl - continued Summary of observations for the objects in the sample. 1 The observed 2-10 keV flux for XIS, 15-50 kcV flux for 
HXD and 20-100 keV flux for BAT, in units 10 -11 ergcm -2 s _1 from the baseline model. The XIS count rates listed are per XIS. We use 
BAT data from the 58 month BAT catalogue (Baumgartncr ct al. 2010). 



Obiect 


IVTi on 


Tn ^tn l m OTit 


Date 


FiVTifwnrp ( ^ 1 

1 J A KJ J Lll K-r V / 


(~^OH Tit" l~£ltp 




Flux 1 


Obs. ID 


lMn nnnnts 

1 ^ \J m V_ yJ LillLrO 




SuzaJcu 


XIS 
HXD 


2006-04-25 


24576 
14510 


2.984 ± 0.006 
0.143 ±0.006 




4.68 
7.50 


701093010 


221506 
6647 




SuzaJcu 


XIS 
HXD 


2006-10-14 


25930 
21880 


3.802 ± 0.007 
0.128 ±0.005 




4.98 
7.89 


701093020 


297900 
9680 


Mrk 509 


SuzaJcu 


XIS 
HXD 


2006-11-15 


24447 
17340 


3.792 ± 0.009 
0.118 ±0.005 




4.74 
7.57 


701093030 


186702 
7260 




SuzaJcu 


XIS 
HXD 


2006-11-27 


33094 
27560 


3.278 ±0.007 
0.110 ±0.004 




4.73 
7.19 


701093040 


218248 
11925 




Swift 


BAT 






(8.4 ±0.2) x 10" 




6.53 








SuzaJcu 


XIS 
HXD 


2006-11-16 


97869 
90500 


927 + 002 
0.025 ±0.002 




1.32 
1.48 


701035010 


1 83736 

i *- ) t • } y.t 

26879 


Mrk 766 


SuzaJcu 


XIS 
HXD 


2007-11-17 


59364 
47660 


0.685 ±0.002 
0.024 ± 0.002 




1.36 
1.67 


701035020 


82749 
8974 




Swift 


BAT 






(2.5 ± 0.2) x 10" 




1.66 








SuzaJcu 


XIS 
HXD 


2007-01-22 


51753 
40710 


819 + 003 
0.039 ± 0.003 




1.40 
2.42 


701084010 


86781 
11274 


Mrk 841 

1V1I.JV UtI 


SuzaJcu 


XIS 
HXD 


2007-07-23 


50925 
44370 


830 + 003 
0.054 ± 0.003 




1.44 
2.82 


701084020 


86626 
13282 




Swift 


BAT 






(3.2 ±0.2) x 10" 


-4 


2.29 








SuzaJcu 


XIS 
HXD 


2008-01-21 


1 60506 
136600 


326 + 001 
0.084 ± 0.002 




1.28 
4.72 


702047010 


107383 

1U ' •)'*>•) 

55436 


NGC 1365 


SuzaJcu 


XIS 
HXD 


2010-06-27 


151613 
114300 


1 59 + 001 
0.070 ±0.002 




0.61 
4.33 


705031010 


51909 
38083 




SuzaJcu 


XIS 
HXD 


2010-07-15 


302175 
231500 


1 16 + 001 
0.060 ±0.001 




0.39 
3.49 


705031020 


77717 

84887 




Swift 


BAT 






(5.2 ±0.1) x 10" 


-4 


5.11 








S uzajcu 


XIS 
HXD 




37503 
29900 


0.372 ±0.002 
0.020 ± 0.003 




1.04 
1.71 


/ UUUUJU1U 


43167 
7137 


NGC 2992 


SuzaJcu 


XIS 
HXD 


2005-11-19 


37494 
31890 


0.486 ± 0.002 
0.042 ± 0.003 




1.37 
2.23 


700005020 


56041 
10370 




SuzaJcu 


XIS 
HXD 


2005-12-13 


46836 
41470 


0.408 ± 0.002 
0.047 ± 0.003 




1.15 
2.44 


700005030 


59084 
15070 




Swift 


BAT 


- 




(2.6 ±0.2) x 10" 


-4 


1.84 






NGC 3147 


SuzaJcu 


XIS 
HXD 


2010-06-03 


150048 
122800 


0.093 ±0.001 
0.009 ± 0.001 




0.17 
0.29 


705054010 


31054 
21785 




Suzaku 


XIS 
HXD 


ZUUo-lu-Zo 


58917 
48070 


1.475 ±0.004 
0.145 ±0.003 




3.91 
5.58 


/ UoUZzUlU 


175651 
21550 




Suzaku 


XIS 
HXD 




53700 
46740 


0.496 ± 0.002 
0.123 ±0.003 




1.85 
7.26 


/ UoUzzUzU 


54717 
18735 




Suzaku 


XIS 
HXD 


2008-11-12 


56572 
46680 


0.707 ± 0.003 
0.117 ±0.003 




2.58 
6.78 


703022030 


81449 
18919 


NGC 3227 


Suzaku 


XIS 
HXD 


2008-11-20 


64568 
43430 


278 + 002 
0.083 ± 0.003 




1.01 
4.82 


703022040 


37518 
17263 




Suzaku 


XIS 
HXD 


2008-11-27 


79430 
37420 


n 4- n nri2 

W.wUO _1_ U.UUii 

0.104 ±0.003 




2.16 
6.26 


703022050 


92162 
14795 




Suzaku 


XIS 
HXD 


2008-12-02 


51411 
36910 


0.413 ±0.002 
0.082 ± 0.003 




1.61 
5.70 


703022060 


43570 
13177 




Swift 


BAT 






(10.3 ± 0.2) x 10 


-4 


7.68 








Suzaku 


XIS 
HXD 


2005/10/12 


134469 
115400 


0.681 ±0.001 
0.125 ±0.002 




2.37 
6.68 


100031010 


289936 
49782 


NGC 3516 


Suzaku 


XIS 
HXD 


2009/10/28 


251356 
178200 


0.456 ±0.001 
0.059 ±0.001 




1.35 
3.32 


704062010 


235894 
58222 




Swift 


BAT 






(10.5 ± 0.2) x 10 


-4 


7.76 
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Table Bl - continued Summary of observations for the objects in the sample. 1 The observed 2—10 keV flux for XIS, 15-50 kcV flux for 
HXD and 20-100 keV flux for BAT, in units 10 -11 ergcm -2 s _1 from the baseline model. The XIS count rates listed are per XIS. We use 
BAT data from the 58 month BAT catalogue (Baumgartncr ct al. 2010). 



\J DJ CCt 


Mission 


Instrument 


uatc 


Exposure (s) 


Count rate 




r lux 


<^/DS. lij 


No. Counts 




SuzaJcu 


XIS 
HXD 


2006/06/24 


75719 
63930 


2.066 ± 0.003 
0.167 ±0.003 




4.59 
9.26 


701033010 


416485 
33592 


NGC 3783 


SuzaJcu 


XIS 
HXD 


2009/07/10 


209503 
162000 


1.942 ±0.002 
0.235 ±0.002 




5.92 
12.04 


704063010 


805999 
87319 




Swift 


BAT 






(16.6 ± 0.2) x 10 


—4 


11.54 








SuzaJcu 


XIS 
HXD 


2005/11/10 


119578 
112600 


0.406 ±0.001 
0.037 ±0.002 




0.87 
1.63 


700004010 


149284 
29700 




SuzaJcu 


XIS 


2008/11/06 


274350 


1.858 ± 0.002 




2.46 


703023020 


1036060 


NGC 4051 


HXD 




204500 


0.062 ±0.001 




3.05 


75619 




SuzaJcu 


HXD 


2008/11/23 


7Qqor; 
I oooO 

58530 


1 r ) r 7fi _i_ n nrtQ 
L.Z i it U.UUo 

0.048 ± 0.002 




1 7Q 

i. ( y 
2.54 


703023010 


21296 




Swift 


BAT 


- 




(4.0 ±0.2) x 10~ 


■4 


2.80 






NGC 4151 


SuzaJcu 


.A. IS 
HXD 


2006-12-18 


iz^you 
123500 


i nQfi _i_ n nno 

i.uoD it u.uuz 
0.293 ± 0.002 




,1 or: 
4.O0 

15.90 


701034010 


ZDZool 

78564 




owirt 


TD AT 






//1Q ft _L n o\ \y in 
(4o.O ± U.z J X 1U 


-4 


OZ.OU 






NGC 4593 


SuzaJcu 


XIS 
HXD 


2007-12-15 


118842 
101600 


0.524 ±0.002 
0.033 ± 0.002 




1.04 
2.15 


702040010 


127503 
30049 




Swift 


BAT 






(7.3 ±0.2) x 10~ 


-4 


5.74 








Suzaku 


XIS 
HXD 


2006-08-08 


47753 
38550 


3.064 ± 0.005 
0.324 ±0.004 




10.05 
17.30 


701030010 


441844 
25416 


NGC 5506 


Suzaku 


XIS 
HXD 


2006-08-11 


53296 
44830 


3.230 ±0.005 
0.343 ±0.004 




10.58 
18.40 


701030020 


519889 
30625 




Suzaku 


XIS 
HXD 


2007-01-31 


57406 
44670 


2.854 ±0.005 
0.321 ±0.004 




9.90 
17.11 


701030030 


330141 
29399 




Swift 


BAT 






(23.2 ± 0.2) x 10 


-4 


15.95 








Suzaku 


XIS 
HXD 


2007-06-18 


31119 
25590 


0.369 ± 0.003 
0.046 ± 0.004 




0.81 
2.70 


702042010 


23857 
9397 




Suzaku 


XIS 
HXD 


2007-06-24 


35915 
31210 


0.650 ±0.003 
0.050 ± 0.003 




1.34 
3.12 


702042020 


47786 
9880 




Suzaku 


XIS 
HXD 


2007-07-08 


30700 
26960 


1.349 ±0.005 
0.064 ± 0.004 




2.54 
4.34 


702042040 


83863 
9824 


NGC 5548 


Suzaku 


XIS 
HXD 


2007-07-15 


30020 
24470 


834 + 004 
0.062 ± 0.004 




1.74 
4.16 


702042050 


51085 
7612 




Suzaku 


XIS 
HXD 


2007-07-22 


28918 
23060 


1.656 ± 0.005 
0.095 ± 0.004 




3.28 
0.22 


702042060 


96936 
8336 




Suzaku 


XIS 
HXD 


2007-07-29 


31810 
27610 


1.075 ± 0.004 
0.075 ± 0.004 




2.19 
4.71 


702042070 


69507 
9377 




Suzaku 


XIS 
HXD 


2007-08-05 


38776 
30380 


0.558 ± 0.003 
0.054 ± 0.003 




1.19 
3.41 


702042080 


44486 
9958 




Swift 


BAT 






(6.1 ±0.2) x 10~ 


■4 


4.56 






NGC 7213 


Suzaku 


XIS 
HXD 


2006-10-22 


90736 
84290 


1.528 ±0.002 
0.064 ± 0.002 




2.41 
3.48 


701029010 


419604 
33910 




Swift 


BAT 






(3.6 ±0.2) x 10- 


■4 


5.13 
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Table Bl - continued Summary of observations for the objects in the sample. 1 The observed 2-10 keV flux for XIS, 15-50 kcV flux for 
HXD and 20-100 keV flux for BAT, in units 10 -11 ergcm -2 s _1 from the baseline model. The XIS count rates listed are per XIS. We use 
BAT data from the 58 month BAT catalogue (Baumgartncr ct al. 2010). 



Object 


Mission 


Instrument 


Date 


Exposure (s) 


Count rate 




Flux 1 


Obs. ID 


No. Counts 


NGC 7314 


Suzaku 


XIS 
HXD 


2007-04-25 


109020 
85780 


0.296 ±0.001 
0.026 ± 0.002 




0.88 
1.43 


702015010 


67434 
25675 




Swift 


BAT 






(4.6 ± 0.2) x 10" 


-4 


3.62 






NGC 7469 


Suzaku 


XIS 
HXD 


2008/06/24 


112113 

85315 


1.091 ±0.002 
0.068 ± 0.002 




2.10 
3.22 


703028010 


248180 
29831 




O Wll l 


BAT 








-4 


A 87 




91 7Q 


PDS 456 


Sii7Pikn 


XIS 
HXD 


2007-02-24 


190600 
64150 


0.237 ±0.001 
0.012 ±0.003 




0.35 
0.25 


701056010 


92499 
20372 


PG 1211+143 


Suzaku 


XIS 
HXD 


2005-11-24 


96324 
78800 


0.268 ±0.001 
0.004 ±0.002 




0.39 
0.50 


700009010 


81584 
29181 


RBS 1124 


Suzaku 


XIS 
HXD 


2007-04-14 


86228 
82970 


0.260 ±0.001 
0.016 ±0.001 




0.49 
0.91 


702114010 


46699 
14206 


SWIFT J2127.4+5654 


Suzaku 


XIS 
HXD 


2007/12/09 


91730 
83321 


1.373 ±0.004 
0.074 ± 0.002 




3.35 
3.33 


702122010 


127435 
27970 




Swift 


BAT 






(4.3 ± 0.3) x 10- 


-4 


3.62 




1689 


TON S180 


Suzaku 


XIS 
HXD 


2006-12-09 


120661 
102400 


0.701 ±0.002 
0.012 ±0.002 




0.56 
0.77 


701021010 


172889 
27241 
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Table B2. Components for the baseline model (without broad or diskline emission) to the observations with Suzaku XIS, HXD and BAT data from Swift. a Unabsorbed POWERLAW 
2 normalisation given in units (10~ 2 phkeV -1 cm~ 2 s _1 ). 6 Flux for COMPTT quoted over the 0.6-10.0 kcV range and REFLIONX over the 2.0-100.0keV range in units 10 — 1:L ergcm — 2 



> 
en 



Ionization parameter given in units erg cms 



Column density measured in units 10 cm 



frozen parameter, f indicates parameters arc tied during the analysis of multiple observations. 



^turb 



1000 km s 1 . f only the intrinsic powcrlaw is absorbed. * indicates a 



•x 

> 
m 

o 
o 
o 



Object 



Powerlaw 

r 



Norm™ 



COMPTT 
kT (kcV) 



Flux 6 



REFLIONX 

Z Fc £ c 



Flux 6 



Warm absorber 
N d 



BAT/XIS 



xl 



1H 0419-577 Obs 
1,2 

1H 0419-577 



1 70+ 001 
i -' u -o.oi 



1.70f 



4C 74.26 
Ark 120 
Ark 564 



1.92 
1.90 
2.34 



+0.01 
-0.01 
0.01 
0.04 
+0.01 
0.01 



3,4,5 



IRAS 13224-3809 2.50 



+0.20 
0.16 



MCG-02- 14-009 1 

MCG-02-58-22 1.64 

MCG-05-23-16 1.84 

MCG-06-30-15 1.87 
Obs 1,2,3 



0.39 



+0.01 
0.01 



0.31 



+0.01 
-0.01 



3C 


111 


1 


5S +o.oi 

oo -0.01 


0.44 


3C 


120 Obs 1,2,3,4 


1 


65+ ' 01 
DO -0.01 


0.98 


3C 


120 Obs 4 




1.65f 


0.86 


3C 


382 


1 


70+0.01 

■' 8 -0.01 


1.16 


3C 


390.3 


1 


70+0.02 

• ,,:> -0.02 


0.75 


3C 


445 


1 


43+0.01 

*°-0.01 


0.03 



+0.01 

-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.02 
+0.01 
-0.01 



1.06 
0.95 
1.10 



+0.01 
-0.01 

0.06 
0.06 
+0.01 
0.01 



Fairall 9 Obs 1 


1 01+0.01 

i - 8i -0.01 





61 


Fairall 9 Obs 2 


1.81f 


0. 


56 


IC 4329A Obs 1,2, 


1 91+0.01 


3 


.51 



+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 



< 0.02 



7 r+0.4 
' -°-0.4 



7.5f 



6 9+°' 3 
6.9f 

< 56.5 

< 55.9 



< 8.5 
4 7+27.3 
^■'-2.4 



9.2 
9.2f 



+0.7 
0.6 



< 14.8 



1.4 



+0.1 
0.1 



1.4t 



1.81 
l.Oj 



1.4 
1.0 
< 25 



< 2.4 

< 2.0 



1.1 



0.1 
0.1 
l.lf 



< 1.5 



0.59 



+0.01 

0.01 



0.46 



+0.01 

0.01 



+0.03 
-0.03 
+0.01 
0.01 
+0.36 
0.31 
< 0.429 



1.82 
0.97 
0.35 



1.05 
3.29 



+0.05 
0.05 
+3.22 
2.53 



0.42 

3.87 



0.33 



1.0* 



l.Of 



1 c+0.3 
i -°-0.2 

1.5f 
1.0* 
5+° 1 
1.0* 



,+0.01 


2 2+ 1 ' 3 
z - -0.4 


'-0.01 


.+0.03 
-0.03 


2.2t 




6+ 01 


,+0.01 




'-0.01 





10.1 



+ 1.9 
3.0 



10. If 



< 19.4 

< 19.4t 
< 2.1 
2 ' 1 

9.9 



1 

+4.5 
7.8 



14.8 



+7.7 
-9.0 
14.8f 

< 1.1 



0.36 



+0.21 

0.14 



0.45 



+0.28 
0.17 



1.92 
2.34 
1.67 
2.37 
2.79 



+0.30 
-0.28 
+0.35 
0.30 
+ 0.14 
-0.93 
+2.44 
■0.65 
10.14 
1.03 



1.0* 


< 5.2 


1 


.23 


1 4+0-5 


< 22 





.93 


1.0* 


< 7.8 





.12 



+ 1.36 
-0.81 
+0.44 
-0.53 
+0.19 
-0.09 



7.94 

7.13 
9.15 



+ 1.83 
-0.61 
+5.47 
3.39 
+ G.18 
1.31 



+0.01 
'-0.04 


12+ 004 
u - ± -0.03 












q O+0.7 


< 21.0 





41+0.09 
'^ 1 -0.17 


+0.01 
-0.01 


94+0-01 
u - 9 ^-o.oi 




1.5: 


-1.4 
-1.5 


1.61+1 


89 
41 


1.0* 


56.6+J-j 


2 


74 +0.47 
'^-0.49 


+0.01 
-0.01 


o 07+O.O3 

z -°'-0.03 


c q+0.3 


16.7 


+ 1.1 
-1.0 


0.45« 


02 

.02 


1.0* 


< 1.0 


2 


7fi +0.21 
■ '°-0.27 


■+0.04 
-0.04 


1 62+ 05 
1-O^_0.04 


< 8.8 


o.sl 


12.9 
0.1 


0.72+.° 


18 
10 


1.0* 


< 12 


1 


oo+O.ll 
ZO -0.35 



15.83±SS;|S 

178.0+™;° 

15.83f 

178.0f 
r, C1 +0.01 
°- 61 -o.oi 

09+ 01 
u - UM -o.oi 

0.09f 

U - JO -0.09 



3.31 



+0.21 e 
0.17 



18.94 



+0.78 



0.27 



0.75 
+0.03 
0.01 



+0.01 
0.01 



0.05 

< 0.02 e 

< 0.02 e 
> 461.0 



< 0.02 e 



+0.03 
0.04 
+0.02 
0.02 



0.35 
0.60 
< 0.0+ 
0.22 
5.00 



+0.01 
0.02 

+2.17 
2.56 



0.12 
1.46 
0.22 



+0.03 
0.02 
+0.01,/ 
-0.01 
0.06 
0.03 



n 47+0-09 

u - 4 ' -0.03 
341 n+ 24 ' 



> 2.87 

2 05+ ' 43 

> 2.87t 
2.05f 

Neutral 

Neutral 

Neutral 
2 74+0.09 
-0.12 

1 93+0.03 
1,M -0 .03 
Neutral 

1 00+ 010 
i.uu_ 07 

Neutral 
1 77+0.18 
i -' '-0.76 

2.34^1 



1.53 



-0.48 



1.84 



+0.05 
0.05 

< 0.01 
> 0.63 
Neutral 

< 0.52 



2.22^ 



-,+0.09 
-0.12 

Neutral 

+0.13 
-0.09 



0.76 

1.76 
2.43 



+0.05 
0.08 
+0.02 
0.05 



100% 

18% 
100% 
15% 
100% 
100% 
100% 
100% 

100% 
89% 
100% 

100% 
100% 
100% 
49% 



100% 

100% 
100% 
100% 
100% 
72% 

100% 
100% 
100% 

100% 
50% 



47+0 11 
u - 4 ' -0.11 



0.47f 



+0.14 

0.14 
+0.04 
0.04 



1.85 
0.83 

0.83f 
0.91 
0.93 
0.85 



+0.06 
-0.05 
+0.04 
0.04 
+0.12 
0.11 



0.72 
1.07 



+0.05 
-0.05 
+ 0.12 
0.12 



0.83Z 



0.63 
0.82 
0.95 



+0.03 
0.02 
+0.03 
0.03 
+0.06 
0.07 



1788.9/1750 

1101.1/1098 

3630.5/3452 

986.6/938 
1499.4/1486 
458.1/422 

1344.8/1303 

741.8/649 
1148.9/1024 



I-0.10 

J -0- 10 3604.6/3276 
0.83f 

°- 95 -0.02 2421.4/2200 



484.86/447 



623.6/543 
2023.7/1892 
1582.8/1466 
2026.8/1823 



a 
s 
o 

Co 

3 

S3 

o 
s 



3 
s 



s 



Co 

CO 



Table B2 - continued Components for the baseline model (without broad or diskline emission) to the observations with Suzaku XIS, HXD and BAT data from Swift. a Unabsorbcd 
POWERLAW normalisation given in units (10 -2 phkcV -1 cm~ 2 s _1 ). b Flux for COMPTT quoted over the 0.6-10.0 keV range and REFLIONX over the 2.0-100.0keV range in units 
10~ 1:L crgcm~ 2 s _1 . c Ionization parameter given in units ergcms -1 . d Column density measured in units 10 22 cm~ 2 . e v tulb = lOOOkms - 1 . ^ only the intrinsic powerlaw is 
absorbed. * indicates a frozen parameter, f indicates parameters are tied during the analysis of multiple observations. 



Object 



Powerlaw 

r 



Norm™ 



COMPTT 
kT (keV) 



Flux 6 



REFLIONX 

Z Fc £ c 



Flux 6 



Warm absorber 
N d 



log(5) c 



BAT/XIS 



MCG+8-11-11 

MR 2251-178 

Mrk 79 

Mrk 110 
Mrk 205 
Mrk 279 

Mrk 335 
Mrk 359 

Mrk 509 Obs 1,2, 
3,4 

Mrk 766 Obs 1 
Mrk 766 Obs 2 

Mrk 841 Obs 1,2 
NGC 1365 Obs 1 

NGC 1365 Obs 2 

NGC 1365 Obs 3 

NGC 2992 Obs 1, 
2,3 

NGC 3147 



1.80 



+0.01 
0.01 



1.56 



+0.01 
0.01 



1.55 



+0.02 
0.02 



i 71 +0.01 
^^-O.Ol 



1.97 
1.76 

2.00 
1.74 
1.66 



+0.02 
0.02 
+0.01 
0.01 



+0.02 
0.02 
+0.01 
0.01 
+0.01 
0.01 



1.99 
1.99f 



+0.01 
0.01 



1.79 
1.69 



+0.01 
0.01 
+0.01 
0.01 



1.69f 
1.69f 



1.58 



+0.01 
0.01 



-, 79 +0.02 
lw -0.03 



,+0.01 
J-O.Ol 



0.90 

0.32 

0.52 
0.31 
0.11 

0.51 
0.12 
0.89 



+0.01 
0.01 



+0.01 
0.01 



0.01 
-0.01 
+0.01 
-0.02 
+0.01 
-0.01 



-0.02 
0.01 
0.01 
+0.01 
0.01 



50+ 01 
u - au -o.oi 

< 0.05 



36+ 01 

03+ 01 
u - u,r> -0.01 



0.02 



+0.01 
0.01 



0.02 



+0.01 
0.01 



0.26 



+0.01 
0.01 



0.04 



+0.01 
0.01 



< 48.7 
10.3±i 2 4 



< 8.3 
< 151.5 
11.7 



+36.9 
5.3 



7.1 



+3.0 
2.3 

7.1f 



< 116.9 



< 16.6 



2.2 



+ 1.4 

l.i 



2.81°:? 

< 2.0 



< 2.2 

< 4.3 
1.2 



1.0 
1.0 



o.4l8:i 

0.4f 



< 3.9 



0.76 



+0.04 
0.10 



0.80 



+0.08 
0.08 



0.31 
0.08 



+0.17 
0.29 

+0.20 
0.03 



o.98l8;8 2 2 

o-o8±8:S? 

2-54l?; 2 8 | 



0.13 
0.15 



+0.02 
-0.01 
+0.02 
0.02 



o.3o; 



-1.44 
-0.23 



1.45 



+ 0.31 
1.14 



1.0* 



1.0* 
1.0* 
1.0* 



2.2 







1.1 

-0.4 

1.0* 
1.0* 



1.0* 
l.Of 

+0.2 
-0.2 

1.0* 
l.Of 
l.Of 
1.0* 
1.0* 



8.0 



+3.5 
2.4 



49.4ll 1 7 5 
< 1.5 
2.8 



+ 1.0 
0.5 



27.011° 
21.8l|f 

21 4+30.9 
zl - 4 -18.4 

o 9 +0.5 

3.2f 



3.9 
24.3 



+2.0 
-1.4 
+ 1.0 
0.9 



24.3t 
24.3t 



54.8 



+7.8 
8.3 



4q 0+23.4 



4.79l 2 -j 2 



0.42 
0.45 
1.14 

0.46 
0.52 
1.51 

0.50 
0.68 

1.43 
1.20 



+0.88 
0.13 
+0.06 
0.17 
+0.29 
-0.30 





14 


l°o 


14 


,+o. 


44 


'-0 


17 


+9. 


50 


-0 


94 


'18: 


11 




09 


'l°a 


14 




12 


.+0. 


86 


'-() 


51 




10 


>1° : 


09 


> + -°o 


11 




10 


i!8: 


06 




.08 


18 


29 




29 


+0. 


15 


-0 


.06 



n 7Q+0-08 

u -' 9 -0.07 



0.06 
0.89 
0.54 
0.06 
0.98 



+0.01 
0.01 
+0.26 
-0.20 
+0.03 
0.03 
+0.04 
0.01 
+0.44 
0.28 



50.94±&™ 

o.o8i8:8| 

> 452.53 



24+ ' 03 
u -^*-0 .02 



n O-.+0.03 
u ' ol -0.02 

2 41+ 010 

z '^ 1 -0.09 
+ 1.40 
-0.98 
+0.09 
-0.10 
+0.69 



5 

0.46 
2.71 



0.74 
+3.35 
3.00 
+ 1.98 
1.92 
+ 7.16 
-6.72 
6.18 
3.33 
+6.29 
15.77 

0.8418:81 



29.57 
36.79 
16.41 
65.38 
17.99 



n 4 o+0.03 
z -^°-0 .03 
Neutral 

9 SS+ ' 07 
A8S -0.08 

i - b8 -0.03 

,+0.13 
-0.14 



2.66Z 



7 tr+0.70 

1 on+0.30 
L -' ii '-0.51 
i 24 +0.96 



2.20 



+0.06 
0.07 



h0.06 



l-20l .07 

i -° i -0.03 
-+0.05 



2.94_ 

2.17: 



0.05 
0.10 



' -0.11 

Neutral 

o.8oi8:^ 

Neutral 
2 ik+0.18 
^• lo -0.39 

Neutral 

< 1.90 
Neutral 



100% 
100% 
100% 
100% 
100% 
100% 

13% 
100% 
72% 



100% 

100% 
100% 
> 92% 
100% 

96% 

96% 

95% 
100% 



0.50 



+0.04 
0.04 



1 04+ ' 06 
1 - u4 -0 .06 



1 40 +0 19 



1.05 
0.80 
1.15 

1.03 
0.72 



+ 0.11 
0.11 
+0.15 
0.15 
+ 0.11 
0.11 



+0.23 
0.23 
+0.20 
0.20 



n 67 +0.04 



0.98 
0.98f 



+0.12 
0.12 



71+0.09 
u - (± -0.08 
,+0.04 
-0.04 



0.76 1 



0.76f 
0.76f 



0.79 



+ 0.13 
0.13 



1049.9/938 

977.5/902 

601.7/545 

483.1/471 
167.7/165 
712.4/669 

842.9/723 
610.1/562 
1976.2/1870 

1057.1/1000 
937.6/857 

2172.5/1979 

1088.2/1079 
280.3/266 



Table B2 - continued Components for the baseline model (without broad or disklinc emission) to the observations with Suzaku XIS, HXD and BAT data from Swift. a Unabsorbed 
POWERLAW normalisation given in units (10 -2 phkcV -1 cm~ 2 s~ x ). b Flux for COMPTT quoted over the 0.6-10.0 keV range and REFLIONX over the 2.0-100.0keV range in units 
10- 



Ionization parameter given in units erg cms 



Column density measured in units 10 22 cm 2 . e v turh = 1000 km s 1 . f only the intrinsic powerlaw is 



absorbed. * indicates a frozen parameter. \ indicates parameters are tied during the analysis of multiple observations. 



Object 


Powerlaw 




COMPTT 






REFLIONX 






Warm absorber 










r 


Norm™ 


kT (keV) 


T 


Flux 6 






Flux 6 




lo g(?) 


Cfrac 


BAT/XIS 


xl 


NGC 3227 Obs 1 


1 Q2+ ' 01 


-, , 1+0.02 


A R+ ' 4 
4 - S -0.4 


1 2+ ' 1 


4 4 2 +0.08 


2+° 1 


n ,+0.1 
z - a -0.1 


1 ^ 77+O.88 
±J - 1 ' -0.89 


, ,7+0.11 
-0.07 

0.43t°°J 


2 nc+0.02 

Neutral 


100% 
100% 


0.82t°;°l 








13+ 001 






n qs+o-03 
u.ao_ 03 






13.274±°j 7 


29 24+ ' 46 


50+ - 18 
U.OU_ 27 


32% 






NGC 3227 Obs 2 


1.79f 


4.8f 


1.2f 


0.2f 


2.3f 


3.37f 
0.43f 
29.24+ 


2.05f 
Neutral 
, co+0.09 

1 - 0,3 -o.io 


100% 
100% 
81% 


0.821 




NGC 3227 Obs 3 


1.79f 


u - zo -o.oi 


4.8f 


1.2f 


76+ 03 


0.2f 


2.3f 


i, i 7 +0.59 


3.37f 
0.43f 
29.24+ 


2.05f 
Neutral 
50+ 03 


100% 
100% 
81% 


0.821 


4585.7/4198 


NGC 3227 Obs 4 


1.79f 


07+ 01 


4.8f 


1.2f 


Q1+ - 02 
u - a -0.02 


0.2f 


2.3f 


o 1O+0.22 
o - 1,3 -0.22 


3.37f 
0.43f 
29.24+ 


2.05f 
Neutral 
1 50+ - 12 

1.0U_o 13 


100% 
100% 
75% 


0.82f 




NGC 3227 Obs 5 


1.79f 


13+ 001 


4.8f 


1.2f 


0.79tr 2 


0.2f 


2.3f 


11 7O+0.48 
14- 'y_0.48 


3.37f 
0.43f 
29.24+ 


2.05f 
Neutral 

o-7o±8:8i 


100% 
100% 
89% 


0.821 




NGC 3227 Obs 6 


1.79f 


n in+ - 01 


4.8f 


1.2f 


0.64t°;«g 


0.2f 


2.3f 


q 4fi+0- 47 


3.37f 
0.43f 
29.24+ 


2.05f 
Neutral 
1 41+0+2 


100% 
100% 
89% 


0.821 




NGC 3516 Obs 1 


i-68±g;°i 


n i c;+o.oi 
u - 10 -o.oi 








1.0* 


6 -J+O- 3 
D - a -0.3 


q cn+0.22 
O.OU_q jo 


13 40+ 9 06 
2 7fi +0.30 
Z ' ,D -0.17 

o-i2±8:8? 


, cc+0.22 
°- oo -0.10 

1 93+ 09 
± - M ' J -0.06 

Neutral 


100% 
100% 
100% 
























24-83±i:« 


60+ 08 
u - DU -0.09 


83% 




1263.2/1126 


NGC 3516 Obs 2 


1.68t 


29+ 01 

u - Z9 -o.oi 


- 


- 


- 


iot 


6.3f 


1 qfi +0 - 12 


n or + 0.09 
u -°°-0.06 


1 02+ - 12 


100% 


0.82f 


















l-54l° i° 


9 9|-,+0.05 
z - zu -0.06 


100% 
























O.OS+Hf 
69.18±l|-f? 


Neutral 


100% 
























2.80±°i°o 


18% 






NGC 3783 Obs 1 


1 -Q + 0.01 

1 -' a -0.01 


i sn+ ' 01 

L - 3U -0.01 


74 +0.5 


1 1+ ' 1 
1 - 1 -0.1 


0.78tg : g 


o-9±8;i 




7.99i°; 2 « 


n -, a+0.02 

0.161,, 02 




100% 






NGC 3783 Obs 2 


1.79f 


-, 7 o+0.01 


7.4f 


lit 


2 ro+0.08 
z -°°-0.08 


0.9f 


4.4f 


8-26±g;Si 


, Q2+0.07 
°- uz -0.08 
n 9 7 +0.03 
u - zl -0.03 


2 (,r+0.01 
z.uo_ 01 
n oq+0.05 


100% 
100% 


0.99f 


2554.2/2305 














01 +0.08 
°- OJ --0.08 


2 nq+0.01 
z - uy -0.01 


100% 






NGC 4051 Obs 1 




0.l0l°;°i 


6.9tg; 2 


- --1-0 1 


u - lcs -0.01 


„ „_i_n i 

o-9±8:i 


9-6to: 5 


^ „, |_0 07 

^U'-O .06 


29+ 06 
u ' zy -0.05 

o-3i±g:8i 


9 Q7+0.05 
z -»'_0.04 
1 or+0.01 
t-OO-0.01 


100% 
100% 


„ — .-.-1-0 O r i 

n 70 ^ u - uo 
u - ' u -0.05 








fil+ - 01 






1 56+ 002 
l.oo_0_02 






1 K4+0.10 

J-^-o.os 




1 QK+0-03 
i.ao_ 03 


67% 






NGC 4051 Obs 2 


1.88f 


6.9f 


1.5f 


0.9f 


9.6f 


0.29f 


2.97f 


100% 


0.70f 


3204.7/2944 














0.31f 
8.96f 


1.85f 
1.95f 


100% 

11% 




NGC 4051 Obs 3 


1.88f 


n S9+ 01 


6.9f 


1.5f 


i nfi+ 01 

l.UD_ Q1 


0.9f 


9.6f 


1 40+ - 12 


0.29+ 
0.31f 
8.96f 


2.97f 
1.85f 
1.95f 


100% 
100% 
27% 


0.70f 





s 
o 

Co 

3 

S3 
^1 . 

o 
s 

Co 



3 
s 



s 

CO 

03 



Table B2 - continued Components for the baseline model (without broad or diskline emission) to the observations with Suzaku XIS, HXD and BAT data from Swift. a Unabsorbcd 
POWERLAW normalisation given in units (10 -2 phkcV -1 cm~ 2 s _1 ). 6 Flux for COMPTT quoted over the 0.6-10.0 keV range and REFLIONX over the 2.0-100.0keV range in units 
10~ 1:L crgcm~ 2 s _1 . c Ionization parameter given in units ergcms -1 . d Column density measured in units 10 22 cm~ 2 . e v turh = lOOOkms - 1 . ^ only the intrinsic powerlaw is 
absorbed. * indicates a frozen parameter. \ indicates parameters are tied during the analysis of multiple observations. 



Object 


Powerlaw 




COMPTT 






REFLIONX 






Warm absorber 










r 


Norm™ 


kT (keV) 


T 


Flux 6 






Flux 6 




log«) c 


Cfrac 


BAT/XIS 


xl 


NGC 4151 


1 56+ 001 


27+ lM)1 
u - zl -0.01 


o 1+0.1 
J - 4 -o.i 


io.6l8:l 


3 6S+0- 07 
o - oo -0.07 


1.0* 


43.6l^ 


12.63±?;S| 


09+ 002 

u - ua -0.01 

14.45+°" 

— U. lo 

45.82+ffl 


< 0.12 
2.16+°°; 
0-52l°?« 


100% 
100% 

58% 


-, ,-7+0.02 
101 -0.02 


647.1/583 


NGC 4593 


l-60±g;gi 


21+ - 01 
u - zi -0.01 


- 


- 


- 


1.0* 


5.4l^ 4 7 


61 + 3 ' 24 


rro+O.OO 


2 .i9i°;i; 


100% 


-, 07+O.12 
1-87_0.11 


293.62/289 


NGC 5506 Obs 1,2 




°-°°-0.02 


- 


- 


- 


5+ ' 1 


10-5±8;1 


9.69l -l 


q qo-1-0 . 55 
^—0.45 

3.26+°' 3 3 i 

151 + nno e 
1-0 —0.02 


< 0.05 
0.59+°'?,° 

1 30+2S? 


100% 
100% 
100% 


u - ' y -0.01 








5 S6+ ' 02 












8.48l ji 


449.0+ 4 °'°, 

— 36.0 


< 1.20 


19% 






NGC 5506 Obs 3 


2.08f 


_ 


_ 


_ 


0.5f 


10. 5f 


2.23f 
3.26f 
1.51 e f 
448. 76f 


0.05+, 
0.59f 
1.30f 
1.20+ 


100% 
100% 
100% 
23% 


0.79f 


NGC 5548 Obs 1,2, 


1 70+ - 01 
L - ' u -0.01 


u - 44 -o.oi 


- 


- 


- 


i n +0.l 
1 - u -0.1 


< 2.2 


9 19+1.27 
z - iz -0.75 


0.59in 9 3 


9 90+O.O7 
z - zo -0.08 


100% 


0.86l°:° 4 4 


1772.0/1734 


3,4,5,6,7 


















0.26l° :°oI 


n qs+ ' 14 
u.ye_ 08 


100% 






NGC 7213 


1 74+0.01 
L - ' 4 -0.01 


o fii+o- 01 


< 61.8 


9 1+0.5 


03+ ' 08 
u ' u,:5 -0.02 


1.0* 


29 fi+ 26 - 9 


n qo+0.63 
u -°°-0 .20 








o.6il°,:°J 


703.5/707 


NGC 7314 


J-o^-o.oi 


22+ ' 01 

u - zz -o.oi 


< 15.8 


< 0.8 


n S4+L03 

u - J *-0.26 








n _ . _i_n 02 


Neutral 


100% 




592.3/545 


IN Kj\j ( 4oy 


L -'°-0.10 


_ c-7+0.04 
0-57-0.06 


< 0.4 


°- 8 -0.5 


n 99+O.O3 
u - zz -0.03 


1 - b -0.3 


^ 1 1 n 
< 11. U 


1 RQ+ 014 
1.D9_q 2g 








1 nQ +0.23 
±.UO_q 23 


Q/i n /q 1 
o4U.Z/ olz 


PDS 456 


n 4 i+0.03 
z - 4i -0.03 


23 +(U)1 
u - z,3 -0.01 


< 33.4 


1.4l°;6 


< 0.204 








421.40 

o 09+1.90 

°- oz -1.09 


3 16+°' 19 
Neutral 


50% 
50%f 




146.9/132 


PC 1211+143 


1 S9+001 
i -° z _0.01 


n i2 ' 01 
u - lz -o.oi 


< 18.3 


est?,;? 


< 0.11 








3.05ir 6 5 7 


2 85+ ' 06 
Z - OJ -0 .08 


100% 




813.62/702 


RBS 1124 


i 71 +0.01 
1 -' 1 -0.01 


12+ 001 
u - lz _0.01 


< 52.9 


0.81J 1 , 3 


< 0.19 


7 +0.5 
u -'-0.3 


< 15.1 


S2+ 4 - 34 










386.8/364 


SWIFT 


o ii +0.03 
Zi±1 -0.02 


, ,- 9 +0.04 








n 5+ - 1 


< 13.0 


2.28l°: 3 I 


0.08l°° 2 


Neutral 


100% 


70+° - 10 
u - ' u -0.10 


871.6/867 


J2127.4+5654 




























TON S180 


2 i 1+0.02 

z - i4 -o.oi 


23+ - 01 

u - zo -o.oi 


6.3+4" 


i.2ti: 8 2 


< 1.14 


2 4+O.8 
Z - 4 -0.5 


269.6tgg;g 


in+ 03 

u - iu -0.03 


02+ 01 
u - uz -0.01 

> 415.20 


^^-l.oe 
1 qq+o-54 
i.yy_ 74 


100% 
60% 




750.0/692 
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Table B3. FeK region properties - distant emission lines and ionized absorption zones in the baseline model. We typically model the 
absorption zones with an XSTAR grid with i>t ur b = 1000 km s — 1 . a v t uib = 3000 km s~ 1 . Line flux given in units 10 — 5 erg cm~ 2 s _1 . 
Column density given in units 10 22 cm -2 , f indicates that parameters are tied in multiple observations. * denotes a frozen parameter. 



Object Emission lines High £ zone 

LineE EW (cV) Flux A\ 2 N H log£ v ou t (kms^ 1 ) A\ 2 



3C 111 


6 - 50 -0.05 
fi 7S+0- 11 
b -< 8 -0.08 


10+23 

19_n 

-11+8 


n t-q + 0.62 
°- 53 -0.30 
n ov+0.21 
0.27_ 21 


-9 

rr- 
-0 


> 3.60 


1 10+I + I 
4 - 4 °-0.37 


a/ifinn+ 4 2600 
84600 _ 42400 


-28 


3C 120 


r- .yo + 0.07 


o+4 
8-4 


0.42l°; 2 J 


-12 










oL/ 44D 


fi fiS+ - 07 

6 - 68 -o.io 


i 7 +402 


00+6.74 

°-2 9 _0.15 


fi 
-O 


\ 99 QC 


A Qo + °' 05 

•*.»J_ gQ 


4000_ 500 


-0 


A S~i T A fie 

4C 74.26 


fi fio+0.09 
6 - 68 -0.08 


m+7 
10_ 7 


0.431JU 


-5 










Ark 120 


o oo + 0.04 
6 - 66 -0.05 
fi qc + 0.03 
6 - 95 -0.03 


oo+7 

20I 7 

qi +9 


0.72^^25 
n no "I - " ■ - ' 


-20 
-31 










Ark 564 


6.63* 


00 + I6 
20Ijg 


0-3llu|o 


-7 










rairall y 


fi v1 +0.03 
6 - 71 -0.02 
fi Qfi+0 03 
6 - 96 -0.03 


20_ 6 
20_ 6 


0.58^q'J 7 

n /i H - 0.13 
0- 48 _0.15 


-32 

-20 










1C 4329A 


c riK+0 05 
6 - 95 -0.06 


n+5 
9 -5 


o-ystoii 


-9 










TTQ A C 1 QOO/1 Qonn 


fi 7o +0.06 
6 - 72 -0.06 


101 -67 


ofi+0. 04 
0-06_ 04 


-6 










MGG— 02-14-0U9 


fi 94+O.05 


qn+24 


0+0.10 
°- 16 -o.io 


-7 










MOCj— o-zo-lu 


fi /ln+0-01 


54_ 6 


6 - 02 to:68 


-82 










MCG-6-30-15 


6.97* 


10l 4 4 




-7 


3-99t?;i 


O Q4+0-08 

d - y4 -0.25 


3200t 4 ~ 


-298 


MGG+o-ll-ll 


o oq+0.03 
6 - 93 -0.03 


01 +6 
21_ 6 


! oo+0. 40 

1.39_ ()40 


-17 










T\/1"R OO^I 1 *7fi 
Ivlrx ZZOl— 1 i o 


fi 40+O.O6 
b - 48 -0.06 


in+9 
12_g 


fic+0.47 
— 50 


1 
-4 


oq+0.31 


10+O.23 

•^-i^-o.os 


qq finn +7000 


qo 
-ay 


Mrk 79 


fi qo+0.01 

6 - 39 -o.oi 


llOljg 


or~ + 0.^8 

2.05+0.28 


-45 












fi fin+0 05 


Ofi+14 


K1 +0.27 
°- 51 -0.25 


-11 












c no+0-06 


on+17 
20 -17 


0.30l ;^ 


-5 










TV/Tv.'U 1 1 A 

ivirk iiu 


fi KK +0.12 

b - bb -o.ii 


11+IO 

n -io 


0.26_ .24 


- 1 










TV T..1 nor 

Mrk ooo 


fi KQ+0.03 
6 - 68 -0.05 


/fn+8 
40 -8 


o,+0.13 
°- 64 -0.13 


-67 












fi Ofi+0-06 
6 - 96 -0.16 

o rr.+O 05 


01+9 
21_g 


00+O.12 
U - 2 °-0.12 


1 A 

-14 










Mrk 359 


00 + I8 

32119 


-1 o+O 11 


-6 










Mrk 509 


o ,o+0.03 

6.43+ 02 

« 07+O.O8 


00+17 
23_ii 
a+12 

8 -6 


1 on + 0.98 
i^S+O.Bl 
, K +0.71 

0.45_ 037 


-9 

-4 










Mrk 7oo 


r> ,.,+0.04 
6 - 64 -0.05 


27 -10 


07+0.14 
0.37_ 14 


-28 


> 3.90 


c- ,,+0.41 
5 ' 44 -1.16 


kooo+2900 
ozuu_ 2900 


-108 


Mrk 766 


6.64f 


27f 


0.37f 






3.68+^ 


5200f 






fi RQ+O-06 


iq+12 
i9 -13 


00+0.20 
U.32_ Q 22 


-0 










NGC 1365 


6 64+ 001 


26l 3 


71+ 014 

u - ' 1 -0.15 


-11 


1 1 o 7 +5.93 a 

±±.U(_ 2 .i g 


74+0.07 
6 - '+-0.06 


< 300 


-2052 




fi flO+0-01 

6.89_ 01 


28_ 5 


,1+0.09 
°- 41 -0.08 


-32 


co oi + 7.43 a 
52 - 84 -13.95 


„ Qo+0.02 

3.96_ 004 


/iann+ 4 00 
4900_ 500 




NGC 1365 


6.64f 


24_jj 


0.71f 




41 co+17.64 a 


4 00+ ' 21 

^■ uu -0.03 


< 300f 






QO.-I- 

D.oyj 


41 -8 


fV y<1 + 

U.41J 




' 70 a 
< D. / 2 


ifi+0.02 
3- 4 O-0.02 


< OUU 




NGC 1365 


6.64f 


3 7+ J 


0.71| 


- 


< 190.40° 


> 4.11 


< 300| 






fi Qn+ 
6.oyy 


A9 + 9 
42 -8 


n a 1 + 
U.41J 




c QO+0.49 a 

5-98_ 48 


11+O.O2 
3- 41 -0.02 


< 2UU 




NGC 2992 


6 40+ - 01 


154+I8 


2 60+ ' 30 
z ' DU -0.30 


-23 












R Q7 +0.05 
b - 97 -0.06 


84 -54 


-1 o+O. 07 
u - ±,:, -0.08 


fi 
-0 










NGC 3227 


6 40+ 01 
°-* u -o.oi 

a oo+O.Ol 
6-»3_ 01 


591! 
28+ -3 


2 98+ ' 17 
• i - ao -0.19 
1 on+0.13 

i^OIo.ia 


-80 
-25 


20.881 2 ; 4 * 


4 35+O.O6 
^■ o,J -0.08 


< 2100 


-67 


JNGC 6221 


6.40f 

fi 8*3-1- 
D.OO | 


1 07+0 

107 7 

5fi+ 6 
t>0_6 


2.98f 




20.88f 


4.35f 


< 2100f 




NGC 3227 


6.40f 
6.83| 


75l 4 

07+4 
37_ 4 


2.98f 
1.20t 


- 


20.88f 


4.35f 


< 2100f 




Tvrpp Q997 


fi zin-f 

0.4U | 

6.83| 


16S+ 9 

»o_io 


9 QS-f 

z.yc | 
1.20t 




on QQ+ 




/ 91 nn+ 
ziuu i 




NGC 3227 


6.40f 
6.83f 


85+5 
42l« 


2.98f 
1.20f 




20.88f 


4.35f 


< 2100f 




NGC 3227 


6.40f 
6.83f 


W9 +6 7 
56l D 


2.98f 
1.20f 




20.88f 


4.35f 


< 2100f 




NGC 3516 


6 42+ 01 
D -^ z -o.oi 


5ll 4 4 


2 19+ - 18 
z - la -0.18 


-15 


10+1. 24 
^• iZ -0.95 


07+0. 18 
3- 8 '-0.14 


< 4500 


-23 




6 69+ - 03 


8 t^ 


39+ 014 
u - oa -0.13 


-8 


qq+2.79 
Z -J'J_0.82 


3- 8 0to:?2 






NGC 3516 


6.42f 
6.69f 


108+8 


2.19f 
0.9f 




< 9600 




NGC 3783 


6 96+ 01 

D - MD -0.02 


30+^ 


07+0. 58 
-0.64 


-21 


5 46+ 21 ' 31 
°'^°-2.02 


^ 00+0.66 
^' zz -0.19 


< 400 


-50 


NGC 3783 


6.96f 


26^7 


2.37f 




2 40+ - 75 
^•^ U -0.58 


3-87±8:ll 


< 400f 
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Table B3 - continued FeK region properties - distant emission lines and ionized absorption zones in the baseline model. We typ- 
ically model the absorption zones with an XSTAR grid with v tU rb = 1000 km s . a Dturb = 3000 km s — . Line flux given in units 
10 -5 ergcm~ 2 s _1 . Column density given in units 10 22 cm -2 , t indicates that parameters are tied in multiple observations. * denotes a 
frozen parameter. 



Object 




Emission lines 




High £ zone 










LincE 


EW (eV) 


Flux 


Ay 2 


N H 


losrf 


v ou t (km s 1 ) 


Ay 2 

*— 1 A 


NGC 4051 


fi 4-J+0.03 


42+ 2 
21±g 


o.53ir i 

0.26l°;H 


-5 
-3 


1 1 7 +0.27 
± - 1 ' -0.26 


3.68l u; ° 5 5 


5800+ 1400 


-51 


NGC 4051 


6.43f 
6.62f 


19+i 
10+ 4 


0.53f 
0.26f 


-5 
-3 


1.17t 


3.68f 


5800t 




NGC 4051 


6.43f 
6.62f 


25+1 

— 1 

12+| 


0.53f 
0.26f 


-5 
-3 


1.17f 


3.68f 


5800+ 




NGC 4151 


6.381°;°! 




8 79+ - 56 
°- ' a -0.56 


-193 


9 ro+41.80 

z.oo_ 2 .53 
145 97+66-56 
i40 - z ' -107.01 


> 3.57 

> 4.41 


128001^ 
12800f 


-164 


NGC 4593 


6.42l° °l 

6 71 +° 13 
0- < ±-0.09 


173+21 

20+"! 
— 10 


9 90+O.27 
—0.28 

0.29+°H 

—0.14 


-24 
-8 










NGC 5506 


6.63* 
6.98+°,° 


< 25 


< 1.10 
0-88+0 68 

— U.DO 


-15 
-36 


- 


- 


- 




NGC 5506 


6.63* 
6.98f 


< 25f 
10f 


< l.lOf 
0.88f 


- 


- 


- 


- 




NGC 5548 








- 


0.88l^ 4 4 


yo+0.23 

°- 1 °-0.29 


< 1700 


-12 


NGC 7213 


fi oq+0.01 

6 61+ 004 
fi Ofi+013 


67 +io 
D '-n 

< 27 

ol -20 


, -0+0.27 

l-'»_0.28 

< 0.75 
71+ - 34 


-45 
-31 

-Ol 










NGC 7314 


fi qo+0.02 
fi QS+0- 05 




1 ic+0.16 
i - ±o -0.16 


-131 












3811? 


33+0.15 


-12 










PDS 456 


6 7^+ 06 
D -'°-0.07 


25+ 17 
— 19 


12+ 008 
U-J - -0.09 


-5 


> 6.31 


4 Q8+ ' 35 


82200^°° 


-27 




7 00+ ' 48 


1Q+19 
ia -17 


09+ 09 
u - ua -0.08 


-4 


22+ - 19 
u - zz _0.14 


3 03+° 14 


< 163200 




PG 1211+143 


6 40+° 10 
°' 4u -0.09 

6 7^+ - 07 
,o -0.07 


nn+16 
zz -16 

4412* 


14+ - 10 
u - ±4 -o.io 

24+ ' 18 
u - Z4 -o.n 


-4 
-3 


- ,o+14.22 
'• ±,3 -3.04 


70+0.29 

°- ' a -0.09 


ISSOOl 4 ^ 


-30 


SWIFT J2127.4+5654 


^ + -°o°ol 
6.98lg;» 


25111 
1611^ 


94+ ' 41 
u - a4 -0.41 
49+0.40 


-16 
-4 
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Table B4. Components for the dual reflector fit to the observations with Suzaku XIS, HXD and BAT data from Swift. COMPTT and 
warm absorber parameters are consistent with those in the baseline model. Here we quote the inner blurred reflector properties, some of 
which are tied to the out unblurred REFLIONX. Fits include partial covering geometries where required as per the baseline model. Where 
a reasonable fit can also be obtained without the use of a partial coverer, both scenarios are tabulated, for those without partial covering 
we quote the change in \ 2 with respect to the dual reflector plus partial covering fit. Note that in some objects accretion disc parameters 
cannot be constrained, denoted by -. a Ionization parameter given in units ergems -1 .* Denotes a frozen parameter, f A\ 2 in relation to 
the dual reflector fit without partial covering (which for these objects has not been tabulated) compared to the fit with partial covering 
tabulated here. b Best-fit parameters to the 2009 NGC 3783 Suzaku data only for the Solar abundance model presented in Section 4.3.4. 



Object 


r 






q 


a 


i° 


Ax 2 


xl 


3C lfl 


1 - DU -0.03 


0-8181 


< 102 


< 2.69 


- 


- 




1087.8/1092 


3C 120 


1 62+ ' 01 


1.0* 


1911 


-1 7+0.3 

1 -'-0.4 


- 


18l 2 




3564.8/3448 


3C 382 


1 7Q+0-03 
L - ' y -0.02 


1.0* 


< 1.4 


< 3 


— 


3018 1 




972.3/933 


3C 390.3 


1 70+O.OI 
L - ,z -0.01 


1.4 


< 2 


i-sli:" 


— 


39l?° 




1478.7/1482 


Ark 120 


1 qs+004 


in 1 
Q rj"|-U.± 

z - z -0.6 


< 21 


j~n 1 

9 (ITU. ± 

z - z -0.6 


— 






715.5/644 


Fairall 9 


1 Qi+0.01 
1 - yi -0.01 


i-7l8:l 


311 


7+O.6 
Z -'-0.5 


< 0.95 


42+1 




3558.3/3271 


IC 4329A 


1 qo+0.02 
^^-O.Ol 


o.8i8;l 


7t\ 


2.4* 


< 0.73 


36* 




2342.8/2199 


MCG-02- 14-009 


1.901°;°! 


n 7+O.5 
u -'-0.3 


< 13 


< 2.0 


- 


> 29 




601.5/539 


MCG-05-23-16 


1.84±°;°i 


1.0* 


< 1.0 


1 fi+°- 5 


- 


24* 




1989.5/1890 


MCG-06-30-15 (p/c) 


2.05i°;°i 


1.0* 


< 11 


9 0+O.2 
z - o -0.1 


o-6il8:i? 


35l 2 




2061.4/1823 


MCG-06-30-15 (no p/c) 


2-0918:81 


1.0* 


< 10 


q+°- 2 

z - y -o.i 


- 


3611 


-41 


2020.2/1823 


MCG+8-11-11 


1 00 +0.01 

1 -82_ .02 


1.0* 


7+3 
'-2 


1.9* 


— 


— 




1034.0/933 


MR 2251-178 


1 r.A+0-06 
1 -°*-0.10 


1.0* 


- 


3.0* 


- 


— 




969.0/897 


Mrk 79 


1-6118:82 


1.0* 


9841^8 


3.0* 


< -0.25 


< 25 




573.2/542 


Mrk 335 


2-04l°;°2 


9 rr + 1.1 

z - J -o.6 


25± 7 3 


9 n +0.2 
z ' u -0.3 


- 


50l? 3 




820.0/721 


Mrk 509 


1 fiq+ ' 03 


1.0* 


23±l 


1 K+0.8 
L -°-1.3 


- 


35* 




1959.4/1867 


Mrk 766 (p/c) 


1 qq+ ' 01 


1.0* 


311 


- 


- 


- 




1055.1/995 


Mrk 766 (no p/c) 


2 01 + 001 
z - ul _0 .01 


1.0* 


411 


3.0* 




80l 2 


+54 


1109.6/1002 


Mrk 841 


1 96+°'° 5 
1 - au -0.04 


1 0+°' 3 


< 2 


9 7 +0.3 
z - ' -0.2 


> —0.40 






91 8/853 

ZS -L\J .O / O * J -J 


NGC 1365 


1 69+°' 01 


1.0* 


25+1 


3.0* 


> 0.70 


8311 


+7995t 


2157.0/1974 


NGC 2992 


1 60+°'° 4 


1.0* 


59111 


< 1.8 








1082.3/1076 


NGC 3227 


1 sn+ - 01 
1-80_ 01 


0-418:1 


611 


9 c+0.3 
z -°-0.3 


< -0.35 


47l 3 


+433t 


4373.9/4075 


NGC 3516 


1 7n+ - 01 

L - ' u -0.01 


1.0* 


611 


2 6+° 1 
z -°-0.3 


< -0.50 




+814f 


1228.6/1121 


NGC 3783 


i.79i8:8i 


1 1+ - 1 

i - 1 -o.i 


411 


n 0+14 
J - a -0.6 


< -0.35 


19l 7 




2502.6/2300 


NGC 3783 6 


i.84i° o: 8? 


1.0 ±0.2 


< 11 


n Q+0.5 


< 0.45 


< 13 




1413.7/1374 


NGC 4051 (p/c) 


l-89_ 01 


i-il8:? 


18+3 


3.0* 








3177.1/2939 


NGC 4051 (no p/c) 


1 93+ 001 


0-818:1 


12+1 


6 l+°' 3 
D - 1 -o.i 


> 0.99 


13l^ 


+ 137 


3321.9/2943 


NGC 5506 (p/c) 


2.081°;°^ 


0.7I8.I 


10+I 


3.0* 




26l^ 




3191.2/2868 


NGC 5506 (no p/c) 


2.o7±8:8i 


0-818:1 


911 


1 y+0.3 
1 -'-0.8 




48lf 


-12 


3179.5/2871 


NGC 7469 


i.sollllll 


u.y_ 2 


< 13 


1 7+0.1 
l -'-o.s 


72+ ' 18 
u - ' z -0.17 


80+ 5 




780.9/805 


SWIFT J2127.4+5654 


2-i9l8:8I 


n q+ - 3 


< 18 


9 9+0.4 
z - z -0.9 




42ll 6 




830.8/866 
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Table B5. List of 15-50 keV component fluxes for each object and the full model 2-10 keV flux. The full model flux includes all model 
components, reflector flux is the the 15-50 keV flux of the REFLIONX component in the HXD whereas the continuum flux is the model 
flux minus any contribution from reflection or partial covering. The hardness ratio is the ratio of the full model 15-50 keV to 2-10 keV 
flux, a powerlaw with slope T = 2.0 gives a hardness ratio of 0.75. The reflection fraction, _Ri5_5o, is the ratio of the 15-50 keV reflector 
flux to the 15-50 keV continuum flux i.e. to the full model minus reflector flux, Equation[T] Flux given in units 10 — 11 ergcm — 2 s — 1 . 



Object 


Full model 


Reflector 


Continuum 


Full model 2-10 keV flux 


Hardness ratio 




i?15-50 


1H 0419-577 Obs 1 


2 


.848 ± 0. 


110 


0. 


237+°,' 


138 
092 


2.156 ±0.055 


1 


.753 


1. 


625 


± 


0. 


.063 





091 


+ 
-0 


053 
08,0 


1H 0419-577 Obs 2 


2 


.315 ± 0. 


.083 


0. 


301+?, 

OOX_g 


187 
114 


1.719 it 0.055 


1 


.374 


1. 


.685 


± 





.060 


0. 


149 


+0. 
-0 


094 

.058 


3C 111 


3 


.292 ± 0. 


.097 




o 




3.292 ±0.097 


1 


.947 


1. 


.691 


± 





.050 










3C 120 Obs 1 


U 


zou it u 


zDy 


n 
u. 


ooy q 


139 
130 


5.347 ±0.055 




.uzy 


1 


*347 


+ 


n 
u 


0^9 


0. 


ioo 


+0. 
-0 


027 
1)20 


3C 120 Obs 2 


5 


.788 ± 0. 


.187 


1 


082+0 


102 
13 9 


4.706 ± 0.055 


3 


.963 


I 


.461 


± 


o 


.047 





.230 


+0. 
-0 


030 
032 


3C 382 


o 


1 o& 4- n 

. 1UO it u 


1UD 


1 


flQQ+0 

uyo q 


092 
000 


4.015 ±0.035 


4 


04^ 
. U40 


1 


zoo 


+ 


n 


uzo 


0. 


272 


+0. 

-0 


025 
157 


3C 390.3 


K 
O 


.oiz it u. 


182 


1 


AAA + 1 

■ 444 -o 


487 
300 


4.368+^ 


Q 
.> 


101 




P,7 1 


+ 


n 

u. 


uoy 


0. 


DDI 


+ 
-0 


359 
090 


3C 445 


2 


.692 ± 


.110 


0. 


qn9+ 3 


278 

333 


1.790 ±0.597 





.699 


3. 


.851 


± 





.157 


0. 


.504 


+2. 
-0 


051 
211 


4C 74.26 


Q 

o. 


yyu it u. 




n 
u. 


821+°' 


541 


3.175 ±0.030 


.> 


1 97 


1 


978 


+ 


n 

u. 


n^n 
uou 


0. 


zoy 


+0. 
-0 


295 
170 


Ark 120 


Q 

o 


.ooo ± U. 


i no 

.1UZ 


i 


.188_ (| 


502 
077 


2.670 ±0.169 


Q 

o. 


.U01 


1. 


ZOO 


1 

± 


n 
U 


.Uoo 


0. 


.440 


+0. 

-0 


231 
278 


Ark 564 


1 


.536 ± 0. 


.146 


0. 


,078i;; 


124 
050 


0.822 ±0.007 


1 


.837 





.836 


± 





.079 


0. 


.053 


+ 
-0 


085 
041 


Fairall 9 Obs 1 


3 


.370 ± 0. 


.076 


1. 


.269^0 


001 
08 j 


2.101 ±0.031 


2 


.317 


1. 


.454 


± 





.033 


0. 


.604 


+0. 

-0 


055 
052 


Fairall 9 Obs 2 


3 


.511 ± 0. 


.180 


1. 


043to 


070 
070 


2.468 ± 0.034 


2 


.165 


1. 


.622 


± 


0. 


.083 





.423 


+ 
-0 


0-14 
044 


IC 4329A 


ie 


i.778 ± C 


|.139 


5. 


.309^0 


104 

760 


11.469 ±0.033 


If. 


1.711 


1. 


.566 


± 





.013 


0. 


.463 


+0. 
-0 


011 
073 


IRAS 13224-3809 




rt9i 4- n 
.UZ1 it U. 


Ull 




u 




< 0.009 


n 


.UOO 


n 
U 


QC9 


1 

± 


n 


9nn 










MCG-02-14-009 





.789 ± 0. 


.093 


0. 


413_ 


001 
171 


u.o/o_ 094 


tJ 


.430 


1. 


.835 


± 


tJ 


.216 


1. 


.098 


+0. 
-0 


451 
728 


MCG-02-58-22 


8. 


.726 ± 0. 


.107 


1. 


"7K9+0 

.752_ 


301 
313 


6.974 ± 0.074 


1 


.872 


1. 


.791 


± 





.022 


0. 


.251 


+0. 
-0 


045 
040 


MCG-05-23-16 


14 


:.518 ± C 


'. 139 


1. 


■797t°' 


137 
321 


12.721 ±0.133 


8 


.929 


1. 


.626 


± 


tJ 


.016 


0. 


.141 


+ 
-0 


011 
020 


MCG-06-30-15 


5 


.013 ± 0. 


.067 


0. 


■584_ 


050 
100 




1 


.160 


1. 


.205 


± 





.016 





.132 


+0.012 
0.036 


MCG+8-11-11 


1C 


i A A n 1 n 

1.440 ± L 


1.160 


3 


1,1+1. 

111_1 


507 
020 


7.329 ± 0.039 


6 


.480 


1. 


.611 


± 


0. 


.025 


0. 


.424 


+ 
-0 


224 
151 


MR 2251-178 


u. 


a&k 4- n 
.400 It U. 


. lUo 




n 
U 




5.365 ±0.108 


A 
4 


.ZoU 


1 

1 


909 

zyz 


1 


n 


.UZO 










Mrk 79 




.875 it 0. 


.139 




o 




1.875 ±0.139 


\ 


.466 


I 


.279 


+ 


o 


.095 










Mrk 110 


2 


.990 it 0. 


.113 





267+2 


559 
083 


2.723 ±0.052 


2 


.126 


I 


.406 


+ 





.053 


0. 


.098 


+0. 
-0 


206 
031 


Mrk 205 


1 


.184 it 0. 


.062 


0. 


299+°, 


040 
113 


n 71 1 +0.023 
u - ' li -0.046 





.934 


I 


.268 


+ 





.066 


0. 


.338 


+0. 
-0 


05:', 
137 


Mrk 279 


2 


.322 it 0. 


.101 


0. 


736+°, 


187 
104 


0.554 ±0.050 





.489 


4. 


.748 


+ 





.207 





.464 


+ 
-0 


183 
138 


Mrk 335 


1 


.502 it 0. 


.125 


0. 


.288 + q 


088 
088 


i 9 i 1+0.024 
i -^ 14 -0.048 


1 


.490 


1. 


.008 


± 





.084 


0. 


.237 


+0. 
-0 


078 
078 


Mrk 359 





.984 ± 0. 


.141 





337+°, 

DDI _Q 


285 
110 


0.647 ±0.054 





.515 


1. 


.912 


± 





.274 





.521 


+ 
-0 


509 
223 


Mrk 509 


7. 


.186 it 0. 


.162 


0. 


966+q 


077 
001 


6.220 ±0.070 


1 


.723 




.521 


+ 





.034 


0. 


.161 


+0. 
-0 


989 
098 


Mrk 766 Obs 1 


1 


.476 ± 0. 


.118 


0. 


329+9, 
Dzy_ 


072 
050 


1.140 ±0.023 


1 


.325 




.181 


+ 


o 


.094 


0. 


.287 


+ 
-0 


072 
061 


Mrk 766 Obs 2 


\ 


846 it 


154 


g 


414+0 
4i4 -0 


001 
078 





\ 


.363 


1 


354 


+ 


o 


113 


0. 


289 


+0. 
-0 


073 
005 


Mrk 841 


2 


.618 ±0. 


.105 





.948+?, 


570 

.338 


1.670 ±0.046 


1 


.416 


1. 


.849 


± 





.074 





.568 


+ 
-0 


394 
285 


NGC 1365 Obs 1 


4. 


.880 it 0. 


.116 


0. 


,752i;; 


203 
050 


0.162 ±0.054 


1 


.284 


6. 


.489 


± 





.154 


0. 


.182 


+0. 
-0 


005 
015 


NGC 1365 Obs 2 


■1 


.428 it 0. 


.127 


0. 


.860+." 


070 
003 


0.143 ± 0.073 





.610 


7 


.259 


± 





.208 


0. 


.241 


+ 
-0 


022 
020 


NGC 1365 Obs 3 


3 


.615 it 0. 


.060 


0. 


909t° 


088 
051 


0.146 ±0.075 





.385 


9. 


390 


± 





.116 


0. 


336 


+0. 
-0 


017 
021 


NGC 2992 


2 


.106 it 0. 


.101 





620t° 


113 
113 


1.486 ±0.057 


1 


.182 


1. 


.782 


± 





.085 


0. 


.417 


+ 
-0 


087 
087 


NGC 3147 





.294 ±0. 


033 


0. 


.072+" 


1)98 
098 


0.223 ±0.056 





.165 


1. 


.782 


± 





.200 


0. 


.324 


+0. 
-0 


407 
049 


NGC 3227 Obs 1 


7. 


.564 ±0. 


.156 


3 


oio+;i 


107 
109 


o 2nq +0.043 


3 


.956 


2. 


.357 


± 





.049 


0. 


661 


+ 
-0 


o-lo 
080 


NGC 3227 Obs 2 


6 


.705 it 0. 


.164 


5. 


.218+0 


147 
144 


365+ 011 


1 


.853 


3. 


.618 


± 





.089 


3. 


.509 


+0. 
-0 


529 
524 


NGC 3227 Obs 3 


6 


.858 it 0. 


.170 


3 


,119+;; 


153 
152 


0.682 ±0.016 


2 


.551 


2. 


.688 


± 





.067 


0. 


.994 


+ 
-0 


079 
079 


NGC 3227 Obs 4 


4. 


.704 ±0. 


.170 


1. 


.127+;; 


113 
111 


0.228 ±0.013 





.996 


1. 


.723 


± 





.171 


7. 


.153 


+2. 
-2 


538 
524 


NGC 3227 Obs 5 


6 


.194 ±0. 


.179 


3 


.128+;; 


129 
127 


0.360 ±0.014 


2 


.138 


2. 


.897 


± 





.084 


1. 


.020 


+ 
-0 


085 
084 


NGC 3227 Obs 6 


5 


.283 it 0. 


.193 


3. 


.250+0 


100 
158 


n 991 +0.011 


1 


.571 


3. 


363 


± 





.123 


1. 


.599 


+ 0. 
-0 


212 
211 
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Table B5 - continued List of 15-50 kcV component fluxes for each object and the full model 2-10 keV flux. The full model flux includes 
all model components, reflector flux is the the 15-50 keV flux of the REFLIONX component in the HXD whereas the continuum flux is 
the model flux minus any contribution from reflection or partial covering. The hardness ratio is the ratio of the full model 15-50 keV to 
2-10 kcV flux, a powerlaw with slope T = 2.0 gives a hardness ratio of 0.75. The reflection fraction, ili5_5o, is the ratio of the 15-50 kcV 
reflector flux to the 15-50 keV continuum flux i.e. to the full model minus reflector flux. Flux given in units 10 _11 ergcm~ 2 s _1 . 



Object 


Full model 


Reflector 


Continuum 


Full model 2-10 kcV flux 


Hardness ratio 




^15-50 


NGC 3516 Obs 1 


7.321 ±0.117 


2 


nqq+° 


132 
.120 


0.920 ±0.061 


2.370 


3 


.089 


-4- 


o 


.049 


0. 


400+° 


029 
026 


NGC 3516 Obs 2 


3.340 ± 0.057 


1 


.198+?, 


073 
073 


0.920 ±0.061 


2.370 


1 


.409 


± 





.024 


0. 


.559+° 


042 
042 


NGC 3783 Obs 1 


9.917 ±0.178 


A 
1 


.460_ () 


.145 
145 


5.452 ± 0.042 


4.586 


Z 




-L- 


n 
u 


uoy 





81Q+ 

»iy_o 


044 
026 


NGC 3783 Obs 2 


12.003 ±0.100 


4. 


650±g 


152 
146 


7.353 ± 0.043 


5.921 


2 


.027 


± 





.017 





.632+° 


026 
025 


NGC 4051 Obs 1 


1.648 ±0.089 





.o92_q 


.039 
.033 


0.339 ± 0.034 


0.873 


1 


.888 


± 





. 102 





t; C1 +0 

■561_ 


063 
059 


NGC 4051 Obs 2 


3.223 ±0.052 





.8871°, 


058 
.046 


2.116 ±0.035 


2.464 


1 


.308 


± 





.021 


0. 


3801q 


028 
023 


NGC 4051 Obs 3 


2.591 ±0.110 


U. 


QOQ + 


.071 
076 


1.322 ±0.034 


1.794 


1 


AAA 

.444 


± 


U 


.U01 


0. 


465_ 


033 
036 


NGC 4151 


15.901 ±0.109 


■7 
( . 


oz0_ ( , 


865 
.917 


4.512 ±0.167 


4.352 


Q 
O 


A 

.004 


i 

± 


U. 


UzO 





yoy _q 


322 
.159 


NGC 4593 


2.152 ±0.130 





■ 388+.Q 


061 


1.764 ±0.084 


1.041 


2 


.067 


± 





.124 





.220+J 


196 


NGC 5506 Obs 1 


17.800 ±0.153 


6 


.279 +0 


'227 
220 


9.755 ±0.036 


10.379 


1 


.715 


± 





.015 





5451° 


024 
023 


NGC 5506 Obs 2 


17.114 ± 0.210 


5 


.952+° 


.218 
218 


9.375 ± 0.035 


9.890 


1 


.730 


± 





.021 





.533+° 


024 
024 


NGC 5548 


3.998 ± 0.087 


1 


.364+.° 


.817 
483 


2.624 ±0.060 


1.841 


2 


.172 


± 





.047 


0. 


■5181° 


350 
207 


NGC 7213 


3.477 ±0.109 





• 208+.°, 


397 
126 


3.269 ± 0.054 


2.410 


1 


.443 


± 





.045 





0641° 


222 
039 


NGC 7314 


1.427 ±0.110 









1.427 ±0.110 


0.879 


1 


.623 


± 





.125 








NGC 7469 


3.478 ±0.102 





.923^° 


076 
158 


2 TO +0.179 


2.102 


1 


.655 


± 





.049 





.3611° 


035 
007 


PDS 456 


0.249 ± 0.062 









0.130 ±0.006 


0.353 





.795 


± 





.176 








PG 1211+143 


0.498 ± 0.249 









0.498 ± 0.249 


0.391 


1 


.274 


± 





.637 








RBS 1124 


0.905 ±0.057 





■2H±o 


862 
138 


0.694 ±0.058 


0.494 


1 


.832 


± 





.115 





304+4 


310 
209 


SWIFT J2127.4+5654 


3.533 ± 0.095 


1 


.276+° 


207 
621 


9 OK7+0.059 
-2-^0 * -0.074 


3.348 


1 


.055 


± 





.028 


0. 


5651° 


108 
317 


TON S180 


0.773 ±0.129 





.052+° 


016 
147 


0.328 ±0.014 


0.564 


1 


.371 


± 





.229 


0. 


.0721° 


026 
072 
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